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in  the  tropical  ambient  and  about  tropical  cyclones,  models  for  the  cyclone-scale 
structure  and  maintenance  of  the  quasisteady  mature  stage  of  a  hurricane  are  presented 
and  evaluated.  Principal  attention  is  devoted  to  the  model  developed  by  Riehl  and 
Malkus  in  the  late  1 950 ' s ,  and  to  the  model  developed  by  Carrier  in  the  early  1970's. 

In  order  to  explain  hurricane  maintenance,  Riehl  and  Malkus  find  it  necessary  to 
postulate  the  existence  of  a  latent  and  sensible  heat  transfer  (from  the  underlying 
sea  to  the  air  in  the  frictional  inflow  layer  within  about  100  km  of  the  storm  center) 
that  is  augmented  by  about  an  order  of  magnitude  over  the  ambient-level  enthalpy  trans 
fer  in  the  autumnal  tropics.  Because  Carrier  permits  some  outward  sloping  of  the  eye- 
eyewal 1  interface  (in  contrast  with  Riehl  and  Malkus),  he  finds  that  hurricane-strength 
winds  can  be  explained  without  an  oceanic  heat  source  postulate;  Carrier  believes  th it 
oceanic  enthalpy  transfer  is  increased  in  a  hurricane  only  modestly  over  ambient  levels 
Further,  Carrier  emphasizes  the  role  of  the  large  mass  of  warm  moist  air  advected 
with  the  storm  from  its  inception  in  understanding  hurricane  lifetimes.  Analyses 
relating  to  maximum  swirl  speed  estimates,  solutions  describing  bound ary- layer  dynamics 
and  energetics,  and  interpretations  of  temperature  measurements  are  criticized  to 
evaluate  the  relative  merits  of  the  two  models  and  their  implications  for  successful 
storm  modification  by  current  seeding  practices.  Attention  is  then  called  to  unresolv 
questions  concerning  tropical -cyclone  intensification;  further  progress  on  the  deepeni 
of  tropical  depressions  will  probably  entail  more  quantitative  description  of  the  pro¬ 
cesses  that  maintain  convective  instability  in  the  tropical  ambient,  along  with  better 
understanding  of  the  spatial  concentration  of  cumulus  convection  by  a  CISK-like  mechani 
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1.  INTRODUCTION 


1 . 1  The  Dangers  and  Benefits  of  Hurricanes 

In  an  average  year  the  Atlantic  and  Gulf  coastal  states  of  America 
suffer  over  $100  million  damage  and  50  -  100  fatalities  owing  to  hurricanes; 
in  a  severe  year  damage  will  exceed  $1  billion  (Meyer,  1971a).  Hurricane 
Agnes  from  June  19-26,  1972  caused  damage  from  Florida  to  Maine  that 
has  been  estimated  at  $3  billion  (Anonymous,  1972),  and  Hurricane  Camille 
in  August  1969  cost  258  lives  (White,  1972).  The  threat  consists  of  high 
winds  [over  200  mph  is  known  (Riehl,  1972a)];  in  rainfall  [up  to  27  inches 
in  24  hours  (Schwarz,  1970)  --  since  1886  hurricanes  have  caused  over  60 
floods  in  the  United  States  (Alaka,  1968),  and  Hurricane  Agnes  is  estimated 
to  have  rained  in  toto  28.1  trillion  gallons  (Anonymous,  1972)];  in  storm 
surges,  especially  along  concavely-curving  coasts  [coastal  ocean  levels 
can  rise  15-20  feet,  with  Hurricane  Camille  setting  the  American  record 
at  24  Feet  along  the  Gulf  coast  (White,  1972)];  and  in  ocean  wave  heiahts 
[wave  heights  of  45  feet  are  known  (Riehl,  1954,  p.  298)].  On  the  positive 
side,  on  a  long-term  basis  the  hurricane-associated  rainfall  over  the 
Eastern  seaboard  states  is  about  one-third  of  the  total  annual  rainfall 
(Penner,  1972);  without  hurricanes,  droughts  seem  inevitable  on  the  west 
coast  of  Mexico,  in  southeast  Asia,  and  elsewhere. 

Globally,  there  are  about  eighty  tropical  storms  annually,  of  which 
about  fifty  intensify  into  hurricanes;  hurricanes  occur  in  all  oceans 
except  the  South  Atlantic  (Atkinson,  1971).  Under  the  current  r*ate  of 
the  art  for  path  predictions  for  either  military  or  civilian  use,  three 
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times  the  area  actually  hit  by  a  hurricane  is  typically  placed  under 
hurricane  warnings  (Meyer,  1971a;  Malone  and  Leimer,  1971).  In  the  greater 
Miami  area,  current  estimates  suggest  that  at  least  $2  million  are  spent  on 
preparations  whenever  hurricane  warnings  are  issued  (Simpson,  1971). 
Unnecessary  preparation  by  the  United  States  Department  of  Defense  instal¬ 
lations  owing  to  false  hurricane  and  typhoon  warnings  costs  $8.3  million 
annually,  exclusive  of  diversion-of-manpower  costs.  Thus,  the  mere  threat 
of  hurricanes  incurs  expenditure  of  prodigious  sums. 

1 • 2  Some  Observational  Facts  on  Hurricanes 

"A  tropical  cyclone  starts  out  as  a  tropical  disturbance  in  which 
there  is  a  slight  surface  circulation  and  perhaps  one  closed  isobar.  When 
the  wind  increases  to  about  20  knots  and  there  is  more  than  one  closed 
isobar  around  the  center,  it  is  called  a  tropical  depression.  When  the 
wind  rises  to  more  than  34  knots,  and  there  are  several  closed  isobars, 
it  becomes  known  as  a  tropical  storm.  If  the  winds  exceed  64  knots  (74 
miles/hour),  it  is  classified  as  a  hurricane  or  typhoon  or  cyclone 
(depending  on  location)"  (Day,  1966,  p.  187),  Thus  any  low-latitude  low- 
pressure  circulation  is  technically  a  tropical  cyclone;  however,  probably 
to  emphasize  the  physical  similarity  of  all  very  intense  tropical  lows 
despite  the  plethora  of  local  appellations,  in  practice  tropical  cyclone 
often  refers  to  the  hurricane  stage  only.  Following  this  practice,  the  name 
tropical  cyclone  will  be  mainly  used  co  refer  categorically  to  the  hurricane 
stage;  hurricane  will  also  be  used  as  a  synonym  for  variation.  The  other 
local  names  will  be  used  mainly  when  geographically  appropriate.  The  local 
designations  for  tropical  cyclones  (Fig.  1)  are  hurricanes  (North  Atlantic 
Ocean),  typhoons  (western  North  Pacific  Ocean),  papagallos  (eastern  North 
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and  cyclo 


Pacific  Ocean),  baguios  (Philippine  Islands),  cyclones  (North  Indian  Ocean), 
trovadoes  (near  Madagascar),  and  willy-willies  [near  Australia  -  a  term 
in  use  in  the  early  twentieth  century,  although  now  it  is  reserved  for 

dust  devils  (Spark,  1972;  Vollsprecht,  1972)]. 

Annually,  especially  in  the  late  summer  and  early  fall,  some  dis¬ 
turbances  over  tropical  oceans  warmed  by  solar  radiation  (usually  at  least 
26° C  -  27°C,  often  28°C  and  higher)  intensify  into  tropical  cyclones, 
which  are  typically  at  least  a  thousand  miles  in  diameter  and  ten  miles  in 
height  (i.e.,  they  extend  from  sea-level  to  tropical  tropopause) .  Not 
only  do  most  typhoons  form  ir  the  autumn,  but  also  the  most  intense  ones 
occur  then  (Brand,  1970a).  These  vortical  storms  are  cyclonic  in  the 
Northern  Hemisphere  and  anticyclonic  (North  Pole  reference)  in  the  Southern 
Hemisphere,  and  take  many  days  to  intensify  --  indicating  that  the  small 
Coriolis  force  is  the  source  of  angular  momentum  and  explaining  why  inten¬ 
sification  within  5°  of  the  geographic  equator  is  very  rare.  In  fact, 
conservation  of  angular  momentum  in  itself  indicates  that  a  fluid  particle 
in  the  tropics  drawn  in  about  five  hundred  miles  will  swirl  at  several 
hundred  miles  per  hour. 

With  satellite  photography,  inspection  of  broad  ocean  exnanses  has 
improved,  and  some  former  estimates  about  the  frequency  of  tropical 
cyclones  have  had  to  be  revised.  About  three-quarters  of  the  annual  global 
total  of  fifty  tropical  cyclones  occur  in  the  Northerm  Hemisphere.  In 
the  North  Atlantic  about  60%  of  the  nine  tropical  storms  that  typically  occur 
annually  intensify  into  tropical  cyclones;  in  the  eastern  North  Pacific, 

33%  of  14;  and  in  the  western  North  Pacific,  80%  of  30  (Atkinson,  1971). 

Since  there  are  daily  disturbances  in  the  tropics  in  autumn,  a  weak  dis¬ 
turbance  has  a  poor  chance  M0%)  of  becoming  a  trooical  storm,  but  any 
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disturbance  that  does  manage  to  become  a  tropical  storm  has  a  good 
chance  to  become  a  h^.ricane  (Palmen  and  Newton,  1969,  p.  503).  Since 
globally  just  over  60%  of  tropical  storms  become  hurricanes,  there  appears 
to  be  no  criterion  such  that  once  a  developing  depression  exceeds  it,  the 
depression  will  definitely  become  a  hurricane.  No  one  today  can  infallibly 
predict  which  tropical  disturbances  will  intensify  into  hurricanes. 

1 . 3  Tropical  Cyclone  Generation 

Hurricanes  form  where  there  is  sustained  local  convective  activity 
over  warm  tropical  seas,  with  surface  temperature  27°C  or  higher  (so  that 
low-level  air  lifted  on  a  moist  adiabat  remains  convectively  unstable 
relative  to  the  undisturbed  ambient  up  to  12  km);  where  there  is  enhanced 
cyclonic  shear  (as  occurs  when  the  Intertropical  Convergence  Zone  lies  at 
a  considerable  distance  from  the  geographical  equator);  and  where  there 
is  weak  vertical  shear  of  the  zonal  wind.  Each  of  these  three  criteria 
warrants  elaboration. 

The  first  criterion  regarding  stability  tends  to  be  fulfilled  only 
in  the  fall,  and  is  discussed  in  detail  in  the  discussion  of  tronical- 
atmosphere  stability  (Section  2.1).  With  regard  to  the  need  for  increased 
cyclonic  shear  as  well  as  warm  ocean  temperatures,  it  is  noteworthy  that, 
globally,  65%  of  the  disturbances  that  later  become  tropical  storms  were 
first  detected  between  10°  and  20°  of  the  equator,  with  only  13%  poleward 
of  this  region  and  only  22%  equatorward  (Atkinson,  1971).  The  absence-of- 
vertical-wind-shear  criterion  may  explain  the  anomalous  cyclone  season 
for  the  northern  Indian  Ocean  (in  the  Arabian  Sea,  the  Bay  of  Bengal,  and 
even  the  South  China  Sea);  in  these  areas,  rather  than  a  single  autumnal 
peak,  twin  peaks  in  cyclone  frequency  occur  in  fall  and  spring,  with  a 
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relati/ely  uneventful  summer  season  (Palmen  and  Newton.  1969).  Gray  (1968) 
suggests  that  unless  cumulonimbi  retain  vertical  integrity,  the  atmospheric 
lightening  associated  with  them  is  dissipated  by  being  advected  in  different 
directions  at  different  altitudes  ("ventilated").  Climatological ly,  regional 
differences  in  vertical  wind  shear  within  the  tropics  become  most  evident 
in  the  upper  troposphere.  Gray  (1968)  also  asserts  that  the  up-  and  down- 
drafts  in  cumulonimbi  themselves  help  suppress  enhancement  of  vertical 
wind  shear  as  baroclinicity  increases  during  intensification  of  tropical 
depressions.  Vertical  wind  shear  at  upper  tropospheric  levels  prevents 
the  higher  structure  of  some  typhoons  to  get  fully  organized,  although 
the  structure  may  be  well  defined  in  the  lower  troposphere;  such  so-called 
"shal low  typhoons  "are  invariably  of  minimal  hurricane  i ntensi ty~TVar*ga r 

Particular  features  of  the  North  Atlantic  hurricane  season  are 
now  enumerated.  As  the  peak  of  the  hurricane  season  approaches,  the  region 
where  tropical  storms  reach  hurricane  intensity  moves  eastward  from  the 
Gulf  of  Mexico  and  the  Caribbean  to  the  Cape  Verde  Islands;  as  the  peak 
hurricane  season  passes,  the  spawning  ground  moves  westward  again  to  the 
Caribbean  (Meyer,  1971a).  Coincidentally,  there  is  a  latitudinal  movement 
northward  in  the  first  half  of  the  season,  then  a  retreat  equatorward 
(Riehl,  1954,  p.  323).  Cyclogenesis  poleward  of  20°N  is  a  particular 
characteristic  of  the  North  Atlantic;  in  fact,  for  the  past  four  years 
the  tropical  North  Atlantic  has  been  mostly  free  of  hurricanes,  which  have 
been  forming  at  25°N  and  higher  (Simpson  and  Frank,  1972). 

At  one  time  discussion  of  tropical  cyclogenesis  inevitably  evoked 
discussion  of  waves  in  the  easterlies  (Riehl,  1954).  But  in  recent  years 
attention  has  focused  on  (1)  the  role  of  the  ICTZ  and  (2)  twice-weekly 
autumnal  disturbances  that  begin  as  large  cyclonic  sandstorms  of  1000  n  mi 
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extent  over  Africa  and  that  drift  westward  at  up  to  10  kts.  There  may 
well  be  two  different  sources  of  hurricane  seedlings,  and  if  so,  the  rela¬ 
tive  roles  of  latent  and  sensible  heat  in  each  may  differ  (Garstang, 

1972). 

1.  Tropical  cyclones  tend  to  form  or>  the  poleward  side  of  the 
equatorial  trough;  in  fact,  80%  -  85%  of  synoptic-scale  tropical  disturbances 
form  within  2°  -  4°  of  the  equatorial  trough  on  the  poleward  side. 
Disturbances  are  rarer,  smaller,  and  weaker  when  the  trough  is  closest  to 
the  equator.  The  trough  is  maintained  by  the  Cl SK  process  (Section  2.3), 
since  low-level  meridional  moist  inflow  sustains  the  persistent  cloudiness 

of  this  perennially  low-pressure  region.  Further,  there  is  low-level 
cyclonic  shear  from  interaction  of  easterlies  poleward  of  the  trough  and 
westerlies  equatorward  of  the  trough.  On  a  zonal ly-averaged  basis  around 
the  globe,  the  equatorial  trough  annually  migrates  from  15°N  to  5°S,  lagging 
the  solar  zenith  by  about  two  months  as  it  does  so  (Riehl,  1972a).  Byers 
(1944)  attributes  the  absence  of  hurricanes  in  the  South  Atlantic  largely 
to  the  failure  of  the  ITCZ  to  become  displaced  south  of  the  equator,  even 
in  February,  at  longitudes  extending  from  the  eastern  Pacific  to  western 
Africa.  Agee  (1972)  has  presented  an  interesting  seauence  of  satellite 
photographs  documenting  a  case  of  tropical  cyclogenesis  in  the  vicinity 
of  the  ITCZ  in  late  July,  1972. 

2.  If  it  is  possible  to  correlate  tropical  North  Atlantic  distur¬ 
bances  with  the  ITCZ,  it  is  also  possible  to  correlate  them  with  disturbances 
first  formed  over  the  mountainous  east  African  bulge,  that  migrate  westward 
(Carlson,  1969).  Half  the  disturbances  over  the  tropical  North  Atlantic 

can  be  so  traced;  further,  half  of  these  disturbances  can  then  be  traced 
across  Central  America  to  the  eastern  Pacific.  Actually  75%  of  eastern 
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Pacific  storms  originate  east  of  Central  America.  A  significant  percentage 
of  midseason  hurricanes  (August  and  September)  have  African  origins; 
analyses  of  dust  samples  taken  on  Caribbean  isles  after  hurricane  passage 
reportedly  confirm  the  African  origin  of  the  storm  systems  (Jennings,  1970). 
An  interesting  case  documented  by  Denny  (1971)  is  an  African  seedling 
identified  on  Septembr'  7,  1971,  which  became  a  depression  on  September 
11  and  Atlantic  Hurricane  Irene  on  September  18.  The  system  crossed 
southern  Nicaragua  on  September  20  and  regenerated  to  eastern  Pacific 
Hurricane  Olivia  (948  mb  central  pressure)  before  dying  on  September  28  near 
Baja  California. 

There  are  many  other  correlations  of  tropica!  cyclone  formation 
that  may  be  attempted;  for  example  Carpenter,  Holle,  and  Fernandez-Partagas 
72)  Trave  leienlly  suggested  a  nrcrjor  peak-in-horr iiarre-'rormat'lbh  hear  a 
new  moon,  a  minor  peak  at  full  moon,  and  minima  at  last  ouarter  and  several 
days  after  first  quarter.  However,  the  correlations  with  the  ITCZ  and 
African  seedlings  just  discussed  have  the  advantage  of  inherent  plausible 
physical  mechanisms.  Discussion  of  modeling  of  tropical  cyclogenesis  will 
be  taken  up  in  Section  4. 

1 .4  Properties  of  Mature  Hurricanes 

Tropical  cyclones  have  a  structure  characterized  in  the  mature 
stage  by  a  relatively  cloud-free  calm  eye  (winds  usually  well  below  15  mph) 
of  about  10-20  mi  radius.  The  eye  is  characterized  by  low  Dressure  at  sea 
level  (often  below  960  mb,  infrequently  below  910  mb)  and  high  temperatures 
aloft  (often  10°C  above  ambient  at  the  same  altitude).  Such  pressure 
drops  are  particularly  spectacular  in  the  tropics,  where  surface  Dressure 
usually  varies  by  little  more  than  0.3%  (Riehl,  1954,  chapter  11).  The 
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eye  is  surrounded  by  an  eyewall,  an  approximately  ten-mile-wide  annulus 
of  intense  convection,  torrential  rainfall,  and  deep  thick  couldiness 
(marked  by  large  numbers  of  cumulonimbi ) .  Outside  the  eyewall  are  convec¬ 
tive  rainbands  that  appear  like  pinwheels  or  logarithmic  spirals  in  some 
satellite  photographs  and/or  radar  displays  taken  from  above  the  storm 
(Figs.  2  and  3). 

The  principal  velocity  component  in  much  of  the  storm  is  azimuthal 
(or  tangential).  Hurricane-force  winds  begin  at  the  eyewall  and  extend 
outward  50  to  70  miles;  winds  usually  fall  to  moderate-gale  level  (^35  mph) 
at  distances  of  100  to  150  miles  from  the  center.  The  vertical  component 
of  velocity  is  largest  in  the  eyewall.  In  the  Northern  Hemisphere,  there 
is.  1ow-1fiy.eX-£.Vclon-fg,,Li7rflOW' '(below  three  km)  and  high-level  outflow  (above 
eight  km,  with  the  maximum  outflow  near  12  km);  toward  the  outer  regions 
(between,  say,  200  to  275  miles  from  the  axis)  the  outflow  turns  from 
cyclonic:  to  anticyclonic  (relative  to  an  observor  rotating  with  the  earth). 
If  one  subtracts  off  the  symmetric  flow  from  the  total  flow,  the  outflow 
layer  possesses  a  horizontal  anticyclonic  eddy  to  the  right  of  the  path 
vector  and  a  cyclonic  eddy  to  the  left;  these  are  seemingly  due  to  ambient 
streaming  around  the  high-level  outflow  (Black  and  Anthes,  1971).  In  the 
midtroposphere  (say,  three  to  eight  km  above  sea  level)  there  is  little 
radial  flow.  While  there  is  much  spray,  the  lowest  few  hundred  feet 
(at  least)  of  the  inflow  layer  remain  cloud-free  in  as  far  as  the  eyewall 
(Riehl,  1954). 

It  is  often  agreed  that  there  is  slow  downward  motion  in  the  eye 
and  in  the  outer  regions  of  the  storm  to  balance  the  strong  updrafts  in  the 
eyewall;  that  radiational  cooling  of  about  1.5°C/day  attends  sinking  in 
the  outer,  less  cloudy  regions;  that,  either  directly  or  indirectly,  the 
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Fig.  2.  NASA  photograph  of  Hurricane  Gladys  west  of  Naples,  Florida  taken 
from  Apollo  7  on  October  7,  1968.  The  maximum  wind  speed  was 
65  kts  at  this  time,  and  the  tropopause  was  at  54,000  ft. 
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Fig.  3.  NOAA  radar  photograph  of  mature  hurricane  revealing  lower-level 
rainbands  outside  the  region  of  overcast  (0/C)  and  upper-level 
cirrus  clouds  of  the  outflow. 


release  of  the  latent  heat  of  condensation  reduces  the  density  in  the 
eyewall  to  p>  tablish  a  large  radial  pressure  deficit  relative  to  ambient 
conditions,  from  hydrostatic  considerations;  that  a  still  further  pressure 
deficit  from  ambient  occurs  in  the  eye  owing  *o  roughly  dry  adiabatic 
recompression  of  air  that  has  arisen  along  the  near-moist-adiabat  vertical 
sequence  of  states  that  characterize  an  appreciable  portion  of  the  eyewall; 
and  that  a  cyclostrophic  balance  (balance  of  radial  pressure  gradient  and 
centrifugal  force),  together  with  an  algebraic  decay  of  the  swirl  with 
increasing  radial  distance  from  the  axis,  yields  on  the  average  a  fairly 
good  estimate  of  the  maximum  swirl  speed. 

Outside  the  eyewall,  the  precipitation  rate  falls  off  roughly 
linearly  with  distance  from  the  center  of  the  tropical  cyclone.  Cloud  cells 
in  the  eyewall  are  typically  5  to  20  km  thick  and  a  few  of  the  smaller  ones  can 
rain  over  15  cm/hr.  Spiral  bands  out  to  150  km  yield  1  cm/’nr,  and  further 
out  0.25  cm/hr  --  though  individual  convective  cells  of  1  km  horizontal 
scale  can  give  much  heavier  rainfall.  In  the  low-rainfall  area  the  precipi¬ 
tation  is  probably  snow  that  turns  to  rain  at  the  melting  level  (Meyer, 

1971a,  1971b).  Incidentally,  clouds  in  a  hurricane  are  probably  strained 
by  the  rapid  swirl  into  the  spiral  band  pattern  seen  on  radar  screens  or 
in  satellite  photographs.  The  spiral  bands  give  visualization  to  parts  of 
the  strain  pattern,  rather  than  streamline  pattern.  A  rainband  persists 
typically  for  one-to-two  hours  (Gentry,  1964). 

In  addition  to  rainfall,  storm  surge,  large  waves,  and  high  winds, 
tropical  cyclones  can  spawn  tornadoes  and  waterspouts  of  moderate  intensity 
(Orton,  1970).  These  are  usually  reported  for  the  region  outside  the 
domain  of  hurricane  winds,  but  this  ma^  only  mean  that  twisters  are  more 
easily  discernible  from  the  general  vortical  intensity  in  the  outer  portions. 


One  would  expect  tornadoes  to  be  most  closely  associated  with  the  intense 
convection  and  large  swirl  of  the  inner  rainbands.  Most  tornadoes  are 
reported  for  the  right  forward  quadrant,  with  respect  to  an  observor  looking 
along  the  direction  of  translation.  This  seems  plausible  because  for 
North  Atlantic  hurricanes  that  quadrant  is  also  the  most  severe  with  rejpect 
to  rainfall,  rainbands,  and  winds  {Hawkins,  1971).  As  previously  noted, 
only  near  the  center  are  tropical  cyclones  axisymmetric  tc  good  approximation; 
near  the  outer  edges  there  is  asymmetry.  The  additive  translational  contri¬ 
bution  to  the  aximuthal  winds  has  been  cited  as  one  plausible  source  of  the 
asymmetry  in  wind  speeds  (Riehl,  1954,  p.  290).  But  Riehl  also  notes  that  the 
largest  radial  inflow  occurs  in  the  right  forward  quadrant  once  translational 
effects  are  subtracted  out.  It  may  be  remarked  that  the  so-called 
beta-plane  effect  (variation  of  the  Coriolis  force  with  latitude  for  fixed 
velocity)  is  not  negligible  over  a  system  as  large  as  a  hurricane,  which  is 
normally  at  least  of  diameter  1000  mi . 

The  longest  recorded  lifetime  for  an  Atlantic  tropical  cyclone  was 
Ginger  (September  5  -  Octooer  5,  1972);  of  this  31 -day  lifespan,  20  were 
spent  as  a  hurricane.  Previously,  Carrie  (1957)  had  lasted  18  days  as  a 
hurricane  and  Faith  (1966)  had  lasted  26  days  as  a  system  (Simpson  and 
Frank,  1972). 

Hurricanes  tend  to  weaken  moderately  rapidly  over  land.  The  central 
pressure  of  Camille  (1969)  rose  from  905  to  990  mb  in  about  13.5  hr  after 
land  fall  (Bradbury,  1972).  Agnes  (1972),  a  minimal  hurricane  with  986  mb 
central  pressure  and  maximum  winds  of  75  mph  with  gusts  to  95  mph  over 
the  Gulf  of  Mexico,  was  downgraded  to  a  tropical  storm  eight  hours  after 
crossing  the  Florida  panhandle  (Anonymous,  1972).  In  an  interesting  report 
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Grossman  and  Rodenhuis  (1955)  cite  Hurricane  Able  (1952)  and  Hurricane 
Diane  (1955)  as  examples  of  hurricanes  which  only  weakly  interacted  with 
their  environment  and  still  maintained  appreciable  energy,  circulation,  and 
precipitation  rates  aftei  passing  inland.  For  example,  24  hours  after 
Diane  passed  inland  over  the  Carol inas  on  August  17,  1955  the  central  pres¬ 
sure  had  risen  from  about  985  mb  to  1J00  mb,  and  the  core  rainfall  had 
decreased  from  9.1  cm/day  to  half  as  much.  Convective  instability  in  the 
core  decreased,  and  maximum  rainfall  occurred  75  mi  from  the  center. 

However,  in  the  next  24  hours  the  central  pressure  began  to  fall  a  few 
millibars  and  the  central  rainfall  at  the  center  of  the  storm  began  to 
increase;  the  convective  instability  in  the  core  increased  again.  Hurricane 
Diane,  although  inland,  encountered  no  major  orographic  changes  during  this 
time  while  moving  northeastward.  Grossman  and  fodenhuis  emphasize  the 
uncertainty  in  numerical  modelina  of  the  role  of  latent-heat  release  in  com¬ 
plicated  situations  in  which  soundings  indicate  layers  of  stable  and  of 
convectively  unstable  air  (Section  3.10).  This  lack  of  knowledge  renders 
detailed  analysis  of  how  tropical  cyclones  die  inland  difficult. 

Matano  and  Sekioka  (1971),  after  examining  several  typhoons  near 
Japan,  seem  to  suggest  that  a  tropical  cyclone  gradually  decays  without 
strong  interaction  with  other  atmospheric  systems  if  it  moves  into  extra- 
tropical  regions  where  there  is  no  marked  midlatitude  barocl ini  city.  If 
there  is  marked  barocl ini  city,  the  tropical  cyclone  will  tend  to  interact 
with  extratropical  cyclones,  either  pre-existing  or  else  induced  on  fronts 
passing  close  to  the  typhoon  core. 
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1.5  Path  Prediction 


Tropical  cyclones  move  westerly  in  the  trades,  usually  at  15  to 
25  mph,  before  turning  poleward,  often  at  greater  translational  speeds. 

Except  when  the  westerlies  are  furthest  north  in  midsummer,  one  frequently 
observes  a  recurvature  eastward  in  the  midlatitudes  along  the  western  side 
of  high-pressure  cells  (Riehl,  1972a);  usually  a  decrease  in  intensity 
follows  recurvature  (Riehl,  1972b).  However,  there  are  so  many  special 
circumstances  that  many  exceptions  could  be  cited  to  virtually  every 
generalization  about  path.  For  example,  coexistent  tropical  cyclones  of 
comparable  size  and  intensity  in  the  same  hemisphere  rotate  about  one 
another  (Fujiwhara  effect)  (Brand,  1970b),  while  binary  systems  in  different 
hemispheres  tend  to  move  parallel  (Cox  and  Jager,  1969)  --  as  suggested  by 
classical  hydrodynamic  potential -flow  theory  for  line  vortices.  Prior 
passage,  within  a  month,  of  a  previous  cyclone  can  also  have  an  effect 
(Brand,  1971).  Tropical  cyclones  can  interact  with  extratropical  cyclones; 
some  typhoons  are  large  enough  to  alter  anticyclonic  highs,  the  result  on 
path  being  similar  to  that  just  cited  for  coexistent  typhoons  in  different 
hemispheres  (Palmen  and  Newton,  1969).  Riehl  (1954,  p.  289)  notes  that 
hurricanes  produce  long  waves  that  travel  three  times  as  fast  as  the  storm 
and  provide  early  warning.  When  hurricane  swell  arrives  at  a  coastline, 
the  normal  wave  frequency  of  10  to  15  per  minute,  is  reduced  to  two  to 
four;  the  wave  direction  can  be  interpreted  to  yield  the  path. 

Tropical-cyclone  p?th  prediction  is  in  an  imperfect  state.  Since 
forecast  errors  involving  the  intensity,  rainfall,  and  movement  of  hurricanes 
can  have  serious  consequence,  the  subject  is  pursued  here  briefly.  For 
1971,  short  term  (12-36  hours)  path  forecasts  by  the  National  Hurricane 
Center  were  best  made  by  the  use  of  past  climatological  and  ana>oque  data 
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(what  did  previous  hurricanes  in  a  similar  situation  do?).  Methods  based 
on  historical  data  are  referred  to  as  "objective."  For  1971,  for  long 
term  (48-72  hours)  NHC  found  methods  based  on  dynamical  principles  superior 
(R.  Simpson,  private  communication).  Today  the  average  error  for  24-hour 
predictions  of  hurricane  movement  is  129  n  mi;  the  average  landfall  error 
for  a  24-hour  prediction  is  ujout  100  n  mi.  The  reduction  is  due  to  the 
closer  monitoring  or  tropical  storms  as  they  approach  the  Atlantic  and  Gulf 
coasts  (Simpson,  1971).  The  less  accurate  forecasts  often  entail  faster 
moving  storms,  and  storms  at  latitudes  poleward  of  the  trades,  where 
recurvature  may  occur. 

Many  dynamical  techniques  treat  the  tropical  cyclone  as  a  point 
vortex  steered  in  a  current,  which  has  been  smoothed  to  remove  the  influence 
of  the  circulation  of  the  storm  itself.  The  steering  current  may  be  the 
flow  at  a  specific  level,  usually  in  the  mid-  or  upper  troposphere  [Byers 
(1944,  p.  447)  cited  10,000  feet].  More  recent  work  emphasizes  the  use  of 
a  steering  layer;  Riehl  (1954,  p.  345)  advocated  the  use  of  a  pressure 
weighted  i  -an  flow  from  the  surface  to  300  mb,  over  a  band  8°  latitude  in 
width  centered  on  the  storm,  to  predict  hurricane  direction  and  speed. 

Today  prognostic  flow  is  sometimes  used  to  predict  path.  Further,  the 
barotropic  model  of  Sanders  and  Burpee  (1968)  averages  over  the  depth  of 
the  troposphere  from  100  to  1000  mb,  and  does  not  involve  reduction  of  the 
tropical  cyclone  to  a  point  vortex.  Despite  such  improvements,  the  out¬ 
look  for  reducing  path  forecast  errors  by  more  sophisticated  dynamical 
prediction  models  is  not  favorable: 

Numerical  models  for  predicting  the  movement  and  develop¬ 
ment  of  hurricanes  remain  a  frail  source  of  guidance  to  the 
hurricane  forecaster  for  three  reasons.  First,  an  error  in 
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direction  of  movement  as  small  as  8  to  10  degrees  --  nominally 
an  acceptable  one  in  a  24-hour  forecast  for  extratropical  storm 
centers  --  can  yield  disastrous  results  in  hurricane  warnings  if 
followed  literally.  Second,  the  performance  of  most  hurricane 
prediction  models  depends  significantly  upon  the  initial  direc¬ 
tion  of  movement  of  the  center,  which  in  turn  depends  upon  an 
exact  knowledge  of  the  current  center  position  and  the  position 
6  and  12  hours  earlier.  The  average  positioning  error  is  more 
than  20  nm,  and  often  leads  to  initial  direction  errors  of 
15-20  degrees.  Finally  the  forecaster  remains  hard  put  to  identify 
and  diagnose  the  frailties  of  numerica1  prediction  models  for 
individual  forecasts.  All  too  often  this  has  led  to  near  abandon¬ 
ment  of  the  guidance  materials  and  the  application  of  empirical 
and  individual  experience  factors  in  decision  making.  (Simpson, 

1971,  p.  1.) 

1 . 6  The  Importance  of  Tropical  Cyclones  in  the  Global  Circulation 

Lorenz  (1966,  p.  409,  418)  suggested  that  hurricanes  were  of 
secondary,  rather  than  primary,  importance  in  the  global  circulation  of 
the  atmosphere  (Fig.  4).  However,  he  noted  that  it  was  unclear  how  great 
a  role  hurricanes  play  in  maintaining  the  currents  of  larger  scale.  He 
also  noted  that  hurricanes  often  do  not  appear  in  numerical  simulations 
of  the  general  circulation  on  advanced  digital  computers.  Of  course, 
appearance  is  hardly  to  be  expected  when  restraints  on  computing  time 
severely  limit  resolution;  for  example,  the  Mintz-Arakawa  model  is  currently 
treated  for  a  five-degree  longitudinal  and  a  four-degree  latitudinal 
spatial  grid  with  a  six-minute  time  step,  such  that  about  twenty  minutes 
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of  computer  time  is  needed  to  simulate  a  day  of  climate  (Ranp,  1970;  Gates, 
Batten,  Kahle,  and  Nelson,  1971).  The  resolution  for  disturbances  of 
dimension  less  than  1000  n  mi  is  poor,  but  two  weeks  of  computer  time  per 
day  of  simulate  climate  would  be  required  on  most  computers  if  satisfactory 
resolution  on  a  100  n  mi  scale  were  sought.  Meteorologists  at  the  Rand 
Corporation  in  Santa  Monica,  California  planned  to  program  the  two-level 
Mintz-Arakawa  model  of  the  general  atmospheric  circulation  for  the  advanced 
ILLIAC  IV  computer.  Even  if  such  plans  are  realized,  a  grid  no  finer 
than  1°  x  1°  seems  practical;  thus,  if  and  when  hurricanes  are  added  to 
the  global  circulation  model,  their  role  may  have  to  be  introduced  by 
parameterization  guided  by  more  fundamental  studies  (such  as  a  subroutine 
with  finer  grid  activated  crt  appropriate  occasions).  After  all,  cumulus 
convection,  radiational  cooling,  and  turbulent  mixing  must  currently  be 
parameterized  into  computer  simulations  of  tropical  cyclones  themselves. 

Manabe,  Holloway,  Jr.,  and  Stone  (1970),  in  their  computer  simula¬ 
tion  of  the  global  circulation,  do  observe  the  two-week  intensification 
of  tropical  disturbances  of  scale  2000  km  -  3000  km  near  +  he  ITCZ  in  regions 
where  tropical  cyclones  are  observed;  further,  there  is  low-level  convergence 
and  upper-level  divergence.  However,  the  resolution  is  not  great,  the 
pressure  fell  only  to  about  980  mb,  and  the  disturbance  appears  more  like 
subtropical  cyclones  (warm  core  above  800  mb,  especially  above  400  mb, 
and  cold  core  below  800  mb).  Appreciable  doubt  remains  whether  these 
disturbances  are  to  be  identified  as  hurricanes;  Bates  (1972)  discusses 
possible  inadequacies  in  the  parameterization  of  cumulus  convection  adopted. 

Some  observational  evidence  has  been  presented  that  in  certain 
regions  during  certain  seasons  hurricanes  do  play  a  role  in  the  global 
circulation.  From  interpretation  of  satellite  photographs  of  tropical 
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storms  and  typhoons  for  1967  -  1969,  Erickson  and  Winston  (1972)  have 
suggested  that  substantial  heat  and  moisture  transfer  from  tropical 
cyclones  to  midlatitude  westerlies  plays  a  role  in  the  autumnal  buildup 
of  the  planetary-scale  circulation.  A  tropical  storm  is  believed  to  trans¬ 
port  energy  from  the  lower  to  upper  tropical  atmosphere;  the  energy  is  then 
carried  to  fronts  in  the  midlatitudes  in  broad  extensive  cloud  bands 
extending  northeastward  over  20°  latitude.  The  warming  so  conveyed  to  the 
extratropics  intensifies  the  midlatitude  zonal  westerlies  around  the  300  mb 
level.  When  the  Hadley-cell  circulation  is  weak  in  late  summer,  perhaps 
hurricanes  do  locally  serve  as  a  substitute  mechanism  for  conveying  energy 
from  the  tropics  to  the  extratropics . 

That  hurricanes  could  locally  and  seasonally  play  such  a  role  seems 

at  least  possible.  First,  some  models  of  hurricane  structure  require 

for  self-consistency  an  export  of  heat  from  the  storm  to  the  rest  of  the 

atmosphere.  For  example,  various  models  require  an  export  of  roughly 

10  ergs/day  (Palmen  and  Riehl,  1957;  Erickson  and  Winston,  1972);  unless 

this  energy  (in  the  form  of  sensible  heat),  produced  by  transformations 

within  the  hurricane,  is  removed  by  other  members  of  the  general  circulation 

via  high-level  advection  on  the  periphery  of  the  storm,  the  models  break 

down  (Riehl,  1954,  p.  338;  Palmen  and  Riehl,  i 9 57 ,  pp.  158-159).  Anthes 

and  Johnson  (1968),  by  applying  the  theory  of  available  potential  energy, 
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estimate  that  about  4  x  10  ergs/day  are  contributed  by  a  hurricane  to 

the  global  scale.  Second,  a  hurricane  has  been  estimated  to  possess  a 

23 

kinetic  energy  comparable  to  that  of  a  hydrogen  bomb,  or  4  x  10'  ergs 
(Battan,  1961;  see  also  Appendix  A),  and  one  has  rained  over  95  inches 
on  one  location  in  four  days  (Silver  Hill,  Jamaica  in  November,  1909) 

(Alaka,  1968).  The  daily  production  of  condensed  water  for  precipitation 


20 


16 

in  a  mature  hurricane  has  been  estimated  at  1.6  x  10  g  (Ooyama,  1969, 
p.  29).  A  weather  system  with  such  energy  and  water  content,  that  extends 
radially  hundreds  of  miles  and  vertically  from  sea  level  to  the  tropical 
tropopause,  and  that  persists  for  weeks,  would  seem  no  local  accident. 

Third,  hurricanes  occur  annually,  mainly  in  the  autumn  after  the  long 
heating  of  the  tropical  oceans  by  solar  radiation,  and  hurricanes  do  turn 
poleward  after  drifting  westward  in  the  trades.  These  facts  invite  the 
previously  stated  speculation  that  hurricanes  are  seasonally  part  of  a 
substitute  mechanism  for  relaxing  energy  poleward  when  the  Hadley-cell 
mechanism  is  not  sufficient.  If  so,  then  the  paradox  discussed  by  Bates 
(1972,  pp.  2,  14)  is  perhaps  resolved:  "The  Hadley  cell  of  the  summer 
heni  sphere  is  weak  or  non-existent.  .  ..  There  seems  to  be  a  paradox  in 
the  fact  that  while  the  Hadley  cell  is  most  intense  in  winter,  the  frequency 
of  oceanic  tropical  disturbances,  which  one  would  expect  to  be  an  important 
contributor  to  its  rising  branch,  is  greatest  in  summer."  Even  if  temporally 
and  longitudinally  varying  eddy  transfer  associated  with  pressure  troughs, 
rather  than  the  zonally  symmetric  mean  meridional  (Hadley)  cell,  is  the 
primary  mechanism  by  which  moisture  and  angular  momentum  are  transported 
from  the  tropics  to  the  midlatitudes  (Riehl,  1954,  chapter  12;  Riehl, 

1969a),  still,  a  significant  portion  of  the  eddy  transfer  occurs  in  the 
upper  troposphere.  The  mode  by  which  quantities  to  be  transferred  are 
convected  to  appreciable  height  in  the  tropics  is  cumulonimbus  clouds,  and, 
as  discussed  later,  these  occur  in  abundance  in  hurricanes.  The  basic 
suggestion  remains:  hurricanes  may  play  some  role  in  the  fall  in  the 
export  of  energy  from  the  tropics  to  the  midlatitudes  by  what  they  export 
themselves,  and  by  what  they  convect  to  the  upper  tropospheric  levels  in 
the  tropics  for  other  mechanisms  to  export. 
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However,  on  an  annual  and  global  basis,  vertical  transfer  in  the 
tropics  is  carried  out  by  smaller-scale,  but  more  numerous  convective 
system,  (Palmen  and  Newton,  1969,  p.  572).  In  fact,  it  is  now  suggested 
that  only  about  one  percent  of  the  total  energy  annually  exported  from 
the  tropics  to  the  midi ati tudes  is  conveyed  by  hurricanes.  For  convenience, 
suppose  in  a  year  there  are  100  hurricanes,  each  lasting  18.25  days  and 
covering  106  mi2.  The  surface  area  of  the  tropics  (with  30°  serving  as 

O  O 

the  limit)  is  about  10  mi  .  The  fraction  of  the  total  annual  total  static 
enthalpy  H  transferred  by  hurricanes  is 

( 100) (18.25) 106  ( 9H/9Z)  . 

R  _ _ OjH  0) 

( 365) 108  (3H/aZ)0ja 

where  the  numerically  unassigned  quantity  in  the  denominator  is  the 
ambient  sea-level  transfer  rate,  and  that  in  the  numerator  is  the  rate 
above  ambient  associated  with  hurricanes  (the  transfer  coefficients  have 
been  cancelled  in  numerator  and  denominator).  According  to  Carrier, 
Hammond,  and  George  (1971),  the  ratio  of  the  two  unassigned  quantities  is 
about  0.2,  so  R  =  0.01.  The  value  of  R  would  be  larger  if  the  model  of 
Malkus  and  Riehl  (1960)  were  used.  These  authors  assert  that  at  the  sea 
surface  in  a  moderate- i ntensi ty  hurricane  (966  mb  central  pressure)  the 
Bowen  ratio  (ratio  of  sensible  to  latent  heat  transfer)  is  about  0.2,  and 
that  the  latent  heat  transfer  is  augmented  by  a  factor  of  10  -  12  over 
the  ambient  level  in  the  trades. 
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2.  ASPECTS  OF  TROPICAL  METEOROLOGY 


2.1  Stability  in  the  Tropical  Atmosphere 

Stability  in  the  tropical  ambient  atmosphere  introduces  consideration 
of  circumstances  rare  in  the  extratropics.  For  convenience  these  properties 
of  the  tropical  atmosphere  curcial  to  tropical  cyclogenesis  are  introduced 
at  this  point. 

Stability  is  conveniently  described  in  terms  of  the  total  static 
enthalpy  H  where 

H  -  c  T  +  gz  +  LY;  (2) 

P 

here  c  is  the  heat  capacity  of  the  atmosphere  (effectively  that  of  dry 
P 

air,  and  independent  of  temperature  for  present  purposes);  T  is  static 
temperature;  g  the  magnitude  of  the  gravitational  acceleration;  z  height 
about  sea  level;  L  the  relevant  latent  heat  of  phase  transition  for  water 
substance;  and  Y  the  mass  fraction  of  water  vapor.  Throughout  this  paper, 
whenever  H  is  (loosely)  described  as  a  temperature,  reference  is  to  (H/cp). 
Actually  the  total  enthalpy  of  a  fluid  particle  is  the  total  stagnation 
enthalpy 

Ht  =  cpT  +  gz  +  LY  +  q2/2  (3) 

where  q  is  the  magnitude  of  the  velocity  vector  (the  wind  speed). 

The  magnitude  of  the  various  contributions  in  the  tropics  to  H  and 
are  worth  mention.  As  shown  in  Fig.  5,  the  static  enthalpy  decreases 
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from  a  sea  level  maximum  of  about  350°K  (L  =  heat  of  condensation  dIus 
fusion)  to  a  midtropospheric  minimum  of  about  330°  at  about  650  mb  (around 
13,000  ft),  before  recovering  the  sea  level  value  of  350°  just  below  the 
tropical  troposphere  (that  is,  at  almost  100  mb  or  50,000  ft).  This  mid¬ 
tropospheric  minimum  is  observed  with  monotonous  regularity  throughout 
the  tropics  in  all  seasons,  although  the  local  sea-level  and  near-tropopause 
maxima  are  reduced  from  the  350°  level  in  the  eouatorial  trough  to  340° 
near  the  subtropics.  Over  the  lowest  few  tens  of  millibars  there  may  be  a 
layer  of  well-mixed  air  in  which  H  is  virtually  constant,  but  this  is  of 
no  major  consequence  here.  The  dynamic  contribution  of  (fl  /2)  is  always 

small,  and  usually  negligible;  for  example,  even  in  a  200-mnh  hurricane, 
o 

(q  /2c  )  =  4°K,  or  just  over  a  1%  contribution  to  (H/cp)  near  sea  level. 

For  stability  one  can  virtually  always  consider  li,  rather  than  H^.  At 
sea  level,  the  static  temperature  contribution  T  is  about  300CK  near  the 
trough,  the  latent  heat  contribution  (LY/cp)  is  about  50°K,  and,  of  course, 
the  gravi tational-potential -energy  contribution  (gz/cp)  is  zero.  The 
moisture  content  of  the  tropical  ambient  falls  off  roughly  exponentially 
with  increasing  height;  the  e-foldinq  distance  is  crudely  650  mb.  Whatever 
modest  variation  in  H  normally  occurs  from  day  to  day  over  a  locale  in 
the  tropics  is  more  likely  due  to  changes  in  water  vapor  content  than  to 
changes  in  temperature  (Johnson,  1969,  pp.  122-123);  the  lapse  rate  alters 
very  little  as  one  moves  from  cloudy  areas  to  clear  areas,  but  there  is 
a  drop  in  the  dew  point  (Gray,  1972b,  p.  14).  In  any  case,  certainly  by 
400  mb  (or  25,000  ft)  for  current  purposes  H  =  0,  where  the  potential- 
temperature-like  quantity  (actually  an  enthalpy)  0  is  defined  by 

0  h  CpT  +  gz.  (4) 
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A  quantity  describing  the  stagnation  potential-like  enthalpy  may  also  be 
def i ned : 

et  =  CpT  +  gz  +  a2/2.  (5) 

At  50,000  ft  the  gravitational  contribution  (gz/c  )  is  about  150°K, 
or  over  40%  of  o  (or  H),  and  hence  no  longer  negligible. 

Clearly  H  plays  a  role  much  like  the  equivalent  potential  temperature 
used  by  most  meteorologists,  but  H  is  preferred  here.  For  one  thing,  what 
a  dry-bulb  and  a  wet-bulb  thermometer  measure  (of  relevance  below)  is 
readily  understood  in  terms  of  H^  and  e^.  However,  unless  L  is  constant, 

H  is  not  rigorously  a  thermodynamic  state  variable.  However,  L  varies 
slowly  with  temperature  and  H  remains  a  very  useful  concent. 

The  condition  that  convection  is  absent  (mechanical  equilibrium)  is 

( 

that  the  entropy  increases  with  height.  If  air  is  approximated  as  a 
perfect  gas,  the  unsaturated  atmosphere  is  stable  if  (3H/3z)  >  0,  convec- 
tively  unstable  if  (3H/3z)  <  0.  If  the  atmosphere  is  saturated,  then  the 
atmosphere  is  stable  if  (3H/3z)  >  0,  statically  unstable  if  ( 3H/3z)  <  0. 
Convectively  unstable  means  a  perturbed  particle  will  resume  its  former 
position  unless  lifted  enough  for  the  onset  of  condensation,  in  which  case 
it  will  rise  until  its  density  discrepancy  relative  to  ambient  is  reduced 
to  *ero  (i.e.,  rise  to  a  new  equilibrium  position).  Statically  unstable 
means  unstable  without  the  requirement  of  a  sufficiently  large  displacement. 
If  the  atmosphere  is  dry,  or  if  for  some  reason  one  wants  to  consider 
stability  excluding  the  role  played  by  condensation  (dry  ascent  as  opnosed 
to  moist  ascent),  then  the  atmosphere  is  stable  if  ( 36/3z)  >  0,  unstable 
f  if  (3B/3z)  <  0. 
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In  the  extratropics ,  H  and  e  both  increase  monotonically  with 
altitude  usually,  so  the  atmosphere  is  stable  to  both  dry  and  moist  ascent 
(Fig.  6).  Only  in  exceptional  circumstances,  as  thunderstorms,  does  will 

penetrative  convection,  with  its  vertical  (as  opposed  to  slantwise)  ascent, 
occur.  In  the  tropics  e  increases  monotonically  with  height,  but  H  (as 
previously  noted)  has  a  midtropospheric  minimum  (Fig.  7);  thus  air  is 
stable  to  dry  ascent,  but  air  in  the  lower  atmosphere  is  unstable  to  moist 
ascent  (convectively  unstable).  How  turbulent  mixing,  cumulus  convection, 
radiational  cooling,  cumulonimbi  (see  below),  and  large-scale  circulation 
maintain  this  condition  of  convective  instability  deals  with  tropical  meteor¬ 
ology  in  general,  and  lies  largely  outside  the  scope  of  this  review  of 
tropical  cyclones  (although  such  distinction  will  probably  soon  prove 
unsound  with  incipient  progress  in  tropical  cyclogenesis).  Actually 
full  understanding  of  how  the  tropical  ambient  is  maintained  does  not 
exist;  the  complexity  is  indicated  by  the  fact  that,  as  will  now  be 
explained,  the  H  profile  generates  different  types  of  clouds,  which  in 
turn  sustain  the  H  profile. 

In  the  convectively  unstable  tropical  atmosphere,  the  larger  the 
lapse  rate,  the  more  suppressed  is  the  cumulus  activity.  Only  as  the  ambient 
lapse  rate  approaches  the  moist  adiabat  is  there  penetrative  convection. 

For  example,  Malkus  (1960)  notes  that  -()T/9z)  6.6°C/km  from  900  to  200  mb 

in  the  1  *opics  normally,  but  -(3T/3z)  ^  6.0°C/km  in  the  inner  rain  area 
of  Hurricane  Daisy  (1958).  There  have  been  many  confirmations  since,  that 
the  less  pronounced  the  midtropospheric  minimum  of  the  total  static  enthalpy 
in  the  tropics,  the  more  convective  activity  is  likely  to  be  present 
(Garstang,  LaSeur,  Warsh,  Hadlock,  Peterson,  1970;  Aspliden,  1971). 
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marked  A  denotes  the  locus  of  the  minimum  of  (H/cn)',  free  convection  can  most  readily  occur  below 
level  A,  and  only  undilute  ascent  can  continue  much  above  A;  undilute  ascent  continues  to  level  B 
where  sea-level  values  are  recovered.  (Ry  permission  from  E.  Palmen  and  C.  W.  Newton,  Atmospheric 
Circulation  Systems,  Academic  Press,  n.  574.) 


given  by  Jordan  (1957).  A  parcel  at  height  zi  will  become  buoyant  at  z0  where 
exceeds  (H/cn)  at  z0(>z-i);  the  ineoualitv  is  satisfied  for  air  in  the  nlanetary  boun 


If  a  glob  of  very  low-level  tropical  air  (buoyant  element)  is  displaced 
vertically  enougi,  for  condensation  to  occur,  ascent  will  continue  many 
kilometers  virtually  to  the  tropical  tropopause,  under  the  stability 
criteria  just  presented,  if  the  ascent  is  rapid  enough  for  no  mixing  or 
radiational  cooling  to  occur  during  ascent.  Near-sea-level  air  would 
rise  dry  adiabatically  until  saturated  (roughly,  50  mb);  thereafter,  enough 
precipi tation  would  fall  out  to  leave  the  air  just  saturated  at  the  local 
pressure  and  temperature,  the  air  retaining  the  latent  heat  of  phase  transi¬ 
tion.  Such  a  locus  of  thermodynamic  states  is  conventionally  referred  to 
as  the  moist  adiabat;  for  rapid  ascent  in  which  ambient  processes  are  too 
slow  to  act,  H  =  const.  (Air  that  has  risen  to  its  level  of  neutral 
buoyancy  will  not  be  unstable  to  descent,  because  any  compressional  heat 
from  work  done  on  the  fluid  by  gravitational  forces  cannot  be  absorbed 
by  the  condensed  water  substance,  which  has  precipitated  out.  However, 
descending  flow  will  accompany  ascent,  end  this  will  be  addressed  in  Section 
2.3.) 

The  core  of  a  towering  cumulonimbus  is  described  by  a  moist  adiabat. 

The  tropical  ambient  must  be  close  to  moist  adiabatic  for  cumulonimbi  to 
be  significant;  this  is  in  fact  the  case  (Riehl,  1969a,  1969b).  If  the 
ambient  were  moist  adiabatic,  cumulonimbi  would  be  unnecessary  because  the 
ambient  would  be  too  unstable  (there  would  be  gross  convective  lifting). 

If  the  ambient  were  far  removed  from  moist  adiabatic,  cumulonimbi  would  be 
nonexistent  because  the  atmosphere  would  be  too  stable  (there  would  be  no 
lifting,  just  heating  up  of  the  air). 
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2.2  Tropical  Cumulonimbi 


About  one  percent  of  the  area  within  the  tropics  is  receiving 

precipitation  on  average;  this  area  is  likely  to  be  in  the  equatorial 

trough,  but  is  not  uniformly  distributed  (Riehl,  1969b).  The  precipitation 

is  concentrated  in  large-scale  disturbances  consisting  of  a  few  cumulo- 

nimbi  immersed  in  a  larger  number  of  cumuli  (Figs.  8  and  9);  Gray  (1972b) 

estimates  that  10%  -  20%  of  the  rain  areas  are  covered  by  active  cumulo- 

nimbi.  Agglomerations  of  cumulonimbi  are  sometimes  described  as  cloud 

clusters.  Cumulonimbi  tend  to  align  themselves  in  east-west  rows,  or 

bands;  the  bands  can  be  tens-to-hundreds  of  kilometers  long  (Kuettner, 

1971),  and  are  typically  20  km  to  50  km  wide  (Charney,  1972).  This  east- 

west  alignment  is  consistent  with  the  low-level  wind  field  of  large-scale 

disturbances  in  the  tropics.  The  ITCZ,  itself  composed  of  one-or-more 

narrow  bands  of  vigorous  cumulonimbus  convection,  is  typically  of  width  300 

km  (Charney.  1971).  It  is  emphasized  that  the  ITCZ,  while  a  region  of 

variable  winds,  appreciable  cloudiness,  large  precipitation,  and  low 

pressure  (Godshall ,  1968),  does  not  consist  of  a  long  unbroken  band  of 

heavy  cloudiness.  Rather,  it  consists  of  intermittent  cloud  clusters 

(with  strong  low-level  convergence,  appreciable  vertical  movement  with 

precipitation,  cyclonic  vorticity,  and  upper-level  outflow  with  possibly 

anticyclonic  relative  voriticity)  with  interstitial  clear  areas  (with 

divergence,  subsidence,  and  anticyclonic  vorticity)  (Holton,  Wallace,  and  Young, 

1971). 4  Gray  (1972b)  also  emphasizes  this  patchwork  nature  of  the  entire 

tropical  belt,  which  consists  of  20%  cloud  area  (which  contains  the 

2 

cumulonimbi,  that  can  rain  at  the  rate  of  2.5  gm/cm  -day),  and  40%  variable 
cloud  area  and  40%  clear  area  (this  80%  of  the  area  yields  negligible 
precipitation).  There  tends  to  be  moist  ascent  in  the  cloudy  areas,  and 
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H/cp  (°K) 

The  total  static  temperature  (H/cp)  over  the  ocean  in  the 
equatorial  trough  and  in  the  tropics  near  the  subtropics,  as  a 
function  of  Dressure.  Both  typical  cumulus  loci  (cu)  and  also 
cumulonimbus  loci  (cb)  are  noted.  (Based  on,  with  permission, 
E.  Palmen  and  C.  W.  Newton,  Atmospheric  Circulation  Systems, 
Academic  Press,  p.  575.) 


NASA  photograph  taken  in  June,  1965  from  Gemini  IV  of  cumulus, 
cumulonimbus,  and  cirrus  clouds  over  the  Pacific  Ocean  off  the 
western  coast  of  Central  America;  the  view  is  northeast  toward 
Mexico. 


Reproduced  from 
best  available  copy, 


dry  descent  in  the  less  cloudy  areas,  though  by  evaporation  of  liquid  water 
and  by  diffusion  of  water  vapor,  water  substance  is  transferred  from  the 
cloudy  to  less  cloudy  regions. 

These  cloud  systems  are  described  here  because  in  some  still 
incompletely  understood  manner,  they  play  a  role  in  tropical  cyclogenesis 
(Pal men  and  Newton,  1969,  p.  585).  Cumulonimbi  are  given  particular 
attention  because  although  they  cover  only  one-thousandth  (Palmen  and 
Newton,  1969,  p.  440)  to  possibly  four-thousandths  (Gray,  1972,  p.  33)  of 
the  area  of  an  equatorial  belt  extending  10°  latitude  on  each  side  of  the 
equatorial  trough,  they  cover  one- tenth  of  the  area  of  the  trough  itself 
(Bates,  1972,  p.  2).  Tropical  cyclogenesis  often  occurs  near  the 
equatorial  trough,  and  cumulonimbi  occur  in  concentration  in  the  hurricane 
eyewall  (Section  4.3).  In  fact,  the  absence  of  cumulonimbi  in  subtropical 
cyclones  occurring  off  western  India,  over  the  Pacific,  and  over  the 
Altantic  serves  as  one  means  of  distinguishing  these  storms  (in  which 
maximum  intensity  is  reached  in  midtroposphere,  and  in  which  a  warm  core 
occurs  at  high  levels  but  a  cold  core  at  low  levels)  from  tropical  cyclones 
(Walker,  1972).  [Simpson  and  Pelissier  (1971)  discuss  various  so-called 
hybrid  storms  in  and  near  the  tropics,  that  may  have  hurricane-force 
winds  at  least  transiently,  but  do  not  have  conventional  hurricane  struc¬ 
ture.]^ 

A  cumulonimbus  (cb)  is  typified  by  a  radius  of  2  km,  a  lifespan 
of  30  to  40  min,  a  height  of  14  km,  a  peak  water  substance  content  of 
4  g/cm3,  maximum  updraft  speed  of  25  m/sec,  a  thermal  anomaly  of  6°C,  and 
a  total  rainfall  of  2.4  cm.  In  contrast,  a  cumulus  (cu)  is  typified  by  a 
r-dius  of  1  km  or  less,  a  lifetime  of  10  to  25  min,  a  height  of  6  km, 
a  peak  water  content  of  1  g/cm3,  maximum  updraft  speed  of  6  m/sec,  a 
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thermal  anomaly  of  4°C,  and  a  total  rainfall  of  0.3  cm  (Lopez,  1972a). 
Cumulonimbi  are  observed  to  overshoot  their  neutral-stability  level  and 
then  oscillate  about  it;  further,  only  the  core  of  a  cloud  rises  undilute 
because  only  the  larger  clouds  reach  upper-troposnheric  levels  (Malkus, 
1960).  Riehl  and  Malkus  (1961)  estimate  that  a  cumulonimbus  conveys  about 
6  x  lO1^  g/sec  and  that  a  fluid  particle  rises  the  full  vertical  extent 
in  about  30  min. 

In  the  summer  tropics  there  is  far  more  local  ascent  and  descent 
than  might  be  suggested  by  the  mean  synoptic-scale  motion  (Johnson,  1969; 
Gray,  1972b);  in  fact,  if  air  in  the  equatorial  trough  rose  only  with  the 
mean  vertical  flow,  the  ascent  would  be  so  slow  that  radiational  cooling 
alone  would  stop  it  in  midtroposphere  (Palmen  and  Newton,  1969)  and  cumulo- 
nimbi  would  not  exist.  Gray  (1972b)  suggests  that  there  is  ten-to-twenty 
times  the  mean  synoptic  convergence  in  the  lower  summertime  tropir.il 
atmosphere,  and  this  is  the  source  of  the  initial  forced  vertical  displace¬ 
ment  that  leads  to  condensation  and  permits  free  convection  to  continue 
undiluted  ascent  in  the  unsaturated'  conditionally  unstable  ambient 
conditions.  Without  one  or  two  orders  of  magnitude  more  up-and-down  recir- 
culatory  local  motion  than  suggested  by  the  mean  synoptic-scale  motion, 
the  water  substance  distribution  in  the  summertime  tropics  (in  which  indi¬ 
vidual  cumulonimbi  can  rain  at  2.5  cm/day  even  though  the  evaporation  rate 
at  the  sea  surface  is  a  uniform  0.5  cm/day)  is  inexplicable.  This  large 
local  recirculation  may  have  implications  for  the  momentum  balance  because 
it  would  seem  to  inhibit  vertical-wirid-shear  development. 

Malkus  and  Riehl  (1961)  described  the  tropical  cumulonimbus  as  a 
hot  tower  because  sensible  heat  from  condensation  was  supposedly  available 
to  warm  the  ambient  surrounding  air.  Lopez  (1972b).  however,  asserts 
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that  cooling  tower  would  be  more  apropos  because  the  condensational  heat 
is  primarily  used  to  raise  the  air  in  the  column  to  its  level  of  neutral 
stability  ( i . e . ,  the  heat  goes  into  potential  energy).  The  air  that  is 
detrained  from  the  cloud  serves  to  cool  the  environment  because  the  heat 
extracted  from  the  ambient  to  re-evaporated  condensed  moisture  in  the 
detrained  air  exceeds  any  sensible  heat  transferred  to  the  surrounding 
environment.  Measurements  around  cumulus  activity  indicate  that  the  static 
temperature  often  falls  (e.g.,  Riehl,  1969b,  p.  588).  Of  course,  the 
ascent  in  a  cloud  cluster  does  indirectly  engender  drying  and  warming 
through  the  subsidence  of  surrounding  air  that  inevitably  accompanies 
ascent  in  the  cloud;  this  heating  counteracts  heat  loss  to  detrainment 
and  radiation.  Thus  on  a  large  scale  cumulonimbus  activity  generates  a 
net  warming,  but  the  effect  on  the  locally  surrounding  air  is  cooling 
(Gray,  1972b).  Most  models  parameterizing  the  role  of  cumulus  convection 
in  large-scale  disturbances  have  the  direct  sensible-heat-transfer  model, 
rather  than  the  indirect  warming-by-subsidence  mechanism,  in  mind  (see 
Section  3.10).  The  Riehl-Malkus  hot-tower  model  is  inconsistent  with  the 
existence  of  large  local  recirculatory  motion  (Gray,  1972b). 

2.3  CISK 


Among  the  most  important  and  fruitful  concepts  introduced  in 
tropical  meteorology  in  several  decades  is  an  idea  due  to  Charney  (Charney 
and  Eliassen,  1964;  Charney,  1971);  the  concept  of  CISK  (conditional 
instability  of  the  second  kind),  in  addition  to  whatever  contribution  it 
may  make  to  extratropical  meteorology,  seems  of  great  utility  in  describing 
tropical  phenomena  of  widely  different  scales,  ranging  from  a  single 
cumulonimbus  to  hurricane  rainbands  to  cloud  clusters  to  tropical  cyclones 
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to  the  ITCZ  (the  intertropical  convergence  zone  is,  for  present  purposes, 
taken  as  identical  with  the  equatorial  trough,  the  doldrums,  and  the 
trade  confluence).  In  all  these  phenomena  there  is  low-level  convergence 
with  large  moisture  content  in  the  presence  of  appreciable  Coriolis  force. 
Phillips  (1970)  notes  that  while  the  physical  processes  involved  in  C IS K 
can  be  qualitatively  described,  much  remains  to  be  done  with  respect  to 
detailed  quantitative  description  (Geissler,  1972). 

The  CISK  process  is  the  feeding  of  convective  activity  in  a  swirling 
flow  by  frictional ly-induced  inflow  in  a  surface  boundary  layer.  The  idea 
is  that  rotating  flows  like  that  in  a  local  low  over  the  ocean  tend  to 
suppress  convergence  except  in  frictional  boundary  layers  (about  a  kilometer 
or  two  thick  and  continguous  to  the  ocean).  In  this  frictional  layer  there 
is  radial  inflow  because  the  Coriolis  and  centrifugal  forces  (which  balance 
the  radial  pressure  gradient  above  the  boundary  layer)  are  diminished  by 
friction;  in  the  frictional  layer  radial  advective  transport  responds  to 
the  "partially  unbalanced"  pressure  gradient.  The  radial  inflow  in  the 
boundary  layer  is  furnished  by  slight  sinking  in  the  swirling  flow  above 
the  boundary  layer,  in  those  regions  where  the  swirl  speed  decreases  with 
radial  distance  from  the  center  of  the  low.  In  the  tropics  the  converging 
boundary  layer  flow  is  warm  and  moist,  because  of  transfer  from  the  warn 
ocean  d.iring  inflow  and  because  the  boundary  layer  air  is  warm  and  moist 
at  the  outset.  The  boundary  layer  flow  erupts  near  the  center  of  the  distur¬ 
bance  in  the  region  where  the  swirl  speed  of  the  fluid  above  the  boundary 
layer  increases  with  radial  distance.  In  the  conditionally  unstable 
tropical  atmosphere  the  core  of  the  erupting,  rotating  column  can  rise 
vertically  without  dilution  to  the  tropical  tropopause.  [Buoyant  elements 
from  higher  in  the  boundary  layer  may  have  somewhat  lower  total  static 
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enthalpy,  and  rnay  rise  to  a  level  of  neutral  buoyancy  somewhat  lower  than 
that  of  sea-level  air.  In  fact,  a  buoyant  element  may  contain  air  from  a 
wide  range  of  the  boundary  layer  thickness;  the  core  of  the  element,  which 
does  not  mix  with  outside  air,  may  not  mix  internally  either,  so  there  may 
well  be  shedding  off  of  fractional  pieces  at  various  upper- tropospheric 
levels.  This  refinement  need  not  be  pursued  here.  But  it  should  be  noted 
that  there  is  good  laboratory  evidence  that  entrainment  into  a  buoyant  plume 
is  reduced  appreciably  with  increasing  rotation  rate  of  the  environment, 
for  all  but  very  small  rotation  rates  (Emmons  and  Ying,  1967).]  The 
hydrostatically  computed  weight  of  the  convective  column,  which  thermody¬ 
namically  is  described  by  a  moist-adiabatic  locus  of  states,  is  taken  to 
be  less  than  the  weight  of  a  column  of  nearby  air.  The  difference  in 
weight  of  the  two  columns  implies  a  radial  pressure  gradient,  which  (by 
the  conservation  of  radial  momentum)  in  turn  implies  swirling.  The  swirling 
leads  to  further  downflux  into  the  surface  frictional  layer,  further  inflow, 
further  moist  adiabatic  ascent,  and  hence  maintenance  (or  even  augmentation) 
of  the  radial  pressure  gradient.  A  crucial  additional  point  about  the 
CISK  process  is  that  it  furnishes  the  basis  of  explaining  how,  instead  of 
competing,  cumulus-scale  activity  can  cooperatively  interact  with  cyclone- 
scale  activity  to  their  mutual  enhancement. 

The  emphasis  in  CISK  on  both  Coriolis  force  and  conditional  insta¬ 
bility  suggests  that  both  the  ITCZ  and  tropical  cyclogenesis  occur  near, 
but  not  at,  the  equator.  The  critical  point  of  the  CISK  process  is  that 
the  weight  of  the  convective  column  is  lighter  than  that  of  the  ambient 
gas.  From  a  facile  view,  the  release  of  the  appreciable  heat  of  condensa¬ 
tion  during  moist  adiabatic  ascent  can  only  lighten  the  column,  and  the 
existence  of  a  radial  pressure  gradient  is  obvious.  However,  Gray  (1972b) 
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has  recently  viewed  this  typical  summary  of  the  CIS'K  mechanism  as  an 
incomplete  description  of  certain  tropical  phenomena,  because  it  fails  to 
emphasize  the  actual  indirect  heating  mechanism  associated  with  cumulus 
activity  and  consequentially  fails  to  emphasize  the  large  local  up-and- 
down  recirculation  above  the  mean  synoptic  motion  (Section  2.2). 

Another  problem  with  quantitative  analysis  of  the  CISK  nrocess  is 
that  the  surface  frictional  layer  over  the  oceans  may  not  be  adequately 
described  even  in  time  average  by  a  linear  quasi -Ekman-layer  theory,  and 
particularly  not  at  low  latitudes.  Gray  (1972a)  notes  that  dissipation 
of  kinetic  energy  in  the  tropical  frictional  layer  so  exceeds  production 
that  indeed  there  must  be  mesoscale  and  synoptic-scale  sinking  of  new 
momentum-carrying  fluid  to  sustain  the  boundary  layer.  However,  Gray 
claims  that,  contrary  to  Ekman's  theory,  measured  planetary  boundary-layer 
thickness  does  not  increase  with  decreasing  latitude;  he  suggests  vertical 
stability  and  turbulence  scales,  rather  than  rotation,  determine  the 
boundary-layer  thickness.  Gray's  general  remarks,  it  may  be  noted,  are 
not  addressed  to  the  exceptional  circumstances  of  a  tropical  cyclone. 

Mahrt  (1972a,  1972b)  has  attempted  theoretically  to  demonstrate  the  sensi¬ 
tivity  of  tropical-disturbance  development  to  the  details  of  the  sub-cloud- 
level  moisture  convergence;  especially  equatorward  of  the  ITCZ  where  the 
Coriolis  parameter  is  small  and  its  latitudinal  variation  is  large,  non¬ 
linear  advective  acceleration  may  so  alter  the  boundary- layer  flow  that 
spiraling  may  be  in  the  opposite  sense  to  that  predicted  by  Ekman's  theory. 
The  thickness,  latitudinal  variation,  and  strength  of  cross-isobar  flow 
may  deviate  from  what  would  be  predicted  by  linear  analysis. 
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3.  MODELS  OF  A  MATURE  TROPICAL  CYCLONE 

3.1  Introduction 

In  sorting  out  the  thermohydrodynamics  of  a  tropical  cyclone,  one 
is  faced  with  understanding  the  interaction  of  two  scales  of  phenomena: 
the  larger  cyclone  scale  (^500  to  1000  mi),  and  the  smaller  cumulus  scale 
(^1  mi).  The  cyclone  must  feed  the  cumulus  scale,  which  in  turn  sustains 
the  cyclone  scale,  in  a  cooperative  interdependence.  Such  ideas  suggest 
that  the  CISK  process  is  not  only  relevant  to  the  maintenance  of  an 
individual  cumulonimbus  cloud,  but  also  (properly  scaled)  to  the  maintenance 
of  the  tropical  cyclone  in  the  large  (Section  2.3). 

The  cyclone-scale  thermohydrodynamics  for  the  mature  stage,  taken 
to  be  axisymmetric  and  quasisteady  for  most  purposes  to  a  lowest  order 
of  approximation,  will  be  set  forth  in  terms  of  a  model  introduced  in  the 
last  three  years  by  Carrier  and  his  co-workers  (Carrier,  1970;  Dergarabedian 
and  Fendell,  1970;  Carrier,  Hammond,  and  George,  1971;  Carrier,  1971a; 

Carrier,  1971b;  McWilliams,  1971;  Dergarabedian  and  Fendell,  1972b). 

These  articles  delineate  overall  dynamics  and  thermodynamics,  and  imply 
(in  contrast  to  other  models  to  be  introduced  later)  no  major  augmentation 
of  ambient  sensible  and  latent  heat  transfer  is  needed  to  explain  hurricanes. 

These  articles  emphasize  that  many  details  and  refinements  in  structure 
still  remain  to  be  quantitatively  treated. 
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3.2  The  Carrier  Model 


Subsequently  (Section  4)  transient  analysis  will  be  devoted  to  how 
a  severe  vortical  storm  may  be  generated  in  the  tropics.  Here,  the  quasi - 
steady  mature  hurricane  is  studied. 

The  Carrier  model,  on  the  basis  of  subdividing  the  tropical  cyclone 
into  segments  where  different  processes  and  scales  predominate,  is  a  four- 
part  analysis.  The  four  regions,  indicated  in  Fig.  10  are  the  throughput 
supply  I,  the  frictional  boundary  layer  II,  the  eyewall  and  efflux  region 
III,  and  the  eye  IV.  Some  of  this  subdivision  is  conventional,  some  not. 
Besides  clarifying  locally  dominant  physical  processes,  subdividing  permits 
retention  of  the  minimal  number  of  terms  in  locally  valid  quantitative 
formulation;  this  procedure  simplifies  the  mathematical  solution  in  a  manner 
unavailable  to  any  direct  finite-differencing  of  uniformly  valid  equations. 

The  Carrier  model  is  closed  for  convenience  —  there  is  no  very 
significant  amount  of  mass  convected  across  any  boundary,  although  mass 
may  be  diffused  across  any  boundary.  The  cylindrical-like  volume  encom¬ 
passing  the  entire  storm  has  the  sea  surface  for  its  bottom;  its  sides 
lie  far  enough  from  the  center  (about  500  to  1000  mi)  so  that  the  winds  are 
virtually  reduced  to  ambient,  and  the  swirl  relative  to  the  earth  is  taken 
as  zero  for  compatibility  with  the  Ekman  condition  for  no  radial  inflow 
across  the  outer  boundary.  The  top  of  the  storm  is  taken  to  be  that 
height  (M50  mb)  at  which  sea-level  air  in  the  outer  part  of  the  storm,  if 
lifted  rapidly  so  that  the  total  enthalpy  of  a  fluid  particle  remained  con¬ 
stant  because  relatively  slow  ambient-maintaining  processes  would  not  have 
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Fig.  10.  This  conjectured  configuration  of  a  mature  hurricane  with  rough 
order-of -magnitude  dimensions  is  not  drawn  to  scale.  The  sub- 
domains  are:  I,  throughput  supply,  a  region  of  rapid  swirl  and 
very  slow  downdraft;  II,  frictional  boundary  layer;  III,  eyewall; 
and  IV,  eye.  Across  the  boundary  layer  there  is  about  a  100  mb 
drop,  and  across  I,  a  further  drop  of  about  200  mb;  the  pressure 
at  the  top  of  the  hurricane  is  about  150  mb,  i.e.,  the  top  is 
near  the  tropopause.  The  eye-eyewall  interface  is  taken  to 
slope  outward,  though  the  effect  may  not  be  so  pronounced  as 
sketched. 
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time  to  act,  would  no  longer  be  unstable  relative  to  the  local  ambient 
air.  Such  rapid  lifting  is,  of  course,  the  moist  adiabatic  ascent  dis¬ 
cussed  in  Section  2.1.  This  "instability  lid"  lies  at  so  great  a  height 
that  there  is  virtually  negligible  swirl,  as  explained  below;  the  ambient 
pressure  and  temperature  at  this  height  is  taken  to  describe  all  radial 
positions  at  this  height,  from  the  center  to  the  outer  edge.  Thus,  the 
top  of  the  storm  is  an  isothermal,  isobaric,  constant-altitude  lid  with  no 
water  vapor  content  for  current  purposes. 

Discussion  now  turns  to  describing  each  of  the  four  regions  com¬ 
prising  the  tropical  cyclone  in  some  detail. 

In  region  I  there  is  warm  moist  air  typical  in  stratification  of  the 
ambient  atmosphere  in  which  the  hurricane  was  generated.  This  air  spun  up 
under  conservation  of  angular  momentum  as  it  moved  in  toward  the  axis  of 
symmetry  during  the  formative  stage.  As  the  mature  stage  was  approached, 
a  gradient-wind  balance  of  pressure,  Coriolis,  and  centrifugal  forces 
choked  off  any  further  inflow;  the  inflow  is  only  enough  to  prevent  the 
eyewall  III  from  diffusing  outward,  and  that  requires  only  an  exceedingly 
small  radial  flow.  The  air  in  I,  then,  is  rapidly  swirling,  the  azimuthal 
velocity  component  greatly  increasing  and  the  pressure  greatly  decreasing 
from  the  edge  to  the  center.  Under  such  a  radial  profile  for  the  swirl, 
there  is  a  small  downflux  from  the  throughput  supply  I  into  the  frictional 
boundary  layer  II.  The  small  downflux  leads  to  a  large  net  mass  flux  into 
the  frictional  layer  because  of  the  large  area  involved.  Furthermore,  the 
downdraft  is  only  a  gross  tanporal  and  azimuthal  average  because  locally 
and  transiently  there  is  intense  convective  activity  by  which  clouds  form 
and  rai n  falls. 
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In  the  frictional  boundary  layer  II,  the  only  region  in  which  angular 

momentum  is  not  conserved  but  is  partially  lost  to  the  sea,  there  is 

appreciable  influx.  In  fact,  the  azimuthal  and  radial  velocity  components 

are  of  comparable  magnitude;  typically,  for  fixed  radial  position,  the 

maximum  inflow  speed  at  any  axial  position  in  the  boundary  layer  is  about 

one-third  the  maximum  azimuthal  speed.  [This  fraction  is  about  the  one 

reported  by  Hughes  (1952)  from  flight  penetration  of  hurricanes  at  altitudes 

of  1000  feet  or  less.]  The  vertical  velocity  component  is  much  smaller. 

The  reason  for  the  inflow  is,  as  previously  discussed,  that  the  no-slip 

boundary  condition  reduces  the  centrifugal  acceleration,  and  a  relatively 

uncompensated  pressure  gradient  drives  the  fluid  toward  the  axis  of  symmetry 

(so-called  "tea  cup  effect").  Far  from  the  axis  in  II  the  classical 

balance  of  the  linear  Ekman  layer  (friction,  pressure,  and  Coriolis  forces) 

suffices;  since  the  downdraft  from  I  to  II  is  probably  fairly  independent 

of  radial  position  (especially  far  from  the  eyewall)  for  swirl  distributions 

of  practical  interest,  three-quarters  of  the  flux  inward  through  II  comes 

from  downflux  across  the  interface  between  II  and  I  where  (rQ/2)  <  r  <  rQ, 

where  r  is  the  radial  extent  of  the  storm.  Closer  in  to  the  axis  the 
o 

nonlinear  accelerations,  especially  the  radial  acceleration  of  radial  and 
tangential  momentum  and  also  centrifugal  acceleration,  must  enter. 

As  the  pressure  gradient  in  I  lets  more  air  sink  into  II,  the  influx 

into  II  drives  boundary  layer  air  moderately  rapidly  up  a  cloudy  eyewall 
III.  In  the  eyewall  hydrostatic  and  cyclostrophic  approximations  hold; 
the  locus  of  thermodynamic  states  is  the  moist  adiabat  based  on  sea-level 
conditions  in  the  eyewall  (this  particular  point  of  the  Carrier  model  will 
be  examined  in  the  critique  t,o  follow).  The  swirl  is  reduced  in  III  owing 
to  boundary  layer  friction.  Near  the  top  of  III  the  flow  moves  further  away 
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from  the  axis  of  symmetry  such  that,  in  the  outflow,  the  air  seems  to  an 
observor  on  earth  to  be  rotating  opposite  in  sense  to  the  rotation  in  I  and 
II  {which  is  cyclonic).  The  air  in  III  slips  over  the  «ir  in  I  with  no 
interaction;  there  is  no  large  radial  pressure  gradient  in  much  of  III, 
unlike  I. 

At  low  altitudes  the  eyewall  flushes  moist  air  out  of  IV,  and  at 
high  altitude  rained-out  air  is  entrained  into  IV.  In  time  the  eye  becomes 
drier  and  better  defined;  it  is  the  central  core  in  which  relatively  dry  air 
sinks,  is  warmed  by  compression,  and  is  entrained  out  into  the  eyewall  or 
recirculated  within  the  ey*;  (stagnation  at  the  base  of  the  eye  would  permit 
transport  processes  to  cool  the  eye).  The  relatively  light  eye  permits 
much  greater  pressure  deficits  from  ambient  and  hence  supports  much  higher 
swirling  speeds.  At  a  fixed  altitude,  the  density  in  the  eye  is  less  than 
that  in  the  eyewall,  which  in  turn  is  less  than  that  of  the  ambient  gas  at 
the  storm  edge.  Since  (except  in  cumulonimbi)  the  hydrostatic  approxima¬ 
tion  is  uniformly  valid,  spatially  and  temporally,  in  a  hurricane,  the 
sea-surface  eye  pressure  is  less  than  the  sea-surface  eyewall  pressure, 
which  is  less  than  the  sea-surface  ambient  pressure.  While  Carrier  has 
not  explicitly  written  so,  it  would  seem  that  as  a  general  trend,  the 
smaller  the  eye  radius,  the  more  spin-up  under  conservation  of  angular 
momentum.  Hence  the  more  intense  storm  might  generally  have  a  relatively 

small  eye,  the  eye  broadening  as  the  storm  weakens.6 

Carrier's  model  thus  pictures  the  tropical  cyclone  as  a  once-through 
process  in  which  a  "fuel  supply"  —  the  warm  moist  air  ini,  which  is  part 
of  the  tropical  cyclone  at  its  inception  and  is  convected  with  the  storm  -- 
is  slowly  exhausted.  The  storm  weakens  because  the  air  drifting  down  into 
the  frictional  layer  toward  the  end  is  typical  of  the  higher  tropical 
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envirorment  and  hence  of  lower  total  static  enthalpy.  Eventually 
the  fuel  supply  is  exhausted,  and  the  boundary  between  III  and  I  sinks 
toward  II.  Some  models  (e.g.,  Eliassen  and  Kleinschmidt,  1957)  picture  a 
recirculation  through  the  storm  of  outflow  air;  the  storm  does  not  survive 
long  enough  for  this,  nor  could  such  recycled  air  maintain  the  storm. 

There  is  one  important  omission  to  the  foregoing  description  that  has 
been  intentionally  deferred:  the  energetics  of  the  surface  frictional 
layer  and  associated  questions  of  air/sea  transfer  of  latent  and  sensible 
heat.  The  relevant  quantity  to  consider  is  the  total  stagnation  enthalpy 
(the  sum  of  static  enthalpy,  the  heat  associated  with  condensible  moisture, 
gravitational  potential  energy,  and  kinetic  energy  contributions).  This 
quantity  is  conserved  at  roughly  its  ambient  stratification  throughout 
regions  I  and  II;  therefore,  it  is  described  by  a  profile  that  decreases 
with  height  from  1000  mb  to  650  mb,  and  then  increases  with  height,  as 
mentioned  above.  The  implication  is  that  the  heat  and  mass  transfer  from 
the  ocean  to  the  atmosphere  is  about  the  same  within  the  hurricane  as  in  the 
ambient.  This  transfer  helps  compensate  for  the  rain-out  in  the  spiral 
bands  and  helps  maintain  the  warm,  moist  nature  of  the  air  in  I.  [Occasion¬ 
ally  the  Carrier  model  is  still  grossly  misrepresented  as  proposing  adiabatic 
conditions  (constant  total  stagnation  enthalpy  in  II  so  the  net  heat  and 
mass  transfer  from  sea  to  air  is  zero);  such  a  solution  cannot  possibly 
satisfy  the  parabolic  boundary-value  problem  describing  the  energetics  of 
the  frictional  boundary  layer  because  it  obviously  violates  the  initial 
condition  at  r  =  rQ,  the  outer  edge.  In  fact,  if  the  supplemental  flux 
from  the  ocean  is  entirely  eliminated,  as  from  passage  over  land,  the 
spin-down  time  is  0(a/v  '  ft  '  )  where  the  eddy  viscosity  v  =  10  mi  hr  , 
the  normal  component  of  the  rotation  of  the  earth  n  =  2  x  10  ^  hr  ^ ,  and 
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the  height  of  the  throughput  supply  a  *  1  mi  —  so  the  spin-down  time 
is  half  a  day  to  two  days.'j  The  model  of  total  stagnation  enthalny  fixed 
at  its  ambient  stratification  breaks  down  in  the  eyewall  III;  there  the 
vertical  velocity  component  is  at  least  one,  probably  two  orders  of  magni¬ 
tude  larger  than  the  relatively  small  downdraft  into  the  boundary  layer; 
the  result  is  that  convection  dominates  the  slow  ambient-sustainina  pro¬ 
cesses,  so  for  at  least  portions  of  the  eyewall  the  total  stagnation 
enthalpy  is  virtually  constant  at  its  local  sea-level  value,  which  (as 
just  explained)  is  roughly  its  ambient  sea-level  value. 

3 •  3  Critique  of  the  Carrier  Model 

A  closed  system  seems  naturally  definable  since  a  single  hurricane 
probably  does  not  interact  significantly  with  the  entire  planetary 
atmosphere,  especially  when  structure  rather  than  path  is  under  considera¬ 
tion.  However,  the  closed  system  is  considered  oDtional  by  Carrier  because 
the  same  conclusions  would  also  hold  for  an  open  system,  although  the 
arguments  would  be  far  more  difficult  to  construct.  Clearly  treatment  of 
asymmetric  effects  should  be  a  future  goal  of  the  Carrier  model  since 
complete  elucidation  of  the  spiral-band  phenomena  and  other  properties 
are  probably  unattainable  otherwise.  Initial  emphasis  on  axisymmetry  has 
been  the  traditional  path  of  development  for  tropical-cyclone  models. 

The  sketch  in  Fig.  10  is  schematic,  but  it  will  appear  below  that 
in  lieu  of  a  greatly  augmented  oceanic  sea-to-air  latent-and-sensible  heat 
transfer,  the  Carrier  model  relies  upon  the  pressure  deficit  achievable  by 
dry-adiabatic  recompression  in  the  eye  to  be  available  to  maintain  dynamic 
equilibrium  in  a  system  with  wind  speeds  known  to  reach  200  mph.  Thus, 
some  small  but  finite  outward  sloping  of  the  eyewall  at  least  down  to 
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midtropospheric  levels,  if  not  lower,  is  required  for  internal  consistency 
in  the  Carrier  model.  Such  ideas  were  present  in  the  work  of  Haurwitz 
(1935)  and  in  the  early  papers  of  Palmen  (1948). 

The  particular  emphasis  placed  on  careful  scrutiny  of  the  dynamics 
and  energetics  of  the  frictional  inflow  layer  is  a  major  contribution  of 
the  Carrier  model.  [Incidentally,  the  thickness  of  this  layer  is  too  great 
to  let  increased  frictional  drag  explain  the  accelerated  decay  after  land¬ 
fall,  unless  the  coastline  is  mountainous.  Even  then,  the  orographic 
effects  are  more  likely  to  be  felt  in  deleterious  premature  rain-out  owing 
to  lifting  of  I,  which  would  result  in  relatively  dry  air  descending  into 
the  boundary  layer  and  running  inward  toward  the  eyewall,  so  that  moist- 
adiabatic  ascent  would  no  longer  be  possible.  The  rate  of  decay  of  inten¬ 
sity  owing  to  sea-to-air  enthalpy-transfer  reduction  after  landfall, 
according  to  the  Carrier  model,  was  discussed  earlier.  The  fact  that  fric¬ 
tion  effects  the  low-level  moisture  convergence  suggests  small  increases 
in  friction  may  even  cause  intensification.]  However,  the  most  novel  con¬ 
tribution  of  the  Carrier  model  is  identification  of  the  "fuel  supply"  in 
I,  which  requires  only  ambient-level  sea-to-air  total-stagnation  enthalpy 
transfer  for  maintenance.  That  the  fluid  in  I  does  sink  down  into  II  is 
suggested  by  the  observation  by  Gentry  (1964,  p.  64)  that  at  lower 
altitudes  the  temperatures  are  lower  in  the  outer  rainbands  than  in  the 
surrounding  air.  There  is  little  undilute  ascent  from  the  inflow  to  the 
outflow  layer.  Rather,  even  the  air  which  ascends  in  the  outer  rainbands 
also  descends  again  in  the  storm  area,  and  is  n£t  immediately  carried 
away  through  the  outflow  layer  in  the  upper  troposphere. 

The  grossest  feature  of  the  Carrier  model  is  the  absence  of  any 
refinement  to  the  eyewall  structure.  While  for  decades  moGelers  have 
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presented  mean  soundings  in  the  eyewall  that  suggest  moist  adiabatic 
ascent  from  the  surface  layers  (Palmen  and  Newton,  1969,  po.  477-482),  Shea 
(1972)  emphasizes  that  the  ascent  is  limited  to  cumulonirnbi  that  cover 
only  10  to  20%  of  the  eyewall  area.  Subsidence  occurs  in  the  regions 
of  the  eyewall  outside  cumulonirnbi.  Shea,  from  flight  data  that  include 
intensifying  and  decaying  as  well  as  mature  tropical  cyclones,  asserts 
that  the  relative  humidity  probably  is  not  100%  throughout  the  eyewall, 
and  there  is  a  small  midtropospheric  minimum  in  mean  vertical  profiles 
of  the  equivalent  potential  temperature  (far  less  pronounced  than  the 
ambient  minimum).  In  fact,  wet-bulb  effects  may  have  led  to  spuriously 
low  temperature  measurements  at  midtropospheric  heights  (Gentry,  1964). 

In  any  case,  the  failure  to  delineate  this  structure  does  not  have  signifi¬ 
cant  repercussions  for  Carrier's  work  on  the  mature  hurricane,  but  does 
have  important  implications  on  his  work  on  intensification,  to  be  discussed 
later  (Section  4).  Shea  also  suggests  that  the  gradient-wind  balance 
does  not  hold  to  good  approximation  in  the  eyewall;  if  so,  this  is  probably 
due  to  a  contribution  from  the  transient  partial  derivative  of  the  radial 
velocity  component,  and  is  not  taken  to  modify  analyses  of  a  mature 
hurricane  appreciably. 

In  the  following  sections  three  specific  analyses  among  the  many 
carried  out  by  Carrier  and  his  co-workers  to  help  substantiate  the  model 
will  be  briefly  reviewed.  The  three  analyses  involve  maximum  swirl  speed 
estimation,  the  dynamics  of  a  nonlinear  Ekman  layer,  and  the  energetics 
of  the  surface  frictional  layer.  The  quantitative  results  achieved  to 
date  concerning  the  Carrier  model  have  been  attained  without  large-scale 
digital  computation.  The  emphasis  on  subdivisional  investigation  of  the 
four  regions  of  the  storm,  with  interfacial  compatibility,  permits  a 
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substantially  analytic  approach.  It  should  be  noted  that  although  certain 
linearizations  are  sometimes  adopted  for  tractability  by  Carrier  and  his 
co-workers  in  the  course  of  their  analysis,  Carrier's  model  for  the  mature 
hurricane  is  definitely  nonlinear,  and  solutions  derived  by  such  lineariza¬ 
tions  are  acceptable  only  if,  a  posteriori,  they  can  be  demonstrated  to 
satsify  the  original  nonlinear  boundary-value  problem  with  acceptably  small 
error.  Without  prior  proof  of  internal  self-consistency  by  subdivisional 
analyses,  a  full  numerical  treatment  of  the  basic  boundary-value  problem 
would  seem  premature. 

3 . 4  Maximum  Swirl  Speed  Estimate  According  to  the  Carrier  Model 

An  upper  and  lower  bound  on  the  central  pressure  deficit  achievable 
in  a  known  spawning  atmosphere  will  now  be  set  forth  by  use  of  the  Carrier 
hurricane  model,  of  hydrostatics,  and  of  the  thermodynamics  of  moist  and 
dry  air.  Specifically,  the  weights  of  various  columns  of  air  in  the 
storm  will  be  determined  in  light  of  different  moisture  content  and  thermo¬ 
dynamic  processes  involved.  The  bounds  on  the  central  pressure  deficit  can 
then  be  translated  into  an  estimate  of  bounds  on  the  maximum  swirl  speed 
through  dynamics  (the  radial  momentum  equation).  Fletcher  (1955)  had 
suggested  use  of  the  cyclostrophic  balance  once  Dressure  deficits  were 
known,  and  Malkus  (1958)  had  suggested  that  pressure  deficits  could  be 
calculated  from  moist  adiabatic  considerations  for  the  eyewall  and  dry 
adiabatic  considerations  for  the  eye.  Here  the  concepts  are  combined  to 
achieve  quantitative  bounds,  but  just  as  important,  to  demonstrate  that 
hurricane  speeds  could  be  achieved  without  requiring  any  enthalpy  transfer 
from  the  ocean  greatly  in  excess  of  the  ambient  transfer,  provided  the 
eyewall  is  not  perfectly  vertical.  Miller  (1958)  performed  somewhat  similar 
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calculations  to  those  to  be  described  and  also  noted  that  the  Riehl-Malkus 
postulate  of  large  oceanic  latent  and  sensible  heat  transfer  to  air  flowing 
in  through  the  frictional  layer  was  not  required  to  explain  low  central 
pressure. 

The  first  step  is  to  neglect  the  frictional  boundary  layer  II,  which 
is  relatively  thin  and  across  which,  except  for  hydrostatic  variations, 
the  pressure  does  not  change  according  to  lowest-order  boundary  layer  theory. 

The  variation  of  pressure  p,  density  p,  and  temperature  T  with  height 
above  the  ocean  z,  for  any  ambient  tropical  atmosphere  in  which  a  hurricane 
forms,  may  be  computed  from 
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(a  =  dry  air); 

(v  =  water  vapor;  a  =  0.622); 
Pa  +  V  pv  =  P(T)  (RH); 
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T  =  f(p),  RH  =  g(p), 


(10) 


where  the  temperature  profile  f(p)  and  the  relative  humidity  (RH)  profile 
g(p)  are  taken  as  known  from  measurement.  The  saturation  pressure  P(T)  is 
well-tabulated  for  vapor  and  liquid  phases  above  freezing,  and  vapor  and 
solid  below  freezing  (Keenan  and  Keyes,  1936);  a  convenient  and  accurate 
expression  for  P(T)  in  mb  was  given  by  Tetens  (Murray,  1967): 


P(T)  =  6.1078  exp 


a (T  -  273.16) 
(T  -  b) 


(11a) 
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a  =  21.8745584 


over  ice; 


a  =  17.2693882 


over  water,  (lib) 


b  =  7.66  )  b  =  35.86  ) 

where  T  is  in  °K.  The  integration  proceeds  from  the  sea  level  upward  in 
altitude  z;  data  typically  extend  from  about  1000  mb  to  150  mb  (Fig.  11). 

The  top  of  the  storm  is  normally  taken  as  the  height  at  which  the  ambient 
total  stagnation  enthalpy  (for  which  the  kinetic  energy  contribution  is 
negligible)  recovers  its  sea-level  value,  as  noted  earlier;  here,  however, 
a  slightly  different  procedure  explained  below  will  be  used.  The  sea-level 
ambient  state  is  henceforth  denoted  by  subscript  s. 

In  a  fully  developed  storm  the  air  rising  up  the  eyewall  ascends  in 
cumulonimbus  clouds,  and  thus  follows  moist  adiabats.  The  initial  states 
to  be  used  for  the  moist  adiabats  are  not  necessarily  known  a  priori. 

If  one  believes  the  total  stagnation  enthalpy  is  constant  along  a  streamtube 
in  the  surface  frictinal  layer  f^om  the  ambient  to  the  eyewall,  then  the 
sea-level  tropical  ambient  state  may  characterize  the  total  stagnation 
enthalpy  of  the  streamtube  rising  in  the  eyewall  closest  to  the  storm 
center  (procedure  A).  [Actually,  mixing  inevitably  occurs  so  perhaps  a 
lower  total  stagnation  enthalpy  characteristic  of  some  height  above  sea 
level  should  be  used,  but  the  small  distinction  is  not  worth  the  effort 
(Carrier,  1971b,  p.  158)].  While  such  a  choice  for  the  initial  state  of 
an  eyewall  moist  adiabat  is  sometimes  made,  and  while  such  results  will  be 
presented  here,  another  preferable  procedure  B  will  also  be  developed.  In 
this  alternative  procedure,  the  temperature  and  relative  humidity  of  the 
sea-level  state  of  the  moist  adiabat  will  be  taken  as  known,  but  the 
initial  pressure  will  be  taken  as  unknown  (to  be  determined  by  iteration  for 
self-consistency  to  be  explained  below). 
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There  is  also  a  comment  worth  noting  concerning  the  equation  that 
describes  the  moist  adiabat.  If  one  simply  takes  dH^  =  0  with  L  held  con¬ 
stant  {Charney,  1971,  pp.  357-358),  where  again  Ht  is  the  total  stagnation 
enthalpy,  then  after  manipulation  with  (6)  -  (9): 
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The  [d(q  / 2 )/dp]  contribution  is  negligible  and  is  henceforth  dropped.  A 
more  careful  derivation  from  basic  principles  (Appendix  B)  yields  a 
slightly  different  expression 
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However,  if  one  takes  p  =  p,  [so  x  =  1  in  (12)]  and  p  =  p  .  which  incurs 

a  a 

an  error  of  only  about  three  percent  at  most  at  any  point  in  the  flowfield, 
then  the  two  expressions  are  equivalent. 

The  procedure  for  computing  the  moist  adiabat  is  to  use  the  dry 
adiabatic  relation  T  ^  p^/ (y-1  ) s  y  =  1.4,  from  the  sea-level  eyewall  con¬ 
ditions  to  that  pressure  at  which  RH  =  1,  at  which  point  one  switches  to 
the  moist  adiabat  and  continues  the  calculation.  The  integration  is 
terminated  at  that  pressure  p-|  for  which  the  temperature  calculated  from 
the  moist  adiabat  and  from  the  ambient  are  equal;  this  temperature  is 
denoted  T-j,  and  the  height  above  sea  level  at  which  T-j  occurs  is  denoted 


z.|  (the  "lid"  on  the  cyclone).  One  then  integrates  (6)  -  (9)  and  the 
differential  equation  for  the  moist  adiabat  from  z  =  z-|  (where  T  =  T-j 
and  p  =  p-j)  to  z  =  0;  RH  =  1  initially  during  this  integration,  but  one 
switches  to  the  dry  adiabat  where  appropriate.  Under  procedure  A  the 
sea-level  pressure  for  the  eyewall  p(z  =  0)  =  pg  is  that  computed.  Under 
the  superior  procedure  B,  an  estimate  of  the  sea-level  eyewall  pressure 
pfi  had  to  be  adopted  to  compute  the  moist  adiabat,  and  this  value  must  be 
recovered  for  convergence.  If  no  eye  existed  in  the  vortex  --  as  seems 
to  be  the  case  for  some  tornadoes  and  waterspouts  --  then  the  just-calculated 
p(z  =  0)  =  pp,  p(z  =  0)  =  pe  would  characterize  conditions  at  the  center 
of  the  vortex. 

The  lid  on  the  storm  as  calculated  here  will  be  higher  in  procedure 
B  than  in  procedure  A,  for  a  given  ambient.  The  taller  the  storm,  the  more 
intense  it  is  according  the  calculations  developed  here.  The  same  relation 
holds  observationally  (Riehl,  1972,  p.  248). 

In  a  mature  hurricane  a  pressure  deficit  in  excess  of  (pg  -  p0)  is 
achieved  by  having  rained-out  air  entrained  from  the  eyewall  sink  in  a 
relatively  dry  eye  under  adiabatic  recompression.  Thus  in  a  hurricane 
(p  -  pg)  is  a  lower  bound  on  the  central  pressure  deficit.  For  an  upper 
bound  on  the  deficit  that  may  be  achieved,  one  may  adopt  the  idealized  model 
that  the  eye  is  completely  dry  (so  no  compressional  heat  is  lost  to  re¬ 
evaporation)  and  that  the  air  entrained  into  the  eye  is  drawn  from  the  top 
of  the  eyewall  (or,  in  any  case,  has  T  =  T-j,  p  =  p-j ,  Y  =  0  at  z  =  z-j ) . 

The  relevant  equations  are  (6)  -  (9),  RH  =  0,  and  T  ^  p^Y  ^Y;  integra¬ 
tion  in  the  direction  of  decreasing  z  yields  p(z  =  0)  =  Pc(<Pe  <  Ps)  and 
p(z  -  0)  =  p  (  <Pe  <  p  )  --  the  density  discrepancies  so  calculated  are  at 
most  25%  and  Fletcher  (1955)  estimates  the  density  does  not  vary  by  even 
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15%,  so  the  density  may  be  held  constant  at  its  ambient  value  throughout 
the  dynamical  calculations  now  discussed. 


If  one  adopts  the  cyclostrophic  balance,  holds  p  constant  at  (say) 
p  ,  and  (since  the  core  is  observed  not  to  rotate)  lets 


(0  0  <  r  <  R 

v(r)  =  {  n  .  (14) 

I.  (v)max  (R/r)P  R  1  r  £  » 


then 


(v) 


max 


(15) 


First,  for  a  one-cell  vortex  (when  there  is  no  eye  so  the  moist  adiabat 
calculation  is  appropriate  all  the  way  to  the  axis),  a  rigid-body-like 
rotation  lies  near  the  core  so  then 


v(r) 


'''lax  <^> 


0  £  r  <  I! 
R  <  r  <  ® 


and  from  the  cyclostrophic  balance 


(v) 


max 


(16) 


(17) 


Clearly,  (15)  holds  for  the  mature  hurricane;  the  maximum  swirl  is  larger 
in  (15)  than  in  (17)  by  the  factor  (n  +  1)^  because  none  of  the  pressure 
deficit  from  ambient  needs  to  be  expended  to  maintain  rotation  of  the  cen¬ 
tral  column.  Next,  although  power-law  decays  of  swirl  with  radial  distance 
are  frequently  adopted  and  suffice  for  current  purposes,  it  will  become 
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evident  that  other  forms  are  at  least  as  plausible,  and  more  convenient,  for 
r  2l  R*  In  any  case,  Miller  (1967)  suggested  from  limited  data  that 
0.5  <  n  <  0.65  usually  suffices,  and  Riehl  (1963)  had  chosen  n  =  0.5; 
earlier,  Byers  (1944,  p.  435)  had  also  recommended  n  =  0.5  and  Hughes  (1952) 
considered  n  =  0.6  -  0.7  suitable  for  an  average  hurricane.  The  upshot 
is  that  estimates  will  be  made  on  the  limits  n  =  0.5  and  n  =  1.0,  although 
n  =  0.5  seems  the  more  realistic.  Finally,  the  gradient  wind  eauation 
would  probably  be  more  appropriate  than  the  cyclostrophic  equation,  but 
the  more  complicated  formula  would  give  maximum  speeds  reduced  by  only  five 
percent  from  those  obtained  from  the  simple  forms  (15)  and  (17). 

First,  results  using  procedure  A  (i.e.,  taking  the  total  static 

enthalpy  H  constant  throughout  the  radial  inflow  as  well  as  e.yewall  region) 

will  be  given.  If  H  is  held  constant  at  sea-level  ambient  conditions 

[pc  =  1014  mb,  T  =  299. 4°K,  (RH)  =  0.84],  then  use  of  (12)  for  the  moist 
s  s  s 

adiabat  gives  condensation  at  972  mb  (T  =  295°K,  z  =  180  ft).  One  finds 
z-|  =  48,150  ft,  p-j  =  138  mb,  T-|  =  203°K.  The  eyewall  pressure  at  sea  level 
pe  =  978  mb  [with  Tg  *  296°K,  (RH)g  =  0.975].  The  eye  pressure  at  sea  level 
p  =  894  mb  (with  T  =  346° K);  this  upper  bound  reflects,  it  is  reiterated, 
the  idealization  of  a  completely  moisture-free  eye.  As  a  variant,  if  sea 
level  ambient  conditions  were  revised  slightly  for  computing  the  moist 
adiabat  only,  [ps  =  1014  mb,  Tg  =  299. 4°K,  (RH)S  =  1.0)],  then  use  of  (12) 
gives  z-j  =  52,500  ft,  p^  =  110  mb,  T-j  =  196. 6°K.  Also,  p0  =  945  mb, 

T0  =  297° K;  pc  =  850  mb,  Tc  =  352. 7°K. 

As  an  illustration  of  the  change  attendant  upon  using  (13)  in  place 

of  (12),  one  finds  for  H  referenced  to  the  unsaturated  sea-level  ambient 

[p  =  1014  mb,  T  =  299. 4°K,  (RH)  =  0.84],  saturation  again  occurs  at 
s  s  s 

972  mb,  z  =  450  ft;  however,  z-j  =  46,500  ft,  p-|  =  150  mb,  T-j  =  205. 4°K; 
further,  pQ  =  988  mb,  T  =  297. 2°K,  (RH)e  =  0.94  (Figs.  11-13).  If  the 
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Fig.  13.  The  density  vs.  altitude  relation  for  each  of  the  three  thermo¬ 
dynamic  loci  of  Fig.  12  are  presented.  The  density  variations 
from  ambient  within  the  storm  are  less  than  30%,  and  speeds  are 
far  below  sonic  everywhere.  < 
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sea-level  ambient  state  is  taken  as  saturated  for  computing  the  moist 

adiabat  only,  z-j  =  50,530  ft,  p-j  =  122  mb,  T  =  199. 6°K;  whence,  pg  =  953  mb, 

T  =  297. 5°K. 
e 

These  results,  to  reiterate,  have  followed  procedure  A  (moist  adiabat 
based  on  sea-level  ambient,  as  given  by  Jordan  (1957)  or  modified  to  be 
saturated  at  the  same  temperature  and  pressure).  The  results  have  been 
presented  here  because  such  calcualtions  have  been  performed  in  the  past, 
and  because  the  size  of  variances  owing  to  use  of  (13)  in  place  of  (12) 


seemed  worth  investigating. 

For  the  superior  procedure  B  (iteration  to  determine  the  consistent 
pressure  at  the  base  of  the  eyewall  moist  adiabat),  only  (13)  will  be  used. 
If  one  adopts  T  =  299. 4°K,  (RH)  =  0.84,  one  finds  convergence  for  po  =  958 
mb  (Figs.  14  and  15);  incidentally,  =  50,100  ft,  p^  =  125  mb,  T^  =  200°K; 
further,  saturation  o:curs  at  918  mb,  z  =  1300  ft.  If  L  is  given  the 


value  of  the  heat  of  condensation  only,  (H/c  )  =  343. 4°K  (sea-level  ambient), 

p  s 

and  ( H/c  )  =  346. 4°K  (sea  level  under  the  moist  adiabat).  [Under  (13), 

P  e 

(H/c  )  is  not  precisely  constant  on  a  moist  adiabat.]  Also,  p  =  875  mb, 

P  C 

Tc  *  349°K.  Thus  it  would  seem  that  the  sea-level  total  static  temperature 
needs  only  to  increase  about  3°K,  with  relative  humidity  and  sea  surface 


temperature  held  fixed,  for  an  appreciably  reduced  pressure  under  eyewall 

cumulonimbi  cores.  According  to  (15),  for  the  deficit  (ps  -  pe)  =  56  mb, 

(v)  v  =  155  mph  for  n  =  0.5  and  (v)m;w  =  219  mph  for  n  =  1.0.  Any  eyewall 

bending  would  permit  appreciably  higher  speeds;  in  fact,  since  (ps  -  pc)  = 

139  mb,  theoretically  sufficient  bending  could  explain  any  speed  up  to 

244  mph  for  n  =  0.5  and  345  mph  for  n  =  1.0.  If  one  takes  again  Tg  =  299. 4°K, 

but  now  lets  RH  =  1.0,  one  finds  convergence  for  pfe  =  890  mb.  Here 

z,  =  57,500  ft,  p,  =  84  mb,  T,  =  191 .8°K.  In  addition,  (H/c  )  =  359. 4°K; 
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Fig.  15. 


In  this  replotting  of  Fig.  14,  explicit  pressure  vs.  altitude 
curves  for  the  ambient,  the  moist-adiabat  taken  to  characterize 
the  eyewall ,  and  the  dry-adiabat  recompression  taken  to  charac¬ 
terize  the  eye  are  presented.  The  moist  adiabat  is  computed  bv 
the  iterative  procedure  described  in  the  cabtion  to  Fig.  14. 


so  an  increase  of  about  16°K  in  the  total  static  temnerature  could  explain 
any  recorded  hurricane  wind  intensity  without  any  eyewall  bending  i_f  the 
sea-level  eyewall  air  were  saturated,  and  if  a  mois*  adiabat  were  taken  to 
be  an  adequate  approximation  to  the  eyewall  sounding.  In  fact,  increases 
above  16°K  appear  impossible  under  these  conditions  as  long  as  the  sea 
surface  is  taken  as  an  isothermal  surface.  Of  course,  use  of  the  moist 
adiabat  for  the  entire  eyewall  and  the  dry  adiabat  for  the  eye  involves 
idealizations  yielding  upper  estimates,  but  the  large  magnitude  of  the 
winds  so  calculated  suggests  the  ideas  behind  the  estimates  are  correct. 

At  this  point  a  brief  summary  of  the  calculated  results  and  their 
implication  seems  apropos.  Even  though  cumulonimbi  occupy  only  10  to  20% 
of  the  eyewall  (Charney,  1971,  pp.  358-359),  results  have  been  presented  for 
moist  adiabatic  ascent  and  hydrostatics  in  the  eyewall  which  suggest  that 
significant  pressure  deficits  relative  to  ambient  may  be  achieved  by  moist- 
adiabatic  ascent  of  air  whose  total  static  temperature  (H/c  )  is  increased  by 

r 

only  one  percent  over  ambient,  and  very  intense  storms  may  be  achieved  by 
ascent  of  air  with  a  five-percent  increase  over  ambient.  The  five-percent 
increase  involves  arbitrarily  increasing  the  eyewall  relative  humidity  at 
shipboard  height  to  100%  from  an  ambient  level  of  84%,  and  it  would  seem 
more  likely  that  the  means  by  which  the  winds  of  very  intense  hurricanes 
are  achieved  is  that  the  eye-eyewall  interface,  as  well  as  it  can  be  defined, 
slopes  at  least  slightly  outward  from  the  top  of  the  inflow  layer  to  the 
top  of  the  storm.  The  calculation  of  wind  speed  for  a  storm  with  an  eye 
implicitly  assumes  such  sloping.  The  additional  dry-adiabatic-compression- 
associated  pressure  deficit  available  to  sustain  swirling  is  reduced  if  the 
eye-eyewall  interface  is  vertical  to  a  certain  height,  and  only  then  slopes 
outward;  the  additional  pressure  deficit  is  unavailable  if  the  eye-eyewall 
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interface  is  perfectly  vertical  virtually  to  the  top  of  the  storm.  (Whether 
or  not  the  eye-eyewall  interface  slopes  outward  does  not  alter  the  conceptual 
thermodynamic  stratification  of  the  eye  nor  the  sea-level  eye  pressure 
deficit  from  ambient.)  Dissenting  publications  will  now  be  discussed. 

Malkus  and  Riehl  (1960)  find  that  moist-adiabatic  ascent  of  ambient 
sea-level  air  to  the  level  of  neutral  buoyancy  produces  a  sea-level  pressure 
of  1000  mb;  insufficient  detail  is  presented  to  ascertain  what  precise 
calculation  was  conducted.  In  any  case,  such  a  pressure  drop  is  insufficient 
to  sustain  even  a  moderate  hurricane,  and  since  these  authors  believe  the 
eye-eyewall  to  be  virtually  perfectly  vertical,  another  means  of  achieving 
and  sustaining  winds  known  to  arise  in  hurricanes  was  required.  The  means 
ultimately  postulated  is  latent  and  sensible  heat  transfer  from  sea  to  air 
greatly  in  excess  of  ambient  transfer  rates;  the  air  ascending  in  the  eyewall 
would  then  possess  large  enough  equivalent  potential  temperature  to  achieve 
the  pressure  deficits  necessary  to  sustain  even  the  highest  known  hurricane 
winds,  lhe  Riehl-Malkus  extra-oceanic-heat-source  postulate  will  be  examined 
below  (Sections  3.7  -  3.9).  Here  it  will  be  noted  that  even  extrapolation 
of  low-speed  transfer  coefficients  to  extreme  hurricane  conditions  could 
not  justify  the  equivalent  potential  temperature  increases  required  by  Riehl 
and  Malkus;  thus  for  very  intense  hurricanes  (only)  an  outward  sloping  of 
the  eyewall  was  conceded  by  these  modelers.  Such  a  procedure  directly  con¬ 
tradicts  the  conclusions  of  Shea  (1972),  based  on  consideration  of  flight 
radar  data;  Shea  finds  that  the  eye-eyewall  interface  slopes  modestly  out¬ 
ward  for  moderate  hurricanes,  but  questions  whether  it  does  so  for  intense 
hurricanes.  However,  radar  data  are  based  on  returns  from  precipitation 
and  err  on  the  side  of  verticality  because  the  strongest  returns  will  come 
from  the  torrential  rains  of  cumulonimbi ,  rather  than  from  weaker  or  decaying 
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cumuli  (Palmen  and  Newton,  1969,  pp.  487-491).  Direct  observational  distinc¬ 
tion  between  an  interfacial  slope  of  a  few  degrees  and  no  slope  at  all, 
especially  when  some  cloudiness  may  occur  in  the  eye,  is  difficult,  and 
observational  evidence  can  be  cited  to  support  either  position.  For 
instance,  Palmen  (1948)  sketches  an  eyewall  with  appreciable  outward  slope 
down  to  nearly  inflow-layer  heights,  based  on  observations  of  a  September 
1947  hurricane  off  Tampa,  Florida.  Palmen  and  Newton  (1969)  construct  a 
schenatic  diagram  of  Hurricane  Helene,  based  on  observational  data  of 
September  26,  1958,  that  indicates  outward  eyewall  slope  commencing  a  few 
kilometers  above  the  sea  surface.  On  the  other  hand,  Riehl  (1954,  pp.  312- 
313)  cites  reports  published  in  1945  that  described  an  observed  eye-eyewall 
interface  that  did  not  have  appreciable  outward  slope  until  30,000  ft. 
Definitive  proof  concerning  whether  or  not  dry-adiabatic  recompression 
contributes  to  the  pressure  difference  available  to  sustain  hurricane  winds 
appears  nonexistent. 

This  section  is  concluded  by  noting  that  maximum  wind  speeds  in 
hurricanes  are  estimated,  in  the  absence  of  aircraft  reconnaissance,  by 
means  of  previously  correlated  formulae  involving  the  area  of  the  overcast 
"circle"  on  satellite  photographs  (Hubert  and  Timchalk,  1969)  (Fig.  3). 

3 • 5  The  Swirl -Divergence  Relation  for  the  Frictional  Boundary  Layer 

Because  the  maximum  speed  achieved  in  a  tropical  cyclone  is  rarely 
much  over  200  mph,  the  Mach  number  rarely  reaches  even  0.3.  Hence,  when 
examining  the  dynamics  (as  opposed  to  the  energetics),  an  incompressible 
constant-property  model  suffices. 

A  steady  axisymmetric  flow  of  an  incompressible  fluid  is  now  studied 
to  confirm  the  crucial  point  that,  under  rapid  swirling,  there  is  downflux 
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from  region  I  to  region  II,  and  sufficient  downflux  enters  the  surface 
frictional  layer  to  account  for  the  mass  flux  up  the  eyewall.  The  analysis 
will  be  carried  out  in  a  noninertial  coordinate  system  rotating  at  the 
constant  speed  of  that  compoennt  of  the  rotation  of  the  earth  which  is  normal 
to  the  local  tangent  plane.  Because  the  boundary-layer  divergence  under 
an  impressed  swirl  (the  major  constraint  furnished  by  the  boundary  layer  on 
the  inviscid  flow  above  it)  is  relatively  small  in  magnitude,  careful 
formulation  and  solution  of  the  coupled  quasi  linear  parabolic  partial 
differential  equations  and  boundary  conditions  describing  the  layer  arp 
required. 

The  relevant  equations  are 

V  •  q  =  0,  (18) 

v(q2/2)  +  (V  x  q)  x  q  +  2fie  x  q  =  -Vp  -  v  vx(v  x  q ) ,  (19) 

where  p  =  (p/p)  +  x  r)2/2  +  gz,  the  gravitational  acceleration  g  =  -gz, 
the  velocity  in  noninertial  coordinates  v  =  x  r  +  a,  the  component  of 
the  rotation  of  the  earth  normal  to  the  local  tangent  plane  =  fiz,  and 
the  kinematic  viscosity  (later  given  eddy-diffusivity  values)  is  v. 

Nondimensionalization  is  effected  by  letting  q'  =  q/Oi^tt) 

P  =  p/ ( H'o^) ,  r‘  =  r/('t'o/fi)^^,  and  E  =  \>/VQ  where  the  Ekman  number  E  «  1 
and  t  characterizes  the  circulation  away  from  the  boundary  (such  as  the 
maximum  swirl  speed  times  the  radius  at  which  it  occurs).  Dropping  primes, 
one  has 

v  •  q,  (20) 

v (q2/2)  +  (v  x  q)  x  q  +  2z  x  q  =  -vp  -  E  vx(v  x  q).  (21) 
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These  equations  are  studies  in  a 'asymmetric  cyclindrical  polar  coordinates 

q=ur+v6+wz,  r  =  r  r  +  z  z.  (22) 

Away  from  the  boundary  (i.e.,  in  region  I)  the  following  expansions 
are  adopted: 

p  =  ^(r,z)  +  v  =  V(r,z)  +  w  =  E^2W(r,z)  +  ....  u  =  o(E1//2). 

(23) 

Substitution  of  (23)  in  (20)  and  (21)  gives  the  gradient-wind  equation: 

"z  =0,  Wz  =  0,  7ir  =  2  V  +  V2/r.  (24) 

Subscripts  r  and  z  here  (and  x  and  c  below)  denote  partial  differentiation. 
The  axially  invariant  solution  is  complete  when  n(r)  or  V(r)  is  specified 
[here  V(r)  will  be  given];  W(r)  is  found  by  matching  the  solution  to 
(19)  to  the  solution  for  the  frictional  layer  II,  and  in  this  sense  W(r) 
is  determined  by  the  boundary-layer  dynamics. 

If  c,  =  z  E-1/2  [which  implies  that  the  frictional  layer  is  0(E^2) 
in  thickness]  and  if  near  the  boundary 

1/2 

U  =  u b(r,0  +  ...,  V  =  vb(r,t)  +  ....  w  =  E  wb(r,t)  +  ...,  p  =  Pb(r,0  +  .. 

(25) 

then  the  axial  component  of  the  momentum  conservation  equation  degenerates 
to  (9pb/9t)  =  0  in  conventional  fashion,  so  the  pressure  field  in  the 
bcundary  layer  is  known  from  (24).  If 

2  %  r\~  1/2  01  /  2  I  0£\ 

ip  =  rvb,  f  =  r  V,  <p  =  rub,  x  =  r,  w  =  2  '  wb,  c  =  2  '  c,  (2b) 
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then  in  terms  of  dimensional  quantities  «f>  =  ru/4^,  ip  =  rv/'i'Q,  v7  =  w/(2fiv)^2, 
t  =  z/(v/2r01/2,  x  =  i^r2/ vi,0 ,  the  boundary  layer  thickness  is  0(v/n)^2. 

In  terms  of  quantities  introduced  in  (26),  one  has  from  (20)  and  (21),  upon 
droppi ng  the  ti Ides , 


*x  +  W(.  =  0, 

(27) 

i|j(,  +  w4  +  (*r  -  i|/  -  if^)/2x  -  (iji  -  t)  -  *  =  0, 

A  L>  GG 

(28) 

Mx  +  wtp r  +  <p  -  iprr  =  0. 

(29) 

Matching  of  expansions  gives 

r,  -*■  »:  4»  -*•  0,  ip  ■+  4'(x)  given, 

(30) 

and  at  z  =  0  no-slip  conditions  are  adopted: 

C  =  0:  4>  =  w  =  ik  -•  0. 

(31) 

Initial  conditions  are  conveniently  given  by  noting  that  at  x  =  xf,  for 
xQ  large  enough,  solution  is  given  by  discarding  all  the  nonlinear  terms 

and  retaining  the  linear  equations  treated  by  Ekman,  in  which  x  enters 

parametrically  only.  The  solution  to  the  balance  of  Coriolis,  pressure, 

and  friction  forces  is  well  known: 

T  =  -[4'(x)]  sin  (2"1/,2c)  exp  {-2~^2z), 

(32) 

+  =  C't'(x)]  [1  -  cos  (2_1/20  exp  (-2“1/2d> 

(33) 

w  =  2'112  [i'x ( x ) ]  |l  -  [sin  (2 :"1/2c)  +  cos  (2’1/2c)]  [exp  (-2~1/20]|. 

(34) 
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Specifically  what  is  sought  is  w(r,H  =  W(r)  for  v(r)  of  interest. 
For  r  large,  from  (34) 

w(r,?-H  =  W(r)  =  2“1/2  4'  (x).  (35) 

Numerical  integration  by  finite-difference  methods  is  formidable 
because  the  flow  component  in  the  time-like  direction,  u,  is,  in  successively 
thinner  strips  lying  parallel  to  the  boundary,  alternately  in  the  direction 
of  integration  (numerically  stable)  and  opposite  to  it  (unstable).  Though 
the  radial  flow  is,  on  net,  in  the  direction  of  integration,  the  integra¬ 
tion  is  marginally  stable.  The  only  finite-difference  results  applicable 
to  the  hurricane  problem  known  to  the  author  were  given  by  Anthes  (1971) 
and  are  discussed  later  (Section  3.10). 

George  (Carrier,  Hammond,  and  George,  1971)  and  Dergarabedian  and 
Fendell  (1972b)  independently  but  simultaneously  applied  the  method  of 
weighted  residuals  to  the  boundary  value  problem,  and  found  that,  except 
near  the  eyewall  where  the  method  was  inadequate,  the  linear  result  given 
in  (35)  sufficed  for  the  nonlinear  problem  as  well  (Fig.  16).  Thus,  for 
a  form  like  4'  =  A(1  -  x/xQ)  [for  which  4>(r0,c)  =  0  according  to  (32)  - 
as  required  by  the  closed-system  model]  w(r,c-*»)  ■+  -(A/xQ),  a  small  negative 
constant  quantity  (equivalent  to  about  0.005  mph  downdraft  for  physically 
interesting  values). 

Incidentally,  the  form  just  mentioned  for  the  rapid  swirl  in  the 
inviscid  flow  above  the  frictional  inflow  layer  is  stable  for  parametric 
values  of  practical  interest,  according  to  Rayleigh's  criterion  (Greenspan, 
1968,  pp.  271-272). 

I 

The  reason  weighted-residual  calculations  fail  near  the  axis,  as 
discovered  by  Carrier  (1971a)  and  McWilliams  (1971)  by  modified  Oseen 
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linearization  and  by  Burggraf,  Stewartson,  and  Belcher  (1971)  by  semi- 
numeri''1  analysis,  is  that  the  structure  of  the  boundary  layer  so  varies 
with  radial  distance  that  adequate  representation  in  terms  of  one  set  of 
orthonormal  polynomials  is  difficult.  Far  from  the  axis  of 
symmetry,  friction  is  important  across  the  entire  layer  of  thickness 
0( near  the  axis,  friction  is  significant  only  in  a  small  sublayer 
near  the  wall  of  thickness  0(r2v/v)1/2,  and  the  remainder  of  the  inflow 
layer  of  0(v/fi)1//2  thickness  is  inviscidly  controlled.  Still,  for  conditions 
of  interest  in  hurricanes,  (35)  is  everywhere  correct  to  within  a  factor 
of  two,  and  often  far  better. 

Still  further  confirmation  of  the  adequacy  of  the  classical  linear 
Ekman  swirl-divergence  relation  for  the  nonlinear  frictional  inflow  layer 
under  a  hurricane  vortex  has  recently  been  given,  in  as  yet  unpublished 
work,  by  Prof.  K.  K.  Tam  of  McGill  University.  Prof.  Tam  has  applied  the 
theory  of  differential  inequalities  to  the  solution  proposed  by  the  Carrier 
group  to  the  boundary-value  problem  (27)  -  (31)  and  indicated  the  satis¬ 
factory  accuracy  of  the  result. 

But  while  the  solution  may  satisfy  the  boundary-value  problem,  one 
may  ask  how  well  the  boundary-value  problem  reflects  the  actual  situation, 
since  a  constant  eddy  viscosity  has  been  adopted  to  model  the  turbulent 
transfer.  It  would  seem  premature  to  adopt  a  second-order  (or  field-closure) 
model  that  reduces  the  Reynolds  time-averaged  formulation  to  a  determinate 
set  by  means  of  a  differential  relationship  between  the  rate  of  change  of 
the  Reynolds  stress  and  diffusion,  production,  and  dissipation  of  turbu¬ 
lence;  still,  an  eddy  viscosity  invariant  with  distance  from  the  air-sea 
interface  seems  outdated,  and  further,  some  small  slip  might  better  model 
the  boundary  condition  at  the  air-sea  interface,  taken  planar. 
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It  should  be  noted,  however,  that  only  one  quanti tative  result  from 
the  nonlinear  Ekman  layer  plays  a  role  in  the  Carrier  model,  and  that  is 
the  swirl-divergence  relation  appropriate  at  the  extreme  outer  edge  of  the 
frictional  inflow  layer.  Furthermore,  the  eddy  viscosity  is  almost  always 
taken  as  independent  of  distance  from  the  "wall"  over  at  least  the  outer 
ninety  percent  of  the  turbulent  boundary  layer  thickness,  aside  from 
interm ittericy  effects,  which  play  very  little  role  in  altering  the  profiles 
of  time-averaged  macroscopic  variables  and  are  very  often  entirely  neglected 
in  calculations.  In  fact,  one  does  not  really  know  how  to  model  slip  and 
effective  roughness  effects  in  the  hurricane  viscous  sublayer,  so  restora¬ 
tion  of  eddy-viscosity  dependence  on  the  normal  coordinate  would  be  a  very 
speculative  undertaking  at  best. 

Other  workers  introduce  a  constant  or  wind-dependent  drag  coefficient 
to  model  frictional  effects  near  the  air-sea  interface  for  computer  solu¬ 
tions.  There  seems  to  be  no  inherent  reason  why  either  of  these  empirical, 
phenomenological  devices  (drag  coefficients,  eddy  viscosity  coefficients) 
is  superior  to  the  other;  either  could  be  made  to  succeed  by  proper  adjust¬ 
ment  (curve-fitting).  Both  face  the  problem  that  empiricism  developed 
for  relatively  small  wind  speeds  and  relatively  smooth  interfaces  must  be 
extrapolated  to  higher  wind  speeds  and  very  rough  seas  (Ooyama,  1969, 
p.  19).  There  exist  attempts  to  study  the  frictional  inflow  layer  entailing 
both  simpler  and  more  detailed  modeling  than  that  just  presented;  for 
example,  Riehl  (1963)  invokes  conservation  of  potential  vorticity  for  the 
frictional  inflow  layer,  and  Riehl  and  Malkus  (1961),  who  assert  that 
dissipation  of  kinetic  energy  in  the  surface  layer  is  independent  of  dis¬ 
tance  from  the  center,  find  that  they  must  adopt  exceptionally  large  eddy 
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transfer  coefficients  to  resolve  problems  that  arise  in  computing  a  refined 
mechanical  energy  budget. 

3 . 6  The  Energetics  of  the  Frictional  Boundary  Layer  and  Throughput  Supply 

For  the  frictional  boundary  layer  II  and  throughput  I,  Carrier, 
Hammond,  and  George  (1971)  take  the  following  approximations  as  adequate 
for  the  quasisteady  mature  phase: 

1.  the  Prandtl  and  Schmidt  numbers  are  unity; 

2.  the  hydrostatic  approximation  holds; 

3.  the  boundary  layer  approximation  holds  (derivatives  normal 
to  the  boundary  exceed  those  tangential  to  the  boundary, 
but  velocity  components  parallel  to  the  boundary  exceed 
those  normal  to  the  boundary); 

4.  the  eddy  transfer  is  adequately  modeled  by  the  laminar  flux- 
gradient  relations  for  diffusion  of  mass,  momentum,  and 
heat  (as  given  by  Fick,  Newton,  and  Fourier,  respectively), 
except  that  the  augmented  kinematic  (eddy)  viscosity  may  vary 
with  radial  position  (but  not  with  axial  position);  and 

5.  the  mixture  of  dry  air  and  water  vapor  may  be  taken  as  a 
perfect  gas  with  constant  heat  capacity  over  the  range  of 
temperatures  of  interest  here. 

In  view  of  the  limited  understanding  of  quantitative  formulation  of 
cumulus  convection  and  turbulent  transfer,  on  the  cyclone  scale  of  interest 
here  these  five  approximations  seem  reasonable.  Clearly  the  familiar 
Reynolds  analogy,  together  with  the  constancy  of  the  eddy  viscosity  over 
most  of  the  boundary  layer,  is  being  adopted.  It  can  then  be  shown  that 
for  L  held  constant,  the  following  equation,  a  generalization  of  ones 
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given  by  Crocco  in  fluid  dynamics  and  Shvab  and  Zel ‘dovich  in  combustion, 
holds  in  the  meteorological  context  of  interest  here: 


Y1  + 

ar 


3H  _  a2H  , 


D, 


(36) 


where  D  denotes  radiational  loss,  taken  by  Carrier,  Hammond  and  George 
(1971)  to  be  representable  as 


•{ 


i 


D  =  -f(z)  H, 

with  f(z)  (to  be  discussed  below)  known.  The  H  used  throughout  this  section 
is  Ht  as  defined  in  (3);  the  subscript  t  is  dropped  to  avoid  double  sub¬ 
scripts  below.  The  boundary-initial  conditions  are  taken  to  be 

z  =  0:  H  =  Hg ;  z  =  z] :  H  »  H] ;  (37) 


r  =  ro:  H  =  Ha(z)> 


(38) 


where,  again,  rQ  is  the  outer  edge  of  the  storm  and  z-|  is  the  top  of  the 
storm.  The  ambient  profile  Ha(z)  is  given;  Hs  =  H-j  =  Hfl(z  =  z^)  =  Hfl(z  =  0). 
Taking  Hr  =  H  (z  =  0)  implies  that  for  r  >  R  (where  R  is  the  radius  of  the 
eyewall),  the  ocean  surface  temperature  is  a  constant  and  the  water  vapor 
mass  fraction  (which  takes  on  its  saturation  value  at  the  nominally  plane 
sea  surface  z  =  0)  is  taken  independent  of  pressure  (see  below).  Since  the 
thermal  conductivity  of  water  greatly  exceeds  that  of  air,  uniform  sea-surface 
temperature  appears  to  be  a  good  approximation.  That  H-|  =  Ha(z  =  0) 
follows  from  the  definition  of  the  lid  on  the  storm. 

The  terms,  from  left  to  right,  in  (36)  represent  radial  advection; 
axial  convection;  turbulent  diffusion  and  cumulus  convection  (both  para- 
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meterized  in  v,  which  will  henceforth  be  treated  as  constant,  though  this 
is  not  necessary);  and  radiation  loss.  Throughout  I,  and  in  II  at  r  =  rQ 
(36)  may  be  approximated  as 

w  -  v  ^  =  -f(z)  H;  (39) 

this  follows  because  <|»(x  )  =  0  [cf.  (29)]  and  u(9H/9r)  is  negligible  in 
I  [cf.  (23)].  In  view  of  the  boundary  conditions  f(z)  must  be  chosen  to 
permit  H(z)  =  H  (z),  where  H  (z)  is  given  (Jordan,  1957);  Carrier,  Hammond, 
and  George  (1971)  take  w(rQ,z)  =  wg(z)  =  0,  but  there  appears  to  be  no  need 
to  require  this.  Thus, 

wa<z>^L-',i4=-f(z)Hai  (40) 

3  Z 

in  fact,  taking  H  =  H  (z)  throughout  I  and  II  satisfies  the  initial  and 
boundary  conditions,  and  renders  u(9H/9r)  =  0.  Further,  the  function  w(r,z) 
is  available  from  the  nonlinear  boundary  layer  dynamics  just  discussed; 
the  dynamics  reveals  (w)  =  w(r ,*;-*»)  =  W(r),  the  value  at  the  outer  edge 

of  the  boundary  layer  --  and  W  is  independent  of  r  over  a  wide  expanse. 

Only  if  W  is  appreciably  different  from  wg  can  H  depart  from  Ha,  and  such  a 
variance  occurs  only  in  the  eyewall,  where  vertical  convection  dominates 
all  other  terms  in  (36)  because  (by  simple  continuity  considerations)  w  is 
increased  two  orders  of  magnitude  over  its  maximum  magnitude  in  I  or  II. 

Hence  for  at  least  portions  of  the  eyewall, 

w  —  =  0,  (41) 

d  Z 
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and  the  relevant  boundary  data  for  this  hyperbolic  suboperator  is  that  given 
at  z  =  0  (37).  ATI  this  discussion,  more  carefully  argued  in  Carrier, 
Hammond,  and  George  (1971),  leads  to  two  significant  results: 

1.  throughout  the  frictional  boundary  layer  and  the  through¬ 
put  supply,  outside  the  eyewall  the  total  stagnation 
enthalpy  is  approximately  fixed  at  its  ambient  stratification. 
Furthermore,  the  small  correction  owing  to  inertial  effects 

in  a  hurricane  is  readily  seen  to  be  a  decrease  of  H  with  i 
such  that  tne  enthalpy  gradient  at  z  =  0  is  increased  slightly 
(about  ten  percent  or  so);  / 

2.  in  portions  of  the  eyewall,  the  total  stagnation  enthalpy  is, 
to  good  approximation,  constant  at  its  sea-level  value, 

i.e.,  columns  of  air  are  very  nearly  rising  on  a  moist  adiabat. 

Numerical  values  of  interest  are  the  ambient  net  sea/air  enthalpy 
transfer  and  the  eddy  viscosity  (Carrier,  Hammond,  and  George,  1971): 


•  5  -2  -1 

-pv[3H(r,z  =  0)/3z]  =  -pv[aH  (r,z  =  0)/3z]  =  1.8  x  10  ergs  cm  sec  ; 


(42) 


v  =  2.7  x  105  cm2  sec”1 


Because  this  result  is  not  compatible  with  many  existing  hurricane 
models,  it  becomes  necessary  to  identify  the  reasons  for  disagreement. 
Hence  in  the  following  sections  the  Riehl-Malkus  theory  is  scrutinized. 

It  way  be  worth  emphasizing  that  (36)  yields  a  steady  solution 
for  the  frictional  inflow  layer  with  all  inertial  terms  considered.  It 
is  not  based  on  obviously  inadequate  models  of  the  inflow  layer,  such  as 
frictionless  or  purely  horizontal  flow. 
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However,  this  result  is  really  not  quite  so  accurate  as  sugqested 
by  Carrier,  Hammond,  and  George,  owing  to  the  omission  of  compressibility 
effects.  A  decrease  in  gas  density  attends  the  decrease  in  gas  pressure 
(of  up  to  12%)  as  one  moves  radially  inward  in  the  frictional  inflow  layer 
from  the  outer  edge  toward  the  eyewall.  This  decrease  in  air  density  is 
negligible  for  the  dynamics,  but  not  the  energetics.  If  the  gas  density 
decreases  as  one  goes  radially  inward,  but  the  total-stagnation-enthalpy 
sea-to-air  transfer  remains  approximately  at  the  ambient  rate  (Carrier, 
1971b,  p.  158  footnote),  then  the  vapor  mass  fraction  (ratio  of  vapor 
density  to  gas  density)  increases  because  there  is  less  mass  to  accept  the 
same  transfer.  Carrier,  Hammond;  and  George  (1971,  p.  162)  acknowledge  the 
effect,  but  fail  to  account  for  it  in  their  calculations.  It  seems  more 
appropriate  to  study  the  vapor  density,  which  at  sea  level  depends  only 
on  sea-surface  temperature,  rather  than  the  vapor  mass  fraction.  A  rough 
first  estimate  is  that  perhaps  a  20%  increase  over  tropical  ambient  levels 
in  latent  and  sensible  heat  transfer  from  the  ocean  to  the  atmosphere  occurs 
near  the  eyewall  of  hurricanes  owing  to  this  density-reduction  effect. 

3 . 7  The  Riehl-Malkus  Postulate  of  an  Oceanic  Heat  Source 

While  the  exposition  of  the  mature  tropical  cyclone  has  been  given 
here  in  terms  of  the  Carrier  model,  most  published  works  adopt  the  framework 
of  the  Riehl-Malkus  theory  alluded  to  earlier  (Riehl,  1954;  Palmen  and 
Riehl,  1957;  Malkus,  1958;  Malkus  and  Riehl,  1960;  Riehl  and  Malkus,  1961; 
Malkus,  1962;  Riehl,  1963).  The  cornerstones  of  this  theory  have  been 
summarized  by  Malkus  (1962,  p.  232): 

The  new  model  relates  core  maintenance  to  mechanism, 

namely,  cumulonimbus  convection  and  sea-air  exchange.  T£ 
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sustain  the  required  pressure  gradients,  two  coupled  pro¬ 
cesses  are  necessary:  first,  a  greatly  magnified  oceanic 
input  of  sensible  and  latent  heat,  and  secondly,  the  un¬ 
dilute  release  of  the  latter  in  concentrated  hot  tower  ascent, 
so  that  air  of  high  heat  content  is  pumped  rapidly  into  the 
upper  troposphere.  The  quantitative  establishment  of  these 
crucial  relationships  was  carried  out  in  a  joint  analytic 
and  observational  framework  (Malkus  and  Riehl,  1960). 

The  important  role  of  large  numbers  of  cumulonimbi  in  the  eyewall 
has  been  confirmed,  although  the  Riehl-Malkus  concept  of  a  direct  hot- 
tower  heating  mechanism  has  since  been  challenged  by  Lopez  (1972b),  as 
discussed  in  Section  2.2.  Here  the  other  major  contention  [".  .  .  it  is 
postulated  that  lowering  of  surface  pressures  in  hurricanes  arises  mainly 
through  an  'extra'  oceanic  heat  source  in  the  storm's  interior"  (Malkus 
and  Riehl,  1960,  p.  12)]  is  questioned.  There  is  undoubtedly  some  augmen¬ 
tation  of  the  ambient-level  sensible  and  latent  heat  transfer  from  sea  to 
atmosphere  over  a  wide  ocean  expanse  under  a  hurricane,  and  the  warm 
autumnal  tropical  seas  are  in  a  nioior  way  responsible  for  the  atmospheric 
total  static  enthalpy  profile  to  begin  with;  however,  an  augmentation  of 
over  an  order  of  magnitude  in  oceanic  heat  transfer  concentrated  according 
to  Malkus  and  Riehl  largely  in  the  region  30  to  90  km  from  the  cyclone 
center  seems  unnecessary  and  unlikely  (Section  3.8). 

The  Riehl-Malkus  theory,  as  noted  earlier  in  Section  3.4,  takes 
the  eye-eyewall  interface  to  be  effectively  perfectly  vertical  so  some 
other  mechanism  than  dry  adiabatic  recompression  must  be  found  to  sustain 
hurricane-force  winds.  Why  this  mechanism  is  stated  to  be  an  "extra" 
oceanic  heat  source  is  understood  by  reconstructing  their  logic.  It  may 


be  worth  noting  at  this  point  that  Riehl  and  Malkus  permit  the  relative 
vorticity  (rather  than  the  relative  velocity)  to  vanish  at  the  outer  edge 
of  the  storm;  workers  using  their  concepts  tend  to  permit  convective,  but 
not  diffusive,  transport  across  the  outer  cylindrical  boundary,  while 
workers  on  the  closed  Carrier  model  permit  diffusive,  but  not  convective, 
transport  across  boundaries. 

Malkus  and  Riehl  (1960,  pp.  3,  7)  acknowledge  the  existence  of  a 
low-level  frictional  inflow  layer,  although  in  a  sample  calculation  there 
is  efflux  from  the  boundary  layer  from  the  outer  edge  of  the  eye  out  to 
500  km  at  an  average  ascent  rate  of  0.7  mph.  Thus  air  flows  in  through 
the  boundary  layer  from  the  outer  edge  of  an  open  system,  rather  than 
sinking  down  into  the  frictional  layer  as  in  Carrier's  model.  The  boundary 
layer  air,  according  to  Malkus  and  Riehl,  undergoes  adiabatic  expansion 
as  it  spirals  inward  toward  the  center,  yet  it  remains  isothermal.  This 
requires  a  vast,  rather  localized  increase  in  sea-to-air  transfer  between 
the  ocean  and  the  contiguous  atmosphere.  That,  the  gradient  normal  to  the 
air/sea  interface  of  temperature  and  of  water  vapor  mass  fraction  is 
large  enough  to  be  consistent  with  vastly  increased  air/sea  transfer  is 
accepted  as  possible: 

In  the  outskirts  of  a  hurricane  the  temperature  of  the 
inflowing  air  drops  slowly  due  to  adiabatic  expansion  during 
(horizontal)  motion  toward  lower  pressure.  It  is  one  of  the 
remarkable  observations  in  hurricanes  that  this  drop  ceases 
at  pressures  of  990  -  1000  mb  and  that  thereafter  isothermal 
expansion  takes  place.  Presumably,  the  temperature  difference 
between  sea  and  air  attains  a  value  large  enough  for  the 
oceanic  heat  supply  to  take  place  at  a  sufficient  rate  to 
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keep  the  temperature  difference  constant.  (Malkus  and 
Riehl,  1960,  p.  9). 

The  actual  transports  [between  sea  and  air],  of  course,  are 

very  large  in  the  hurricane  compared  to  the  trades.  Sensible 

2 

heat  pickup  is  720  cal/cm  /day,  and  increase  by  a  factor  of 

2 

50  over  the  trades.  .  .;  latent  heat  pickup  is  2420  cal/cm  /day, 
higher  by  a  factor  of  12-13.  (Malkus  and  Riehl,  1960,  p.  12). 

Similar  arguments  are  made  elsewhere  by  Riehl  (1954,  pp.  286-287): 

Many  published  records,  notably  those  by  Deppermann.  .  ., 
have  proved  that  the  surface  temperature  outside  the  eye  is 
constant  or  decreases  very  slightly  toward  the  center.  The 
implications  of  this  remarkable  fact  passed  without  notice 
until  Byers.  .  .  drew  attention  to  it.  The  temperature  of 
the  surface  air  spiraling  toward  a  center  should  decrease  if 
adiabatic  expansion  occurred  during  pressure  reduction.  For 
instance,  air  entering  the  circulation  with  the  average  nron- 
erties  of  the  mean  tropical  atmosphere  should  reach  the 
930  mb  isobar  with  a  temperature  of  20.5°C  and  specific  humidity 
of  17  g/kg.  Because  of  condensation,  a  dense  fog  should 
prevail  at  the  ground  inward  from  the  970  mb  isobar.  But  this 
is  never  observed.  It  follows  that  the  potential  temperature 
of  the  surface  air  increases  along  the  inward  trajectories. 

We  also  know  that  the  specific  humidity  increases  and  that  the 
cloud  bases  remain  between  a  few  hundred  and  1000  feet. 

The  surface  air  thus  acquires  both  latent  and  sensible 
heat  during  its  travel  toward  lower  pressure  .... 

i 


so 


A  source  for  the  heat  and  moisture  increment  is  obvious. 

The  ocean  is  greatly  agitated,  and  large  amounts  of  water 
are  thrown  into  the  air  in  the  form  of  spray.  It  is  hard  to 
say  where  the  ocean  ends  and  where  the  atmosphere  begins! 

As  the  air  moves  toward  lower  pressure  and  begins  to  expand 
adiabatically,  the  temperature  difference  between  ocean  and 
air  suddenly  increases.  Since  the  surface  of  contact  between 
air  and  water  increases  to  many  times  the  horizontal  area  of 
the  storm,  rapid  transfer  of  sensible  and  latent  heat  from 
ocean  to  air  is  made  possible.  In  the  outskirts,  say  beyond 
the  990  mb  isobar,  the  turmoil  is  less  and  the  process  of 
heat  transfer  is  not  operative.7 

The  Riehl-Mal kus  theory  that  greatly  augmented  heat  and  mass 
transfer  sustains  the  tropical  cyclone  has,  in  fact,  been  parameterized 
into  all  existing  computer  simulations  (see  Section  3.10).  For  example, 
Rosenthal  (1971b)  closely  reflects  the  Riehl-Malkus  theory  and  notes 
similar  logic  in  the  work  of  another  computer  modeler  of  tropical  cyclones 
(Ooyama,  1969): 

Air-sea  exchanges  of  sensible  and  latent  heat  have  long  been 
considered  important  ingredients  In  the  development  and  main¬ 
tenance  of  tropical  storms.  Palmen  (1948)  showed,  on  a 
climatological  basis,  that  tropical  storms  form  primarily  over 
warm  ocean  waters  (Tgea  >  26°C).  Malkus  and  Riehl  (1960)  showed 
that  the  deep  central  pressures  associated  with  hurricanes  could 
not  be  explained  hydrostatically  unless  the  equivalent  potential 
temperature,  eg,  in  the  boundary  layer  was  10°  to  15°K  greater 
than  that  of  the  mean  tropical  atmosphere.  Byers  (1944)  pointed 
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out  that  the  observed  near -isothermal  conditions  for  inward 
spiraling  air  in  the  hurricane  boundary  layer  required  a 
source  of  sensible  heat  to  compensate  for  the  cooling  due  to 
adiabatic  expansion.  .  .. 

Ooyama  (1969)  found  drastic  reductions  in  the  strength  of 
his  model  storm  when  the  air-sea  exchanges  of  sensible  and 
latent  heat  were  suppressed.  He  pointed  out  that  at  sufficiently 
large  radii,  the  boundary  layer  is  divergent  (the  so-called 
Ekman  layer  "sucking").  .  ..  This  subsidence  tends  to  decrease 
the  boundary  layer  9  since  99  /sz  <  0  in  the  lower  troposphere. 
Ooyama  argued  that  unless  the  energy  supply  from  the  ocean 
can  again  raise  the  9g  of  the  boundary  layer  air  to  sufficiently 
large  values  before  the  inflowing  air  reaches  the  inner  region, 
the  convective  activity  will  diminish  in  those  reqions  and, 
hence,  the  storm  will  begin  to  weaken. 

Ooyama 's  line  of  reasoning  can  be  extended  to  show  that 
evaporation  is  far  more  important  than  sensible  heat  flux.  The 
air  sucked  into  the  boundary  layer  has  a  higher  potential 
temperature  than  the  original  boundary  layer  air.  The  sub¬ 
siding  air  has  a  smaller  9g  only  because  it  is  relatively  dry. 
(Rosenthal,  1971b,  p.  772). 

3.8  The  Intensity  of  a  Tropical  Cyclone  and  the  Underlying  Sea  Surface 
Temperature 

> 

Palmen  and  Riehl  (1957,  p.  156)  state:  "Evidently,  a  cyclone  will 
decay  rapidly  if  it  encounters  thermally  unfavorable  conditions.  The 
generation  term  in  [the  equation  for  conservation  of  kinetic  energy]  will 
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then  decrease  rapidly;  it  will  become  negative  if  the  air  ascending  in  the 
core  is  cooler  than  the  surroundings  in  spite  of  release  of  latent  heat. 

This  would  happen  when  a  storm  moves  over  a  relatively  cold  ocean  surface, 
over  a  continent,  or  when  colder  air  masses  invade  the  cyclone  near  the 
surface." 

Iri  discussing  the  Riehl-Malkus  theory,  one  should  distinguish 
between  what  these  authors  themselves  assert,  and  what  others  citing  their 
work  have  added.  For  example,  Shuleykin  (1970;  1972)  goes  beyond  Riehl 
and  Malkus  and  asserts  that  the  "...  power  of  the  hurricane.  .  .  increases 
sharply  with  the  temperature  of  the  underlying  water  surface"  (Shuleykin, 
1972,  p.  1).  However,  attempts  to  correlate  central  pressure  deficit  and 
sea  surface  temperature,  without  sufficient  account  of  other  factors,  seems 
vnplistic.  This  particular  point  will  be  now  discussed. 

First,  ocean  temperatures  are  difficult  to  determine  accurately  from 

currently  available  records  and  by  currently  available  methods  (Perlroth, 

1967;  Hidy,  1972,  p.  1091).  Next,  Gentry  (1969a,  p.  406)  presents  data 

relating  "...  the  maximum  intensity  of  several  tropical  cyclones  to  the 

temperatures  of  the  sea  beneath  them  and  shows  that  both  severe  and  weak 

tropical  cyclones  occur  when  ocean  temperatures  are  relatively  high.  This 

suggests  that  variations  in  parameters  other  than  the  transfer  of  heat 

from  the  ocean  to  the  atmosphere  also  influence  the  storm's  intensity 

•  •  ••  "  Actually  Gentry's  data  (Fig.  17)  indicate  that  some  intense 

tropical  cyclones  lie  over  relatively  cold  water  (<28.0°C).  While 

Brand  (1971)  cites  a  supposed  correlation  of  central  pressure  deficit  with 

sea-surface  temperature  throughout  the  lifetime  of  Hurricane  Ester  of 

September  9-26,  1961,  it  may  be  noted  that  during  one  twelve-hour  period 

of  constant  sea  temperature,  the  central  pressure  rose  15  mb;  for  three 

days  while  the  sea  temperature  hovered  about  84°F,  the  central  pressure 
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MINIMUM  SEA  LEVEL  PRESSURE  (mbs) 


RELATION  OF  MINIMUM  SEA  LEVEL  PRESSURES 
IN  TROPICAL  CYCLONES  TO  SEA  TEMPERATURES 
(EACH  DOT  REPRESENTS  RELATION  FOR  AN  INDIVIDUAL 
STORM  SELECTED  FROM  1955-67  SEASONS) 
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Fig.  17.  Measurements  of  minimum  sea-level  pressure  vs.  local  sea-surface 
temperature  for  several  tropical  cyclones  (from  Gentry,  1969a, 
p.  406;  reproduced  with  permission  of  the  American  Meteorological 
Society) . 
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nonmonotonically  rose  from  930  mb  to  955  mb;  and  at  various  times  when  the 
sea  temperature  was  at  86°F,  the  central  pressure  was  at  values  as  low 
as  927  mb  and  as  high  as  953  mb.  Perlroth  (1967)  shows  that  for  Hurricane 
Ginny  of  October  21-28,  1963,  for  five  days  while  the  sea  temperature 
hovered  near  80°F,  the  central  pressure  fell  from  995  mb  to  970  mb  (Fig.  18). 
Perlroth  (1969)  cites  the  importance  in  tropical  storm  generation  of 
considering  not  only  high  sea-surface  temperature  (above  26°C)  but  also 
small  oceanic  temperature  variation  with  depth  (less  than  7°C  within  200  ft 
of  the  surface).  Perlroth  emphasizes  that  the  oceanic  environment  is  only 
one  factor  in  tropical  storm  intensity,  and  the  sea-surface  temperature  is 
only  one  factor  in  the  oceanic  environment. 

Finally,  the  fact  that  hurricanes  leave  cold  wakes  in  the  ocean 
might  suggest  that  they  do  drain  the  upper  sea  layers  of  heat.  However, 
the  reason  for  the  sea-surface  temperature  dropping  about  two  degrees 
Centigrade  after  hurricane  passage  is  the  upwelling  of  lower,  colder  water 
owing  to  convection  currents  induced  in  the  ocean  by  the  pressure  deficits 
of  the  hurricane  (Perlroth,  1967).  Slow-moving  intense  tropical  cyclones 
can  cause  upwelling,  from  depths  of  40  to  65  m5  within  the  radius  of  hurri¬ 
cane-force  wi  nds  ;  compensatory  downwelling  of  as  much  as  80  to  100  m  occurs 
45  to  110  n  mi  from  the  path  of  the  storm  (Revesz,  Jr.,  1971;  Leipper,  1967). 
The  warm  ocean  water  near  the  center  of  the  hurricane  flows  radially  out¬ 
ward,  and  is  replaced  by  colder  water  from  lower  depths;  outside  65  n  mi 
from  the  center,  below  the  outflow,  there  is  compensatory  inflow  toward  the 
center  of  the  hurricane  to  complete  the  convection  cell  (Leipper,  1967; 
Wright,  1969). 

Leipper  finds  the  temperature  readings  at  depths  in  the  outer  portion 
of  the  hurricane  path  to  be  higher  than  the  temperatures  before  the  storm. 
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Thus  a  significant  portion  of  the  heat  of  the  upper  layers  would  seem  to 
be  transported  radially  outward,  vertically  downward,  and  then  radially 
inward.  Leipper  and  Wright,  however,  postulate  that  all  the  heat  of  the 
relatively  warm  upper  oceanic  level  in  excess  of  that  of  the  colder  upwelling 

g 

is  transferred  to  the  atmosphere  as  an  augmentation  of  the  ambient  rate; 
Leipper  suggests  a  transfer  rate  of  4500  cal/cm^-da.y  under  Hurricane  Hilda 

p 

(1964),  and  Wright  suggested  7200  cal/cm  -day  under  Typhoon  Shirley  (1965). 

The  previously  cited  figure  of  Malkus  and  Riehl  (1960)  for  a  moderate 

2 

hurricane  was  2420  cal/cm  -day.  Leipper  acknowledges  that  the  transfer 
so  predicted  is  so  much  larger  to  the  rignt  of  the  path  that  even  observed 
asymmetries  in  intensity  are  inadequate  for  his  estimation.  Further 
reasons  for  doubt  concerning  Leipper's  calculation  "or  Hurricane  Hilda  will 
be  given  in  Section  3.9.1. 

Leipper  and  Volgenau  (1972)  later  introduced  the  concept  of  a 
hurricane  heat  potential,  which  they  applied  to  oceanographic  data  for  the 
Gulf  of  Mexico.  The  idea  is  that  all  the  heat  content  of  the  upper  oceanic 
layer  (top  500  cm)  implied  by  water  temperature  in  excess  of  26°C  is 
available  to  sustain  a  hurricane.  The  authors  acknowledge  that  their 
theory  implies  that  in  a  hurricane,  the  air  controls  the  sea  temperature, 
although  in  almost  every  other  physical  situation,  the  sea  is  regarded  as 
characterizing  the  air  temperature.  If  all  this  heat  were  really  available 
for  transfer,  why  would  a  hurricane  be  required  (rather  than  just  the  nor¬ 
mally  operative  turbulent  transfer  mechanisms)  to  effect  the  transfer? 

One  wonders  how  such  a  large  potential  could  be  developed,  and  why  hurricanes 
are  not  always  formed  over  waters  with  surface-layer  temperatures  above 
26°  C. 
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3 •  9  Critique  of  the  Riehl -Mai kus  Model 


A  critique  of  the  Riehl-Malkus  reasoning  that  an  extra  oceanic 

source  of  latent  and  sensible  heat  maintains  the  tropical  cyclone  might 

y/(y-1 ) 

begin  by  noting  that  use  of  the  adiabatic-expansion  relation  p  ^  T 
is  inappropriate  for  a  (turbulent)  frictional  layer.  Such  a  relation  implies 
that  the  entropy  is  constant  for  each  fluid  element.  Riehl  (1954,  p.  26) 
is  correct  that  neglect  of  diffusion  gives  almost  a  6°C  drop  (with  conden¬ 
sation)  in  passing  from  1014  mb  (autumnal  sea-level  ambient)  to  930  mb. 

In  fact,  the  inflow  layer  is  largely  cloud-free,  and  the  same  adiabatic 
expansion  taken  as  dry  produces  a  7.3°C  drop.  Such  a  large  temperature  drop 
does  not  occur  (though  whether  the  static  temperature  is  really  constant 
as  Riehl  believes  is  uncertain),  and  should  never  have  been  expected. 

But  Riehl  (1954.  p.  286)  was  surprised  at  measurements  he  interpreted  to 
imply  constant  static  temperature  during  inflow  in  the  frictional  layer, 
and  passed  immediately  from  the  (clearly  inappropriate)  isentropic  expan¬ 
sion  to  the  postulate  of  latent  and  sensible  heat  transfer  from  the  ocean 

9 

greatly  in  excess  of  ambient  transfer  level:. 

3.9.1  Evaluating  the  Arguments  for  Greatly  Augmented  Enthalpy  Transfer 

There  are  no  indisputable  semi-empirical  theories  or  direct  measure¬ 
ments  of  the  heat  and  mass  transfer  from  the  ocean  to  the  atmosphere  in 
a  hurricane:  "unfortunately  there  is  little  information  on  C^-  (  the 
nondimensional  coefficient  for  the  air-sea  energy  interchange  )  under 
hurricane  conditions,  other  than  the  semispeculative  guess  that  the  exchange 
coefficients  of  latent  heat,  sensible  heat  and  momentum  are  probably  of 
the  same  magnitude"  (Ooyama,  1969,  p.  15).  Hidy  (1972,  pp.  1086,  1097) 


88 


notes  the  lack  of  quantitative  progress  on  air-sea  transfer,  cites  the 
paucity  of  data  on  the  drag  coefficient  at  sea  level  for  winds  over  20  mph, 
and  mentions  the  difficulties  in  extrapolating  data  accurately;  in  fact, 
one  of  few  known  results  casts  doubt  on  the  validity  of  the  Reynolds 
analogy  relating  turbulent  transfer  of  heat  and  momentum  at  the  air-sea 
interface.  There  have  been  attenr  ts  to  extrapolate  to  hurricane  conditions 
empirical  laws  relating  ocean  evaporation  rates  to  wind  speed.  The  actual 
measurements  refer  to  low  wind  speed  and  involve  relatively  dry  air,  far 
from  saturation;  extrapolation  to  hurricane  winds  and  nearly  saturated 
low-level  tropical  air  is  unjustified  and  misleading.  Garstang  (1967) 
emphasizes  that  a  critical  assumption  in  the  theory  of  Riehl  and  Malkus 
is  that  the  sea-to-air  transfer  coefficients  increase  linearly  with  wind 
speed  indefinitely  into  the  hurricane  regime.  Zipser  (1969,  p.  813)  found 
that  using  such  linearly-extrapolateci  hulk  turbulent  transfer  relations, 
which  imply  an  order-of-magnitude  increase  in  latent  heat  transfer  under 
moderately  strong  tropical  disturbances  over  the  ambient  rate  of  transfer, 
would  have  produced  clouds  in  an  observed  tropical  disturbance  with 
unsaturated  convective  downdrafts  on  a  time  scale  an  order-of-magnitude 
too  short.  Zipser  also  notes  that  Riehl  and  Malkus  expected  downdrafts 
around  cumulonimbi  to  be  saturated,  but  the  observed  downdrafts  were 
unsaturated  and  so  details  of  the  Riehl-Malkus  predictions  on  sea-to-air 
transfers  near  cumulonimbi  must  be  modified.  In  the  absence  of  more 
quantitative  evidence  Riehl  (1954,  p.  287)  asserts  that  much  sea  spray  is 
tossed  into  the  air  in  a  hurricane.  This  does  not  in  itself  assure  augmented 
net  enthalpy  transfer  from  sea  to  air,  since  heat  must  be  drawn  from  the 
air  to  evaporate  the  drops  for  later  condensation  of  the  water  vanor  in 
the  eyewall  and  inner  rainbands.  evaporation  of  spray  leads  to  a  temperature 


decrease  as  well  as  a  dew-point  increase,  such  that  the  total  stagnation 
enthalpy  remains  unchanged. 

Possibly  the  closest  approximation  to  a  direct  measurement  of  the 
Riehl-Malkus  postulate  of  an  augmented  oceanic  heat  source  in  a  hurricane 
are  the  tritium-tracing  measurements  by  Ostlund  in  Hurricane  Faith  (1966). 

The  tritium  deposited  initially  in  the  stratosphere  following  American 
and  Russian  fusion  bomb  tests  in  1961  and  1962  temporarily  created  trace 
amounts  in  the  troposphere  that  exceeded  the  amounts  in  the  ocean;  deter¬ 
mining  tritium  concentrations  could  theoretically  determine  the  source 
of  water  substance  in  a  hurricane.  While  the  experiment  cannot  be  repeated 
without  the  resumption  of  atmospheric  nuclear  testing,  and  while  the  present 
author  does  not  regard  the  experiment  as  highly  reliable,  Hidy  (1972) 
reported  the  result  without  reservations .  Ostlund  found  that  the  ambient 

^  rate  of  latent  heat  transfer  from  sea  to  atmosphere  in  the  tropics  (which 

is  equivalent  to  evaporation  of  about  0.5  cm/day  of  water)  was  increased  by 
a  factor  of  about  three  within  the  inner  100  km  radius  of  the  hurricane, 
as  opposed  to  the  order-of -magnitude  increase  predicted  by  the  theory  of 
Riehl  and  Malkus. 

There  seems  to  be  no  indisputable  direct  evidence  for  the  internal 
ocean  heat  source  postulate  of  the  Riehl-Malkus  theory;  the  question  does 
remain  whether  or  not  Carrier  and  his  co-workers  have  obviated  the  ind-rect 
justification,  specifically,  the  assertion  by  Riehl  and  Malkus  that  the 
oceanic  heat  source  is  necessary  to  satisfy  the  laws  describing  conservation 
of  energy  and  momentum  iri  the  storm. 

In  this  respect  the  earlier  quotation  from  Rosenthal  (1971b,  p.  772), 
citing  Ooyama  (1969,  pp.  14-15),  should  be  discussed.  Rosenthal  asserts 
that  higher-level  air  sinking  down  into  the  surface  frictional  layer  will 
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be  drier  and  of  lower  equivalent  potential  temperature  than  the  original 
boundary-layer  air.  If  the  quasisteady  approximation  of  the  Carrier  model 
is  justified,  then  in  the  mature  stage  the  total  stagnation  enthalpy  is 
a  function  mainly  of  position.  According  to  Carrier  and  his  co-workers, 
the  plume  and  radiative  transfer  mechanisms  that  carry  heat  and  moisture 
up  into  the  ambient  tropical  atmosphere  persist,  with  neither  augmentation 
nor  diminution  to  any  appreciable  degree,  in  the  presence  of  the  storm. 

These  transfer  mechanisms  help  compensate  for  rain-out  in  the  outer  spiral 
bands  such  that  by  the  time  a  particle  has  slowly  sunk  into  the  boundary 
layer  at  0.005  mph,  it  has  been  enriched  in  total  stagnation  enthalpy  to 
the  ambient  value  at  its  current  height.  In  other  words,  much  of  the  sea- 
to-air  transfer  of  latent  and  sensible  heat  continues  to  pass  across  the 
boundary  layer  with  little  diminution,  just  as  in  the  ambient,  such  that 
the  enthalpy  of  air  originally  in  the  700-900  mb  strata  increases  to  that  of  air 
originally  in  the  boundary  layer  owing  to  enrichment  of  its  total  static 
enthalpy  as  it  slowly  descends.  Only  as  the  hurricane  leaves  the  tropical 
oceans  is  this  ambient  transfer  reduced;  the  air  entering  the  boundary  layer 
eventually  is  of  lower  total  static  enthalpy,  since  it  comes  from  air  which 
is  originally  higher  in  the  tropical  ambient  (hence  colder  and  drier),  and 
since  that  air  is  not  appreciably  enriched  as  it  descends.  In  this  way 
traverse  over  ocean  patches  of  varying  temperature  can  heln  cause  the  well- 
known  nonmonotonic  perturbations  in  hurricane  intensity  within  the  general 
level  of  strength  computed  from  the  spawning  ambient  as  discussed  earlier 
(see  also  Section  4.4). 

In  this  critique  of  the  Riehl-Malkus  model,  one  should  note  a  later 
modification  of  earlier  statements.  Riehl  and  Malkus  (1960,  p.  17)  write: 

"The  total  added  heat  energy  from  the  ocean  is  on  the  order  (for  the  moderate 
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storm)  of  2.5  cal/gm,  while  the  average  normal  heat  content  (latent  dIus 
sensible)  of  tropical  air  is  about  80  cal/gm."  Thus  the  enthalpy  contribution 
added  by  the  ocean  to  ambient-air  enthalpy  is  about  three  percent  for  a 
moderate  storm.  "The  heat.  .  .  gained  by  the  air  is  only  a  minute  fraction 
of  that  carried  inward  through  the  cylinder  at  the  distance  of  the  1000~mb 
isobar.  .  ..  Nevertheless,  it  is  these  small  increments  that  produce  the 
strong  inward  warming.  .  ..  They  are  thus  of  utmost  importance  for  genera¬ 
tion  and  maintenance  of  tropical  storms,  even  though  they  may  be  wholly 
neglected  in  a  general  energy  balance.  .  .  "  (Palmen  and  Riehl,  1957, 
p.  156).  That  such  a  marginal  amount  of  heating  is  critical  to  both  genera¬ 
tion  and  maintenance  seems  intuitively  surprising.  Later  Malkus  (1962, 
p.  247-249)  limited  the  role  of  greatly  augmented  enthalpy  transfer  strictly 
to  the  maintenance  of  the  hurricane,  since  during  development  the  high 
winds  that  supposedly  permit  the  increased  turbulent  transfer  do  not  exist. 

Hawkins  and  Rubsam  (1968)  have  computed  the  structure  and  budgets 
for  Hurricane  Hilda  for  October  1,  1964,  when  it  was  an  intense  storm 
over  the  Gulf  of  Mexico.  An  unprecedented  five-level  collection  of  data 
was  available  for  exhaustive  analysis,  and  the  total  static  enthalpy  budget 
was  calculated  very  similarly  to  the  approach  used  by  Malkus  and  Riehl 
(1960).  Hawkins  and  Rubsam  found,  however,  that  (upon  inclusion  of  the 
postulated  latent  and  sensible  enthalpy  flux  from  the  sea)  the  amount 
of  energy  estimated  to  be  radiated  away  from  the  top  of  the  storm  was 
twenty  times  what  seems  physically  possible.  The  kinetic  energy  budget 
computed  for  Hilda  was  also  at  significant  variance  f^om  the  one  computed 
by  Riehl  and  Malkus  (1961)  for  Hurricane  Daisy  (1958).  In  view  of  this, 
Leipper's  (1967)  suggestion  that  the  oceanic  heat  transfer  to  Hilda  would 
be  underestimated  using  Malkus's  formulae,  by  50-90%  seems  incomprehensible. 
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3.9.2  Temperature  Measurements  in  Hurricanes 


If  numerous  accurate  temperature  measurements  in  rapidly  swirling, 
heavily  raining  porticns  of  the  hurricane  could  be  readily  carried  out, 
resolution  of  doubt  about  the  existence  of  a  large  ocean  heat  source  within 
hurricanes  could  be  quickly  accomplished.  However,  such  measurements  are 
still  subject  to  controversy  and  inaccuracy.  Aircraft  measurements  tend 
to  be  taken  above  the  inflow  layer  for  safety  reasons:  "...  the  preferred 
altitude  for  reconnaissance  traverses  of  a  hurricane  is  from  10,000  to 
14,000  feet"  (Meyer,  1971a,  p.  58).  Currently  used  instrumentation  may 
suffer  from  slow  response  time,  from  wet-bulb  effects,  and  from  the  need 
to  deduce  the  static  temperature  from  a  measurement  that  includes  the 
dynamic  contribution  (Meyer,  1971b,  pp.  19-37).  Riehl  and  Malkus  (1961, 
p.  188)  discuss  possible  unresolved  errors  in  midtropospheric  temperature 
measurements  of  a  few  degrees  Centigrade  in  magnitude. 

First,  to  a  reasonable  degree  of  accuracy  a  dry-bulb  thermometer 
measures  [T  +  (q2/2  c  )],  and  a  wet-bulb  thermometer  (with  knowledge 
of  the  pressure)  yields  [I  +  (q2/2  c  )  +  (LY/cp)]  (Hess,  1959,  p.  61). 

At  very  low  altitudes,  the  dry-bulb  temperature  is  the  total  potential - 
like  temperature  (e./'c  ),  and  the  wet-bulb  temperature  is  then  related  to 

L  p 

the  total  stagnation  temperature  (Ht/cp),  because  the  miss>ng  gravita¬ 
tional  contribution  (gz/c  )  is  negligible.  Hert  the  recovery  factor 
(usually  denoted  r)  for  the  thermometers  is  taken  as  unity;  for  the 
relevant  Mach  and  Reynolds  numbers  this  seems  satisfactory,  since  for 
a  well -designed  instrument  r  >  0.9  for  turbulent  flow  in  air.  The 
belaboring  of  the  dynamic  contribution  would  normally  be  unjustified  in 
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meteorology,  but  a  200-mph  wind  is  equivalent  to  a  contribution  of  about 
4°C;  this  amount  is  in  excess  of  one  percent  of  both  the  static  temperature 
(•^300°K)  and  also  the  total  stagnation  temperature  (^3bO°K)  for  the  sea- 
level  autumnal  tropical  ambient.  Whether  temperature  measurements  made 
several  decades  ago,  when  hurricanes  generally  were  not  believed  to  be 
capable  of  such  speeds,  are  properly  corrected  for  this  dynamic  contribution 
is  not  always  clear.  It  should  also  be  recalled  that,  while  the  theory  needs 
some  improvement,  the  still-quite-useful  result  of  Carrier,  Hammond,  and 
George  (1971)  implies  that  the  wet-bulb  temperature  would  remain  approxi¬ 
mately  radially  invariant  in  a  traverse  of  a  hurricane  along  a  ray  and  at 
constant  height  up  to  about  700  mb  altitude,  from  the  outer  eyewall  to  the 
outer  edge. 

Deppermann  (1937,  p.  5),  whose  work  is  frequently  cited  by  Riehl, 

writes:  ".  .  .  remarkably  uniform  diurnal  temperature  oscillations  are 

maintained  both  before  and  after  a  typhoon;  but  when  the  station  is  under 

typhoon  influence,  a  reduction  of  the  maximum  temperature  by  about  3°C 

occurs,  while  the  minimum  temperature  remains  the  same."  If.  by  maximum 

and  minimum  temperatures  Deppermann  means  dry-bulb  and  wet-bulb  temperatures, 

respectively,  and  if  one  assumes  Deppermann  correctly  accounted  for  the 

o 

dynamic  contribution  (q  /2  c  ),  then  Deppermann  finds  that  the  static 
temperature  T  decreases  by  about  3°K  in  a  typical  typhoon,  while  the  total 
stagnation  temperature  increases  by  a  few  degrees  Centigrade  (5°C-6°C  at 
very  most).  In  any  case,  there  is  no  evidence  of  large  increases  in  (H/cp),  as 
an  oceanic  heat  source  would  purportedly  cause.  Also,  if  (H/cp)  is  roughly 
radially  constant  in  the  inflow  layer,  then  a  decrease  in  T  of  about  3.7°C 
from  sea-level  autumnal  tropical  ambient  is  necessary  for  condensation,  and 
this  is  reportedly  not  achieved  in  a  typical  typhoon,  so  a  cloud-free 


frictional  inflow  layer  is  no  surprise  in  Carrier's  model. 

Arakawa  (1954,  p.  119)  reported  that  the  air  temperature  dropped 
1°C  from  28°C  to  27°C,  and  the  wet-bulb  temperature  remained  about  constant 
radially  at  approximately  2S.5°C,  during  the  passage  of  an  898-mb  typhoon 
over  a  Japanese  weather  ship  in  October  1944;  only  in  the  eye  did  the 
two  temperatures  rise  sharply.  Winds  rose  to  about  120  mph  in  the  eyewall . 
Again,  this  seems  compatible  with  the  Carrier-Hammond-George  predictions, 
which  include  no  oceanic  heat  transfer  within  hurricanes  much  above  ambient 
level.  Palmen  and  Newton  (1969,  p.  478)  assert  with  no  elaboration  that  the 
measurement  implies  an  increase  of  the  equivalent  potential  temperature 
from  360°K  to  385°K,  and  thus  is  a  confirmation  of  the  internal  heat  source 
postulate;  the  author  does  not  understand  how  this  result  was  achieved. 

To  avoid  difficulties  in  boundary  layer  measurements,  Riehl  (1963, 
p.  277)  studied  equivalent  potential  temperatures  deduced  from  aircraft 
data  taken  at  the  245-250  mb  level  at  the  inner  edge  of  the  eyewall.  Little 
complication  from  moisture  is  expected  at  such  heights,  and  the  variation 
of  equivalent  potential  temperature  with  height  is  not  expected  to  be 
large.  Riehl  reports  equivalent  potential  temperature  of  370°K  for  a 
hurricane  of  central  pressure  of  960  mb,  but  provides  no  further  details 
on  the  instrumentation  or  data  reduction  whatever.  If  the  result  is  cor¬ 
rect,  very  significant  evidence  exists  for  an  oceanic  heat  source.  However, 
Dr.  R.  H.  Simpson  of  the  National  Hurricane  Center  in  Miami,  Florida  has 
informed  the  author  that  a  vortex  thermometer  was  used  (private  correspon¬ 
dence).  If  one  recalls  that  250  mb  corresponds  very  roughly  to  approximately 
35,000  ft,  then  the  remarks  of  Gentry  (1964,  p.  12)  are  illuminating:  "The 
wind  tunnel  tests  of  the  vortex  probe  under  icing  conditions  (Ruskin  and 
Schecter.  .)  indicated  that  the  recorded  temperature  might  be  too  high 
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due  to  the  change  of  state  on  impact  at  probe  entry.  ...  No  evidence  is 
readily  available  as  to  whether  much  ice  accumulated  on  the  vortex  probe 
of  the  aircraft  flying  at  about  35,000  ft  or  higher.  Ordinarily,  little 
liquid  water  would  be  encountered  at  the  higher  elevations  and  the  ice 
accumulation  should  be  minor.  Some  data  collected  by  the  research  air¬ 
craft  in  1961  and  1962  indicated,  however,  that  at  times  the  ice  accumula¬ 
tion  might  be  significant.  There  is  the  possibility,  therefore,  that  the 
temperatures  recorded  for  the  upper  troposphere  are  too  high  by  an  unknown 
amount.  There  will  probably  always  remain  a  question  as  to  the  accuracy 
of  the  temperature  measurements  until  some  absolute  standard  of  comparison 
is  developed." 

3.9.3  On  the  Magnitude  of  Possible  Increases  in  the  Wet-Bulb  Temperature 
at  Sea  Level 

Some  insight  into  possible  increases  in  the  total  stagnation  tempera¬ 
ture  (H/Cp)  =  T  +  (gz/Cp)  +  ( q 2/ 2  cp)  +  (LY/cp)  at  sea  level  within  a 
hurricane  will  now  be  sought.  Very  near  z  =  0,  q  =  0.  Also,  the  air 
temperature  assumes  the  temperature  of  the  sea;  since  the  sea  temperature 
is  held  spatially  constant  in  all  hurricane  analyses  known  to  the  author 
[with  one  experimental  exception  (Anthes,  1972,  pp.  473-474)],  and  since 
the  sea  temperature  for  an  intense  slowly  translating  hurricane  decreases 
toward  the  center  of  the  storm  if  it  changes  at  all,  the  conservative 
aoproximation  for  current  purposes  is  to  hold  T  everywhere  constant  at  the 
ambient  sea-level  value  for  the  tropical  autumn.  Then  the  only  changes  in 
(H/Cp)  are  due  to  changes  in  ( LY/Cp) ,  or  more  specifically,  to  changes  in 
Y  since  L  and  cp  are  taken  constant  for  current  purposes.  Further,  with 
only  one-percent  error, 
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o(RH)  P(T)  •  a(RH)  P(T) 
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so  the  vapor  mass  fraction  Y  increases  as  the  relative  humidity  RH  increases 
or  as  the  total  pressure  p  decreases,  for  temperature  T  fixed.  In  the 
autumnal  tropical  ambient,  if  L  represents  latent  heat  of  condensation 
plus  glaciation,  at  sea-level  (LY/Cp)  =  50CIK.  The  ambient  sea-level  relative 
humidity  is  0.84  and  the  pressure  is  1014  mb.  [The  relative  humidity  is 
unity  at  the  air-sea  interface  but  sea-level  measurement  generally  means 
readings  at  shipboard  height,  or  more  phenomenologically,  it  means  measure¬ 
ment  in  the  postulated,  roughly  thirty-foot-thick  constant-stress  layer 
near  the  air-sea  interface  (Hidy,  1972,  p.  1084).]  For  a  moderate  hurricane 
with  central  pressure  of  966  mb,  the  increase  in  (H/c  )  is  2.5°K,  and 

r 

for  an  intense  hurricane  with  central  pressure  of  910  mb,  the  increase  in 

(H/Cp)  is  5.7°K  —  provided  RH  is  held  at  0.84.  If  RH  increases  from  0.84 

to  unity  under  the  eyewall,  then  (H/c  )  increases  by  12.5°K  for  the  moderate 

P 

hurricane  and  16.3°K  for  an  intense  hurricane  --  although  the  justification 
for  adopting  such  an  increase  in  RH  appears  unestablished,  and  Anthes 
and  Johnson  (1968,  p.  297)  cite  evidence  that  the  relative  humidity  of  the 
boundary  layer  tends  to  remain  between  80°/  and  90%.  Malkus  and  Riehl  (1960, 
p.  16)  claim  the  equivalent  potential  temperature,  which  is  clearly  closely 
related  to  (H/c  ),  increases  by  12.5°K  from  ambient  to  eyewall  in  the 

r 

moderate  hurricane  [although  Riehl  (1963,  p.  277)  seems  to  augment  this  to 
16°K],  and  by  35°K  in  the  intense  hurricane.  By  the  reasoning  adopted 
here  the  physical  basis  of  the  12.5°K  increase  becomes  evident,  but  the 
35°K  increase  for  an  intense  hurricane  represents  a  physically  impossible 
relative  humidity  of  over  100%.  Succinctly,  while  Carrier  and  his  co-workers 


discuss  increases  in  the  first  normal  derivative  oF  H  at  sea  level  in 
hurricanes  over  ambient  on  the  order  of  10%,  Riehl  and  Malkus  discuss 
occasional  increases  in  H  itself  on  the  order  of  10%.  Carrier  and  his 
r.o-wcrkers  neglected  the  full  effect  of  decreased  air  density  as  the  eyewall 
is  approached,  and  found  H  to  be  radially  constant;  the  conjecture  here  is 
that  the  error  in  H  so  incurred  is  about  one  percent  or  so. 

3.10  Numerical  Simulation  of  Hurricanes  on  Digital  Computers 

In  recent  years  there  has  been  a  proliferation  of  computer  models 
of  the  tropical  cyclones  (Yamasaki,  1968;  Ooyama,  1969;  Rosenthal,  1970; 
Sundqvist ,  1970;  Kurihara,  1971;  Anthes,  1972).  These  models  warrant  a 
review  in  themselves,  and  only  a  limited  discussion  is  undertaken  here. 

First,  it  must  be  understood  that  these  models  are  by  no  means 
exact  solutions  of  the  full  conservation  equations,  subject  to  appropriate 
boundary  and  initial  conditions.  In  fact,  to  paraphrase  Lorenz  M967), 
one  might  remark  that  detailed  reproduction  of  the  hurricane  lifespan 
by  numerical  methods  might  not  appreciably  increase  physical  insight,  because 
the  total  behavior  is  so  complex  that  the  relative  impor Lance  of  various 
features  might  be  no  more  evident  from  examination  of  the  computer  solutions 
than  from  direct  observations  of  the  real  atmosphere.  In  any  case,  com¬ 
puter  solutions  are  today  far  from  exact  solutions;  in  fact,  today  no 
computer  solution  even  attempts  to  treat  observed  initial  data  (Ooyama, 

1969,  pp.  35-37;  Rosenthal,  1971b,  p.  35),  but  rather  describes  how  a 
general  hurricane  grows  and  dies  --  in  fact,  works  in  computer  simulation 
emphasize  results  for  the  later  stages  of  intensification  and  the  peak- 
intensity  portion  of  the  lifetime.  The  adopted  initial  conditions  almost 
invariably  have  pronounced  circulations ,  and  only  major  alteration  of 


98 


parameterizations  prevent  inevitable  hurricane  formation. 

Mention  of  parameterizations  introduces  another  area  of  compromise 
with  exact  solution.  Because  many  phenomena  on  the  cumulus  scale  and 
smaller  are  crucial  to  the  hurricane-scale  flow,  computer  solution  today 
is  feasible  only  if  the  roles  of  turbulent  diffusion,  radiative  transfer, 
and  (especially)  cumulus  convection  are  simulated  by  parameterization.  No 
one  knows  with  certainty  how  to  describe  the  role  of  these  phenomena  in 
terms  of  variables  occurring  in  the  hurricane-scale  solution,  or  even  if 
adequate  description  of  this  type  is  possible.  Reviewing  the  parameteriza¬ 
tions  of  cumulus  convection  currently  used  in  numerical  simulations  of 
hurricanes,  Ooyana  (1971,  p.  744)  acknowledges  that  .  .  it  is  generally 
agreed  that  these  methods  are  essentially  stopgap  measures  which  are  used 
for  want  of  better  alternatives."  Gray  (1972,  p.  57,60)  writes:  "...  the 
authors  (Gray  and  Lopez)  feel  that  it  is  very  difficult  (or  impossible)  at 
the  present  time  for  any  numerical  model  of  tropical  motion  to  come  to 
grips  with  the  real  problems  of  incorporating  the  cumulus  convection  in 
terms  of  the  broader-scale  flow.  ...  In  that  the  cumulus  cloud  is  such  a 
distinctive  physical  unit,  it  would  appear  that  it  should  be  independently 
treated.  ...  The  physics  of  the  cumulus-broadscale  interaction  may  not 
necessarily  be  overcome  by  applying  the  primitive  equations  to  even  smaller 
grids  and  time  steps  without  additional  insight  into  the  character  of  the 
individual  convective  elements."  While  there  are  some  physical  concepts 
that  guide  currently  employed  parameterizations  (such  as  weak  integral 
constraints  on  energy  and  water  substance),  the  parameterizations  remain 
rather  arbitrary  and  experimental  (for  a  summary  of  them,  see  Haltiner, 

1971;  Bates,  1972).  Often  the  parameterizations  are  freely  adjusted  until 
the  cyclone-scale  results  are  consistent  with  measurements  and/or  intuition. 
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The  computer  results  are  very  sensitive  to  the  assumptions  made  concerning 
the  formulation  of  cumulus  convection  and  eddy  transfer  (Ramage,  1971, 
p.  40;  Haltiner,  1971,  p.  262;  Garstang,  1972,  p.  619).  Without  any 
disparagement,  such  adjustment  of  parameterizations  may  be  termed  curve- 
fitting.  What  is  learned  from  current  computer  experiments  is  mute.  Of 
course,  no  model  can  today  entirely  escape  such  curve-fitting.  The  ouestion 
to  be  pondered  is  how  refined  a  solution  is  currently  warranted  in  view 
of  the  limited  physical  understanding  of  certain  critical  phenomena? 

A  further  complication  in  evaluating  numerical  models  is  that  some¬ 
times  the  finite  differencing,  intentionally  or  inadvertantly,  adds  to  the 
stated  formulation  in  the  sense  that  effects  not  in  the  differential  equa¬ 
tions  are  present  in  the  difference  equations.  The  most  accurate  differencing 
of  the  differential  equations  is  not  necessarily  preferred,  but  rather 
numerical  techniques  are  rated  ".  .  .on  intuitive  meteorological  inspection 
of  results"  (Anthes,  Rosenthal,  and  Trout,  1971,  p.  747).  For  instance, 
errors  introduced  by  less  accurate  upstream-difference  schemes  were  retained 
to  simulate  lateral  mixing  when  more  accurate  central -difference  schemes 
gave  less  realistic  results  (Rosenthal,  1970,  pp.  657-658).  Later  the 
differential  formulation  was  modified  so  that  an  accurate  differencing  of 
the  modified  expression  for  lateral  diffusion  gave  the  numerical  contribu¬ 
tion  desired.  However,  the  modified  expression  ".  .  .  is  not  very  satisfying 
from  a  physical  point  of  view"  (Anthes,  Rosenthal,  and  Trout,  1971,  p.  747). 
What  can  be  learned  from  such  procedures  aqain  seems  mute. 

Even  if  the  parameterizations  for  cumulus  convection,  turbulent 
mixing,  and  radiative  transfer  were  adequately  known,  and  the  finite- 
differencing  were  adequately  accurate,  the  question  of  feasibility  in 
seeking  a  uniformly  valid  solution  for  the  entire  cyclone  would  remain. 
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Large  gradients  occur  over  relatively  small  scales  in  important  subregions 
of  the  storm  (e.g.,  the  frictional  boundary  layer  and  the  eyewall),  while 
small  gradients  occur  over  relatively  large  scales  in  the  bulk  of  the 
storm  (the  rapidly  swirling  regions  and  the  outflow  layer  aloft).  On 
current  computers  usually  a  more-or-less  fixed  grid  of  10  km  or  20  km 
radial  resolution  and  at  most  thirteen  layers  of  vertical  resolution  is 
adopted  by  practical  considerations;  the  overwhelming  bulk  of  the  grid 
points  then  lie  outside  the  eye,  the  eyewall,  and  the  frictional  boundary 
layer  --  where  important  processes  are  occurring.  The  domain  sizes  are 
quite  limited  (typically  440  km)  so  relatively  weak  boundary  conditions 
(requiring  only  that  purely  advective  influx  occur  at  the  side  boundaries) 
are  employed  and  much  of  the  storm  lies  outside  the  domain  of  computation. 
Garstang  (1972,  p.  619)  states  numerical  results  are  often  very  sensitive 
to  the  treatment  of  the  outer  limits.  An  alternative  to  direct  numerical 
solution  of  the  uniformly  valid  equations  is  to  subdivide  the  storm  into 
natural  portions  where  different  gradients  and  phenomena  are  operative. 

This  is  the  approach,  with  provision  for  interfacial  compatibility,  used 
by  Carrier  and  his  co-workers. 

Some  evidence  for  uncertainty  concerning  current  numerical  models 
for  the  hurricane  is  provided  by  results  for  the  frictional  inflow  layer. 
First,  with  his  model,  Rosenthal  (1971b)  performs  numerical  experiments  in 
which  the  latent  heat  transfer  from  the  sea  is  greatly  augmented  over  ambient 
but  the  sensible  heat  transfer  is  set  to  zero.  The  goal  is  to  show  that 
the  latent  heat  transfer  is  the  important  portion  of  the  total  enthalpy 
transfer.  Of  course,  Rosenthal  is  well  aware  that  the  sensible  heat  trans¬ 
fer  is  required  to  permit  the  latent  heat  transfer  (Riehl,  1954,  p.  336). 

The  computation  is  cited  as  evidence  of  how  incompletely  the  full  conservation 
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equations  are  represented  in  current  models.  Thus  the  current  computer 
programs  can  not  prove  the  physical  validity  of  the  Riehl-Malkus  oceanic- 
heat-source  postulate;  the  solutions  which  the  current  computer  programs 
generate  merely  reflect  the  incorporation  of  the  Riehl-Malkus  postulate 
in  their  formulation.^  Next,  the  study  by  Anthes  (1971)  of  the  incomnressibl 
ax^symmetric  flow  in  a  frictional  layer  under  an  intense  Rankine-like 
vortex  furnishes  an  opportunity  to  test  how  well  the  finite-difference 
techniques  used  in  computer  simulations  of  hurricanes  succeed  in  a  problem 
of  reasonable  complication  to  which  the  answer  is  apparently  known. 
Interestingly,  Anthes  introduces  an  initial-value  problem  as  a  convenient 
device  for  achieving  the  desired  steady  solution,  without  claiming  validity 
for  the  transient  phase  necessarily  computed.  There  are  three  points 
about  the  solution  that  will  be  noted.  (1)  The  multilayer  solution  reveals 
that  the  one-layer  treatments  of  Ooyama,  Rosenthal,  and  Yamasaki  under¬ 
estimate  the  tangential  and  (particularly)  the  radial  velocity  components 
of  the  frictional  layer,  and  give  a  premature  eruption  from  the  Ekman  laver 
at  incorrectly  large  radial  distances  from  the  center.  (2)  While  both 
diffusion  coefficients  are  held  constant,  Anthes  takes  the  radial  diffusion 
coefficient  four  orders  of  magnitude  larger  than  the  axial.  Either  by 
this  means,  or  by  large  truncation  errors  through  forward-differencing  of 
advective  terms,  for  plausible  solution  radial  diffusion  plays  a  role  in 
regions  of  the  flow  where  boundary-layer  theory  suggests  radial  diffusion 
should  be  negligible.  (3)  The  Ekman  layer  divergence  changes  from  a  weak 
downdraft  to  a  large  updraft  about.  125  km  from  the  axis  of  symmetry.  The 
spuriousness  of  this  result,  and  the  actual  two-part  structure  of  the 
frictional  layer  at  distances  from  the  axis  at  which  nonlinear  terms  are 
significant  [as  revealed  by  analytic  and  more  sophisticated  numerical 
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techniques  (Carrier,  Hammond,  and  George,  1971;  Carrier,  1971a;  Burggraf, 
Stewartson,  and  Belcher,  1971;  Dergarabedian  and  Fendell,  1972b)]  remain 
obscured  (Section  3.5). 

Nevertheless,  more  elaborate  computer  simulations  are  continually 
developed.  Anthes  (1972)  has  attempted  to  account  for  departures  from 
azimuthal  symmetry.  The  axis  of  the  simulated  hurricane  vortex  undergoes 
a  curious  anticyclonic  rotation  at  approximately  eight  mph  about  the  center 
of  the  storm.  Anthes  suggests  the  path  of  a  cyclone  is  not  influenced  by 
horizontal  sea-surface  thermal  gradients,  although  Brand  (1971)  claims 
observational  evidence  to  the  contrary. 

3-H  Implications  of  Hurricane  Models  on  Seeding 

The  current  mode  of  attempting  artificial  modification  of  tropical 
cyclones  is  to  introduce  silver  iodide  crystals  in  the  supercooled  water 
believed  to  exist  in  appreciable,  naturally  unfrozen  quantities  high  in 
the  eyewall,  and  perhaps  also  in  rainbands  and  cumulus  clouds  both  close 
to  and  far  from  the  eyewall.  In  a  typical  field  experiment,  over  two 
hundred  rocket  cannisters,  each  containing  about  150  g  of  silver  iodide, 
are  released  at  about  33,000  feet  over  a  30  km  distance  (entirely  outside 
the  eye),  along  a  ray  90°  clockwise  to  the  direction  of  storm  translation; 
such  an  experiment  takes  two-to-three  minutes  to  execute.11  In  laboratory 
studies,  a  gram  of  silver  iodide  produces  1012  to  1014  crystals,  which 
serve  as  ice  nuclei  and  which  grow  to  three-micron  diameter.  The  goal  is 
that  there  will  be  10  nuclei/liter-of-space  over  a  depth  of  several 
kilometers  a  few  minutes  later;  this  concentration  is  reduced  to  negligible 
levels  in  a  couple  of  hours  (Gentry,  1971a;  Penner,  1972).  The  silver 
iodide  will  hopefully  cause  the  supercooled  water  droplets  to  freeze,  and 
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thus  to  release  the  heat  of  fusion  (Battan,  1969). 

The  proposed  mechanism  by  which  this  heat  release  (and  supposed 
attendant  temperature  rise  of  a  couple  of  degrees  Celsius,  and  density 
decrease)  is  efficacious  has  been  altered  on  several  occasions  (Simpson 
and  Malkus,  1964;  Rosenthal,  1971a;  Gentry  and  Hawkins,  1971),  and  at  times 
several  alternative  possible  mechanisms  have  been  set  forth  (Gentry,  1971a). 
A  concise  history  of  the  various  suggestions  has  been  given  by  Riehl  (1972a, 
pp.  267-270).  One  accepted  point  is  that  seeding  in  the  nascent  eye  of  a 
developing  tropical  storm  should  be  avoided  since  this  procedure  would 
probably  abet  intensification  (Rosenthal,  1971a).  The  most  recent  proposal 
has  been  to  seed  cumuli  outside  the  eyewall  and  above  the  0°C  isotherm 
(roughly,  600  mb  or  15,000  ft)  with  silver  iodide.  The  aim  is  not  to 
capitalize  directly  upon  the  heat  of  fusion  released  by  supercooled  water 
droplet  freezing,  but  rather  to  cause  these  cumuli  uufside  the  eyewall 
to  grow  several  thousand  feet  in  height  from  midtropospheric  to  upper 
tropospheric  levels.  This  would  hopefully  initiate  a  chain  of  events 
from  which  an  order-of-magnitude  more  heat  than  that  associated  with 
freezing  will  indirectly  be  realized.  This  added  heat  would  come  from 
some  low-level  inflowing  warm  moist  air  erupting  not  up  the  eyewall,  but 
prematurely  up  the  growing  cumulus,  with  associated  condensation  and 
additional  glaciation.  If  entirely  successful,  proponents  have  suggested 
that  perhaps  a  substitute  eyewall  further  from  the  center,  and  hence  (from 
conservation  of  angular  momentum)  with  reduced  swirling  speeds,  could  be 
realized  from  this  premature  eruption  of  boundary-layer  air  (Gentry  and 
Hawkins,  1971;  Rosenthal,  1971c). 

Fukuta  (1972)  suggests  that  the  field  experiments  conducted  to  date 
have  employed  overseeding  of  supercooled  water  droplets  to  form  small  ice 


104 


particles;  the  energy  release  would  be  small  and  the  effect  over  in  less 
than  half  an  hour.  Instead,  Fukuta  suggests  reducing  the  seeding  material 
to  1%  or  even  0.1%  of  that  used  in  past  field  experiments,  to  form  large 
precipitable  ice  particles.  The  latent  heat  of  fusion  these  falling 
particles  could  extract  from  the  inflow-updraft  stream  would  supposedly 
raise  the  central  pressure  for  hours,  reduce  swirling,  and  induce  snreadinq 
out  of  the  storm.  Detailed  analysis  fully  substantiating  either  the  con¬ 
cepts  of  Gentry  and  Hawkins  or  the  concepts  of  Fukuta  has  not  appeared. 

Clearly,  like  numerical  simulation  of  hurricanes,  the  seeding  of 
hurricanes  warrants  an  entire  review  unto  itself.  If  fact,  with  the 
proposed  extension  of  the  U.S.  seeding  effort  [Project  Stormfury,  a  joint 
effort  of  the  Department  of  Commerce  (NOAA)  and  the  Department  or  Defense 
(Navy)]  from  the  western  North  Atlantic  to  the  western  North  Pacific  Oceans 
(Mallinger,  1971),  this  might  be  an  apt  time  for  such  a  review. 

Attention  here  will  be  confined  to  a  few  observations,  especially 
what  the  hurricane  models  suggest  about  seeding  procedure. 

Even  proponents  of  seeding  anticipate  only  a  10  to  15%  decrease  in 
maximum  winds.  In  the  one  case  in  which  larger  decreases  were  noted, 
the  anomaly  is  now  attributed  to  synoptic  peculiarities  relating  to 
upper-level  outflow  (Hawkins,  1971).  This  is  not  to  disparage  such  decreases; 
statistical  treatment  of  a  model  suggests  multimillion  dollar  annual 
savings  in  damage  from  hurricane  seedings  that  would  reduce  peak  winds  by 
fifteen  percent  (Boyd,  Howard,  Matheson,  and  North,  1971).  The  problem 
is  that  the  storm  is  naturally  oscillating  in  intensity  by  the  same  order 
of  magnitude,  so  it  is  very  difficult  to  distinguish  natural  and  artificially 
induced  changes.  In  fact,  the  gradual  increase  of  the  maximum  winds  of 
Hurricane  Ginger  by  15%  during  the  day  after  initial  seeding  on  September 
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26,  1971  was  ascribed  to  natural  forces  by  Project  Stormfury  personnel, 
as  was  the  11%  decrease  after  the  second  seedinq  on  September  28  (Lieb, 
1972).  Fujita  ( 1972}  could  find  no  evidence  that  seeding  altered  Ginger; 
since  the  central  pressure  of  this  weak  hurricane  was  only  980  mb  at  the 
times  of  seeding,  he  suggests  that  perhaps  the  cloud-top  heights  were  too 
low  for  the  silver  iodide  to  be  effective.  The  threat  of  litigation  has 
constrained  the  number  of  seeding  experiments  severely. 

The  National  Hurricane  Research  Laboratory  (Gentry,  1969a;  Gentry, 
1969b;  Gentry,  1970)  has  discussed  a  six-to-twelve  hour  cycle  of  ameliora¬ 
tion  after  seeding;  physical  basis  for  this  time  scale  has  yet  to  come 
forth  (see  below).  Further,  if  the  central  pressure  deficit  is  reduced 
as  reported,  eyewall  seeding  must  alter  the  eye  in  a  still  unidentified 
manner,  although  Black,  Senn,  and  Courtright  (1972)  suggest  that  one  hour 
after  seeding  outside  the  eye,  the  eye  may  expand  more  than  is  consistent 
with  natural  variability.  More  complete  post-seeding  probing  of  the  tropi¬ 
cal  cyclone  v«ould  be  helpful  in  evaluating  these  claims. 

Rather  similar  computer  models  have  produced  different  guidance  with 
regard  to  current  seeding  practice.  The  reasons  for  the  discrepancy  are 
not  fully  available  because  some  details  remain  unpublished.  Rosenthal 
(1971c)  notes  that  both  the  magnitude  and  the  duration  of  the  heating  taken 
to  simulate  seeding  in  his  computer  experiments  seem  excessively  large. 
Nevertheless,  in  those  numerical  experiments  believed  to  most  closely 
simulate  the  field  experiments  performed  on  Hurricane  Debby,  August  18  and 
20,  1969,  Rosenthal  (1971a)  predicts  an  increase  after  seeding  in  the 
maximum  wind,  and  over  wide  radial  extents  an  augmentation  in  wind  level, 
at  700  mb;  the  comparable  measurements  at  12,000  ft  indicated  a  reduction 
in  winds  (Hawkins,  1971).  Rosenthal  predicts  a  decrease  in  the  maximum 
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wind  at  sea  level  after  seeding  (though  in  substantial  portions  of  the 
hurricane,  surface  winds  are  predicted  to  increase);  no  post-seeding  measure¬ 
ments  were  made  at  sea  level.  Quantitative  detail  is  omitted  here  because 
Rosenthal  (1971a,  p.  415)  notes:  .  .  at  best,  the  results  should  be 

considered  qualitative  guidance  material."  Sundqvist  (1971)  states  flatly 
that  seeding  will  intensify  a  hurricane,  Estoque  (1971,  p.  4),  noting  the 
tentative  nature  of  his  conclusions,  remarks: 

We  have  [studied]  the  properties  of  a  hurricane  model  which 
is  suitable  for  simulating  the  effects  of  artificial  seeding. 

The  main  improvement  of  this  model  over  previous  ones  is  the 
explicit  inclusion  of  cloud  micro-physical  processes,  including 
the  prediction  of  the  liquid  water  distribution.  Thus,  the 
model  is  able  to  simulate  more  realistically  the  release  of 
the  latent  heat  of  fusion  due  to  artificial  seeding.  The  model 
has  been  used  to  simulate  artificial  seeding  over  three  dif¬ 
ferent  radial  locations  in  the  vicinity  of  the  maximum  in 
the  surface  wind.  The  expected  reduction  and  outward  displace¬ 
ment  of  the  maximum  radial  surface  pressure  gradient  did  not 
materialize.  Instead,  the  pressure  gradient  and  also  the  wind 
increased  in  intensity.  However,  the  seeded  storm  intensity 
decreased  after  the  seeding  stopped.  The  greatest  effect  occurs 
when  seeding  is  done  over  a  region  just  inside  the  surface  wind 
maximum  location. 

An  interesting  feature  of  the  unseeded  hurricane  is  the 
occurrence  of  periodicity  in  time  of  the  hurricane  intensity. 

The  magnitude  and  the  period  of  this  periodicity  is  about 
10  m  sec’1  and  8  hours,  respectively.  This  magnitude  is  about 


the  same  order  as  that  of  changes  induced  artificially  by 
seeding.  If  this  periodicity  is  real,  it  should  be  taken 
into  account  in  the  interpretation  of  seeding  experiments  in 
the  actual  atmosphere. 

The  effectiveness  of  current  seeding  practices  remains  unresolved. 

If  the  Carrier  model  is  valid,  silver  iodide  seeding  can  only  transiently 
upset  the  stable  hurricane  configuration.  Under  this  model,  warm-fog  dis¬ 
persal  methods  would  have  to  be  applied  to  the  entire  "throughput  supply" 
layer  of  warm  moist  air  to  achieve  the  significant  goal  of  premature  rain- 
out  in  outer  spiral  bands.  The  layer  of  warm  moist  air  is  so  spatially 
extensive  that  such  attempts  seem  somewhat  impractical.  Dergarabedian  and 
Fendell  (1971)  have  discussed  use  of  warm-fog  dispersal  methods  for  lower 
levels  in  a  hurricane  for  premature  rainout  of  the  throughput-flux;  about 
simultaneously,  so  did  Mathews  (1971,  p.  H-8),  albeit  more  optimistically: 
The  feasibility  of  warm  cloud  modification  prior  to  cold 
cloud  modification  should  be  examined  because  warm  cloud 
modification  may  permit  growth  of  small  warm  clouds  to  tempera¬ 
tures  at  which  cold  cloud  modification  will  be  effective. 

The  combined  use  of  warm  cloud  and  cold  cloud  modification 
techniques  would  permit  selective  seeding  in  all  regions  with¬ 
out  cloud  top  temperature  restrictions. 

Seeding  to  divert  path  seemingly  holds  little  better  promise  than 

seeding  to  alleviate  intensity,  since  there  appears  to  be  no  way  to  discern 

what  path  alterations  were  due  to  human  intervention  under  current  under - 
1  ? 

standing. 
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4.  THEORY  OF  TROPICAL  CYCLONE  INTENSIFICATION 


There  exists  no  satisfactory  theory  explaining  how,  annually,  a 
small  number  of  tropical  disturbances  intensify  to  become  hurricanes,  and 
why  in  contrast  most  disturbances  do  not.  Charney  (1971)  notes  that, 
unlike  in  the  midlatitudes,  tropical  temperature  gradients  are  weak,  and 
tropical  cyclones  most  often  form  over  the  tropical  oceans  where  the  sur¬ 
face  temperatures  are  particularly  uniform.  No  theory  of  how  pronounced 
barocl inicity  develops  in  the  barotropic  tropics  is  likely  to  emerge 
until  the  detailed  local  balances  and  large-scale  circulations  of  the 
tropical  ambient  are  better  understood  quantitatively,  because  early 
stages  of  cyclogenesis  will  probably  involve  perturbation  about  that 
ambient. 

Carrier  (1971b)  has  presented  an  intensification  process  by  which 
his  quasisteady  mature  model  evolves  in  time  from  a  tropical  depression. 
Tracing  back  to  even  earlier  evolution  seems  premature  at  the  current 
state  of  understanding.  Although  Carrier's  theory  has  been  reproduced 
by  Penner  (1972)  without  objection,  Carrier's  theory  of  intensification 
is  not  currently  entirely  satisfactory.  Thus,  it  is  used  here  only  as  an 
interesting  vehicle  to  raise  some  significant  points  about  the  later, 
briefer  stage  of  intensification,  as  opposed  to  the  longer,  earlier  stage. 
Ooyama  (1969,  pp.  35-37),  amoig  many  others,  has  written  about  an  initial 
very  modest  rate  of  pressure  fall  in  intensifying  tropical  cyclones, 
followed  by  an  appreciably  accelerated  rate  of  fall  —  it  is  this  latter 
stage  which  is  principally  addressed  here. 
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4.1 


Carrier's  Outline  of  Intensif icat ion 


Carrier  begins  with  a  schematic  view  of  an  axisymmetric  tropical 
depression  in  which  there  is  a  weak  Rankine-vortex-like  swirling,  the 
maximum  azimuthal  speed  lying  relatively  far  from  the  axis  of  rotation. 

There  is  radial  inflow  in  region  I  during  the  transient  phase  (Fig.  19).  The 
swirling  would  quickly  establish  a  shear  layer  beneath  it;  there  is  weak 
upflux  out  of  (downflux  into)  the  surface  layer  where  the  swirling  speed 
increases  (decreases)  with  radial  distance.  Since  linear  theory  correctly 
predicts  the  downflux  in  the  mature  stage  with  rapid  swirl,  linear  theory 
probably  suffices  to  predict  the  boundary- layer  divergence  during  intensi¬ 
fication.  However,  since  equilibration  of  the  boundary  layer  requires 
times  of  or  roughly  sixteen  hours  in  the  tropics,  and  since  this 

is  possibly  an  appreciable  fraction  of  the  later-stage  rapid-intensification 
time-scale,  a  transient  linear  theory  will  be  required.  It  must  yield  a 
radial  influx  through  the  boundary  layer  II  in  excess  of  the  radial  inflow 
speed  in  I.  [Incidentally,  the  times  characterizing  significant  change 
in  the  theomodynamic  state  of  the  lower  tropical  ambient  are  probably  on 
the  order  of  100  hours,  a  time  span  partly  related  to  the  typical  magnitude 
of  eddy  viscosity  appropriate  for  the  tropical  atmosphere.  Hence,  unless 
the  intensification  time  from  depression  to  cyclone  greatly  exceeds  100 
hours  --  and  this  seems  dubious  —  the  thermodynamic  state  of  the  ambient 
may  possibly  be  justifiably  taken  as  fixed  throughout  the  later  stages  of 
intensification.  This  statement  by  no  means  implies  that  the  ambient 
can  be  taken  as  fixed  throughout  the  entire  intensification  process.  The 
very  fact  that  descent  speeds  in  the  outer  portions  of  even  a  mature 
hurricane  are  roughly  comparable  to  ambient  tropical  descent  speeds  suggests 
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that  the  tropical  ambient  probably  cannot  be  held  fixed  throughout  intensi¬ 
fication.] 

The  air  erupting  from  the  boundary  layer  begins  to  displace  the 
air  initially  in  the  cenxral  core  of  the  developing  storm.  The  air 
initially  present  in  the  ore  is  of  slightly  lower  pressure  than  the 
ambient  air  at  the  edge  of  the  storm,  but  not  vastly  different  in  vertical 
stratification.  The  air  erupting  from  the  boundary  layer  under  the  Rankine- 
vortex-1 ike  swirl  displaces  the  air  initially  present  in  the  core  vertically 
upward;  since  there  is  a  "lid"  on  top  of  the  storm,  the  vertically  dis¬ 
placed  initial  air  is,  near  the  top  of  the  core,  squeezed  radially  outward. 

The  following  competition  develops.  The  new  air  rising  out  of  the 
boundary  layer  is  drawn  entirely  from  relatively  warm  moist  air  near  the 
bottom  of  the  atmosphere.  Thus,  displacing  the  air  initially  present  in 
the  core  is  relatively  light  air.  On  the  other  hand,  the  convective  motion 
of  new  air  is  small,  especially  at  early  times,  and  the  ambient  processes 
(turbulent  diffusion,  radiational  cooling,  cumulus  convection)  try  to 
maintain  the  original,  near-ambient  stratification  in  the  core.  If  the 
convective  displacement  wins  out,  then  the  core  becomes  lighter  and 
lighter,  relative  to  a  column  of  air  at  the  outer  edge  of  the  storm.  This 
paragraph  contains  those  points  of  Carrier's  intensification  model  that 
seem  most  unsettling,  and  will  be  subject  to  close  scrutiny  in  the  critique 
below. 

However,  Carrier’s  basic  proposition  is  that  the  swirling  in  I  has 
led  to  a  downflux  in  II,  a  spiraling  inward  in  the  boundary  layer  and  an 
upflux  into  the  core,  and  a  lightening  of  the  core  by  hydrostatic  consider¬ 
ations.  For  dynamic  consistency,  the  centrifugal  force  (anticipated  to 
be  the  dominant  inertial  effect)  must  increase  to  balance  the  augmented 


radial  pressure  gradient.  Since  angular  momentum  is  conserved  in  I,  where 
friction  is  negligible,  the  fluid  particles  must  necessarily  move  in  closer 
to  the  axis  of  symmetry  (axis  of  rotation).  The  result  is  that  in  time, 
in  the  Rankine-vortex-like  swirl  distribution,  the  maximum  azimuthal  speed 
increases  in  magnitude  and  the  position  of  the  maximum  lies  closer  to  the 
axis  (Figs.  19  and  20).  Hence,  the  more  the  pressure  falls  in  the  core, 
the  more  fluid  sinks  into  the  boundary  layer  to  spiral  inward,  erupt 
upward,  and  cause  further  pressure  reduction  in  the  core.  If  the  crucial 
early  competition  is  resolved  in  favor  of  the  organized  convection,  ulti¬ 
mately  the  particles  erupting  out  of  the  boundary  layer  rise  so  quickly 
that  they  lie  on  a  moist  adiabat,  and  the  greatly  lightened  core  is  entirely 
flushed  of  its  original  fluid. 

While  most  descriptions  (Palmen  and  Newton,  1969)  tend  to  picture 
the  eye  as  formed  gradually  as  the  pressure  deficit  develops,  in  the  Carrier 
model  the  central  core  may  well  be  completely  flushed  of  ambient-like  air, 
so  that  the  air  in  the  core  lies  on  a  moist  adiabat,  before  much  trace  of 
an  eye  is  to  be  found.  At  first,  a  Rankine-vortex-like  swirl  holds  every¬ 
where.  From  this  fully  developed  one-cell  structure  with  at  most  only 
the  rudiments  of  an  eye,  a  well-defined  two-cell  structure  with  a  calm 
center  region  emerges  rapidly,  probably  in  much  less  than  an  hour,  owing 
to  inertial  oscillation,  in  the  following  way. 

As  the  pressure  falls  in  the  core  relative  to  ambient,  the  particles 
in  I  necessarily  move  in  closer  to  the  axis  to  permit  a  compensating 
centrifugal  force  to  develop.  Once  the  core  is  flushed  and  moist  adiabatic 
ascent  characterizes  the  full  height  of  the  core,  no  further  pressure 
deficit  can  be  generated.  By  inertia,  the  spinning  particles  continue  to 
move  in,  a  dynamic  imbalance  is  created,  and  a  radial  acceleration  develops 
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Fig.  19.  In  this  schematic  view  by  Carrier  of  the  flow  configuration  and 
circumferential  velocity  distribution  in  an  intensifying  tropical 
depression  at  some  early  time  t  «  t-|  (say),  the  interface  C-C 
between  the  new  and  initial  air  in  the  core  is  idealized  as 
horizontal  for  convenience.  However,  it  is  easy  to  show  that  for 
all  but  very  modest  circulations  (i.e.,  excent  for  angular 
velocities  less  than  three  times  ambient),  a  rigidly  rotating 
core  (which  implies  uniform  updraft  velocity  at  the  ton  of  the 
nonlinear  frictional  layer)  could  not  accept  all  the  fluid  pumped 
through  the  Ekman  boundary  layer.  Thus  a  uniform  ascent  in  the 
core  is  dynamically  impossible,  and  in  face  ascent  in  the  eye- 
wall  is  concentrated  in  cumulonimbi.  For  the  gross  balances  being 
discussed,  many  points  can  be  made  without  accounting  for  such 
refinements  in  eyewall  structure,  though  ultimately  details  of 
the  structure  are  crucial.  (From  Carrier,  1971b,  n.  146). 
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Intensification  from  tropical  depression  to  hurricane  has  pro¬ 
gressed  to  a  more  advanced  stage  in  this  schematic  diagram, 
holding  at  t  =  t2  >  tls  according  to  Carrier's  model.  The  maqni 
tude  cf  the  maximum  swirl  is  increased  and  its  position  lies  ' 
closer  to  the  axis  of  symmetry.  (From  Carrier,  1971b,  p  146  ) 


to  force  the  particles  away  from  the  axis  of  symmetry.  This  reverse  motion 
creates  a  rarefaction  at  the  center,  and  relatively  dry  warm  motionless  air 
sinks  down  the  axial  column  to  form  an  eye.  This  air  may  be  air  from 
above  the  storm  or  rained-out,  slowly  swirling  air  entrained  out  of  the 
top  of  the  moist  adiabatic  column  (Fig.  21).  Because  there  is  no  appreciable 
swir'  (hence  no  associated  pressure  gradient)  in  the  eye,  there  is  no 
frictional  boundary  layer  under  the  eye.  The  moist  adiabatic  column 
becomes  an  annulus  displaced  from  the  axis,  i.e.,  the  eyewall;  the  inertial 
oscillations  of  the  eyewall  eventually  damp  in  time. 

4 .2  Critique  of  the  Carrier  Model  of  Intensification 

Carrier  has  offered  more  than  the  conventional  linearized  stability 
analysis  to  determine  what  scales  of  motion  grow  in  what  time  interval 
(Yanai,  1964;  Ogura,  1964;  Yamasaki,  1969).  [Bates  (1970)  did  study  the 
growth  to  finite  amplitude  in  the  ITCZ  of  nonzonally  symmetric  eouatorial 
wave  disturbances  of  2000-km  longitudinal  scale,  which  are  observed  to  move 
westward  at  13  kts.  The  time  scale  for  early  growth  owing  to  barotronic 
instability  was  about  two  days;  frictional  dissipation  eventually  curtails 
further  growth  effected  by  conversion  of  condensational  heat  to  eddy  avail¬ 
able  potential  energy.  Because  Charney's  concept  that  condensational  heating 
in  the  tropics  is  dependent  on  convergence  of  moisture  in  the  planetary 
boundary  layer  was  incorporated,  the  disturbance  had  a  warm  core  with 
maximum  intensity  in  the  lower  atmosphere.  Holton  (1972)  Furnishes  further 
lucid  insight  into  equatorial  wave  disturbances,  and  in  fact  concisely 
discusses  many  other  topics  arising  in  tropical  meteorology.]  Carrier  has 
also  emphasized  the  role  of  lightening  of  the  core  rather  than  (say)  large- 
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Schematic  picture  by  Carrier  of  the  flow  configuration  which 
prevails  when  the  radius  of  maximum  swirl  R  is  increasing  and 
an  incompletely  formed  eye  is  being  filled  with  relatively  drv 
and  motionless  air,  which  sinks  down  from  the  top  of  the  storrr 
under  dry-adiabatic  compression.  With  the  formation  of  an  eve 
lighter  in  weight  than  the  eyewall,  the  terminal  stages  of 
intensification  and  the  beginnings  of  ouasisteadv  mature-staqp 
structure  are  realized.  (From  Carrier,  1971b,  p.  150.) 


scale  horizontal  turbulent  diffusion.  However,  there  are  problems  con¬ 
cerning  the  nature  of  the  convective  lifting  in  the  core,  as  now  discussed. 

The  ambient  distribution  of  total  stagnation  enthalpy  is  such 
that  when  a  column  stratified  like  the  ambient  is  lifted  as  a  whole 
adiabatica ily,  warm  air  is  squeezed  radially  outward  at  the  top,  even  if 
warm  air  is  added  at  the  bottom  (Fig.  7).  The  net  effect  for  small  con¬ 
vective  motion,  according  to  a  method-of-characteristics  numerical  integra¬ 
tion  of  the  resulting  hyperbolic  boundary-value  problem,  is  not  a  lightening 
of  the  developing  core  (Dergarabedian  and  Fendell,  1972a,  pp.  62-67). 13 

A  closely  related  observation  is  to  note  Carrier  (1971b),  in  his 
discussion  of  intensification,  seemed  more  concerned  with  how  a  hurricane 
is  formed  than  whether  it  is  formed.  Hence,  he  adopted  effectively  the 
following  proportionality: 

_ _ pressure  fall  from  ambient _ 

total  pressure  fall  for  a  moist  adiabatic  core 

_  _ new  air  into  core _ 

total  air  required  to  flush  core 

This  statement  is  true  at  the  outset,  and,  if  a  hurricane  forms,  at  the 
end,  of  flushing  the  core.  In  fact,  however,  such  a  simplification  lets 
a  hurricane  always  be  formed.  In  truth,  the  pressure  fall  for  small  times 
is  less  than  that  predicted  by  the  fractional-volume  proportionality.  As 
the  quasisteady  mature  state  approaches,  the  initially  slow  pressure  fall 
probably  accelerates  to  a  rapid  decrease,  and  the  thermodynamic  stratifi¬ 
cation  of  the  emerging  eyewall  asymptotically  approaches  a  near-moist- 
adiabat  state. 

The  problem  cannot  be  satisfactorily  buried  in  parameterizations 
of  turbulent  transfer,  radiative  cooling,  and  cumulus  convection  in  the 
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core.  What  appears  needed  is  explicit  recognition  that  a  broad  uniform 
ascent  in  the  core  does  not  account  for  the  actual  mechanism  of  convection, 
which  would  seem  to  entail  ascent  in  a  concentrated  number  of  intensified 
cumulonimbi.  The  somewhat  tentative  nature  of  this  remark  stems  from 
statements  like  that  by  Palmen  and  Newton  (1969,  p.  490):  "It  seems 
clear  from.  .  .  observations  and  from  descriptions  of  turbulent  encounters 
that,  in  some  cyclones,  ascent  in  the  inner  region  is  dominated  by  cumulo¬ 
nimbus,  .  .  .  while  in  others  a  broad-scale  and  more  uniform  ascent  takes 
place." 

If  cumulonimbi  do  play  an  essential  role  in  the  core  ascent  during 

intensification,  some  semblance  of  an  eye  would  probably  form  at  an  earlier 

1 4 

stage  of  hurricane  development  than  Carrier  suggests.  The  eye  would  probably 
not  be  entirely  created  on  a  time  scale  of  20  min  at  the  end  of  intensifi- 
'  cation,  as  Carrier  suggests,  although  inertial  oscillations  at  the  end  of 

intensification  might  rapidly  delineate  a  previously  ill-defined  eye. 

The  tritium-tracing  analyses  by  Ostlund  (1968)  of  Hurricane  Hilda 
(1964)  and  Hurricane  Betsy  (1965)  suggest  that  while  the  eye  may  have  par¬ 
tially  consisted  of  air  entrained  from  the  eyewall,  there  had  been  substantial 
subsidence  of  stratospheric  air  down  into  the  eye.  Thus,  Carrier  seems 
correct  in  citing  two  possib  e  sources  for  the  air  in  the  eye. 

4.3  The  Distribution  of  Cumulonimbi  during  Intensification 

Riehl  (1972,  p.  249)  notes:  "Very  fine  weather,  with  hardly  a 
cloud  in  the  sky,  often  precedes  the  arrival  of  a  hurricane  by  1  day,  when 
the  storm  is  moving  westward  in  the  tropics.  In  the  past  this  deceptively 

4  beautiful  weather  has  been  the  cause  of  many  disasters  with  heavy  loss  of 
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human  life  in  the  Caribbean  and  elsewhere."  In  a  possibly  related 
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observation,  Lopez  (1972b)  calls  attention  to  the  work  of  Oliver  and 
Anderson  (1969),  in  which  satellite  photographs  are  used  to  indicate  that 
the  majority  of  hurricanes  form  around  a  circulation  center  in  the  clear 
area  ahead  of  (rather  than  directly  under)  a  convective  region.  The 
concentration  of  cloudy  convective  regions  amid  clear  calm  regions  in 
tropical  cyclogenesis  may  be  related  to  characteristics  of  subsidence 
that,  by  continuity,  must  accompany  convection. 

However,  perhaps  also  there  is  a  concentrating  of  pre-existing 
dispersed  cumulonimbi  that  is  a  concomitant  of  large-scale  tropical  cyclo¬ 
genesis.  This  admittedly  speculative  subject  is  included  because  it  may 
be  worth  pursuing.  Furthermore,  Malkus,  Ronne,  and  Chafee  (1961)  observed 
that  for  the  inner  rain  area  of  radius  200  n  mi  in  Hurricane  Daisy  (1958), 
in  the  pre-hurricane  stage  there  were  60  cumulonimbi  that  covered  one 
percent  of  the  area;  during  intensification  two  and  one-half  percent  of 
this  area  was  covered  by  cumulonimbi.  In  the  mature  stage  there  was  a 
persistence  of  cumulonimbus  patterns  relative  to  the  center,  and  there 
were  200  cumulonimbi  that  covered  four  percent  of  the  area.  By  studying 
measurements  in  Hurricanes  Carrie  (1957)  and  Cleo  (1957),  as  well  as 
Daisy,  Gentry  (1964)  estimated  five-percent  coverage  of  the  inner  pighty- 
nautical-mi le-radius  area  by  cumulonimbi  (total  cloud  coverage  in  this  area 
was  under  25%).  Gentry  adds  that  both  aircraft  penetration  and  radar 
returns  suggest  that  the  area  of  strong  convective  activity  is  larger  in 
more  intense  hurricanes  and  is  located  closer  to  the  center. 

It  might  also  be  speculated  [in  contrast  to  conjectures  by  Gray 
(1972,  p.  64)  on  conditions  for  tropical  cyclogenesis]  that  for  a  synoptic- 
scale  pressure  difference  to  develop  not  all  the  air  ascending  in  the 
cumulonimbi  should  locally  subside,  but  some  of  the  descent  should  be 


119 


deferred  to  a  larger  horizontal  scale. 


4 . 4  The  Time-Dependent  Flowfield  during  Intensification 

Although  the  energetics  of  the  core  may  need  revision,  Carrier's 
modeling  of  the  dynamics  in  the  outer  portions  during  the  transient  phase 
should  remain  substantially  unchanged.  Thus,  this  portion  is  developed 
here  in  some  detai 1 . 

For  axisymmetric  intensification,  a  set  of  uniformly  valid  equations 
for  the  dynamics  is  now  given  in  cylindrical  polar  coordinates.  If  sub¬ 
scripts  r,  z,  and  t  denote  partial  differentiation,  for  a  noninertial 
coordinate  system  rotating  at  the  speed  of  the  normal  component  of  the 
rotation  of  the  earth  (i.e.,  normal  to  a  plane  locally  tangent  to  the 
earth  at  about  15°  latitude  so  ri  =  0.06  rad  hr  ^ ) : 

2  - 1 

ut  +  u  ur  +  w  uz  -  ?MV  -  (v  /r)  +  P  Pr  *  V  uzz, 

(rv)t  +  u(rv)r  +  w(rv)z  +  2ftru  =  v  (rv)zz, 

Pz  +  P9  =  0, 
and 

{ru)^  +  rwz  =  0.  (47) 

These  equations  hold  i n  0  <  r  <  r  ,  0  <_  z  <_  z-^  where,  again,  z-j  is  the 

top  of  the  storm  (a  slippery  lid)  and  rQ  is  the  outer  edge.  The  boundary- 

initial  conditions  of  relevance  here  are 


rQ:  u  =  v  =  0;  z=0:  u  =  v  =  w  =  0; 

(48) 

t  :  u,v  small . 

(49) 
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(44) 

(45) 

(46) 


There  is  angular  momentum  in  the  initial  flowfield  in  an  inertial  frame 
of  reference,  and  this  is  conserved  with  the  anticipated  radial  inflow 
toward  the  axis  of  symmetry,  the  picture  of  the  flow  is  that  given  in 
Fig.  19. 

For  convenience  the  following  notation  is  adopted.  The  radial, 
azimuthal,  and  axial  velocity  components  will  be  denoted  by  capitals 
(U,  V,  W)  in  I  and  by  lower-case  letters,  u,  v,  w  in  the  frictional  boun¬ 
dary  layer  II. 

From  results  given  earlier  for  the  mature  stage  boundary-layer 

dynamics,  and  from  other  considerations  developed  below  for  the  transient 

boundary  layer,  it  is  anticipated  that  f  =  r  V  is  a  linear  function  of 
2 

r  ,  and  independent  of  z;  then  there  is  a  downflux  into  the  frictional 
boundary  layer  II  from  the  supply  region  I  which  is  independent  of  r  to 
lowest  order.  It  follows  from  (45)  and  (47),  together  with  (48)  and  (49): 


and 


Y  =  r  V(r,t) 


<J>  =  r  U(r,t)  =  R  R 


9 


(50) 


(51) 


W(r,z,t) 


2  R  R 


z  , 


(52) 


where  R( tQ )  =  RQ.  The  associated  pressure  field  from  (44)  and  (46)  is 
discussed  later.  This  solution  holds  in  I,  i.e.,  in  R  i  r  <_  r  .  At 
t  =  tQ,  ini,  i.e.,  in  RQ  <.  r  <_  rQ,  f  =  0.  The  particles  at  r  =  Rq  at 
t  =  t  .  at  a  later  time  t  >  t,  lie  at  r  =  R  <  R  .  Typically,  r  = 

0  u  0  0 
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0(1000  mi),  R0  =  0(200  mi),  and  (RQ/12)  <  R(t)  <  RQ. 

The  next  step  is  to  develop  the  boundary  layer  in  II  that  lies  under 
the  swirling  flow  in  I  given  by  (50)  -  (52).  The  equations  for  the 
viscous  layer  are  taken  to  be  the  linear  transient  parabolic  eauations 
appropriate  for  low-Rossby-number  flow.  While  the  inequality 
($2  +  y2)  «  a2r4  may  not  formally  hold  everywhere  in  the  flowfield,  it  is 
widely  valid  and,  in  addition,  from  experience  with  the  mature  quasisteady 
storm,  the  key  output,  the  boundary-layer  divergence,  is  adeauately  given 
by  linearized  theory.  Hence, 


v  $Z7  ■  =  ~2o.ii, 

(53) 

v  i^zz  ■  =  2n<|>, 

(54) 

where 

$  =  r[u(r,z,t)  -  U(r,t)],  *  =  r[v(r,z,t)  -  V(r,t)], 
w  =  W(r)  +  w‘(r,z,t)  . 


(55) 


If  x  =  4>  + 


V  Xzz  -  Xf  =  2nix» 


$r  +  4>r  +  r  w'z  =  0  , 


(56) 

(57) 


subject  to 

x(r,z=0,t)  -*•  - ( <t>+i ) ,  xlr.z-^.t)  -*■  0,  x(r*z»0)  =  0  •  (58) 


if  (r  2  _  r2)  -+  r  an  excellent  approximation,  then  it  is  convenient  to 
'  o  ’  o 

define 
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(59) 


F(r.t)  =  F-| ( t )  (l  ~  ^4)=  [R02  -  R2(t)]  (l  -  ^4). 

ro  '  ro 

2 

Carrier  (1971a)  studied  the  special  form  F](t)=  PQ  exn(at).  For  the 
more  general  form  (59),  by  Laplace  transformation  (Dergarabedian  and  Fendell , 
1972,  pp.  56-57),  the  relevant  boundary-layer  divergence  is 


w(r  ,z-**,t) 


J'2 

~F 


Re 


1 

5rrT 


C  +  ia 

/ 

C-loo 


Ft  (S) 


2a  i 


(s  +  2ai) 


Yj2  CexP (st ) ]  ds 


‘  (a)  77(7 It  ■  a)  (Ro2  -  r2) 


(60) 


for  the  R(t)  of  practical  interest.  The  divergence,  a  function  of  time 
only,  is  largest  in  magnitude  in  the  steady  state;  if  R  -*■  (RQ/12)  typically. 


w(r,z^»)  =  -  &)  (vn)1/2  =  0(10"3 


mph) 


The  net  volumetric  flux  down  into  II,  and  hence  (in  an  incompressible  model) 
up  into  the  core,  Q(t)  is  given  by 


« ■  Qo  (rtk  ■  ')  [’  -  (it-)2] 


(61) 


where  the  quasi  steady  discharge  into  the  base  of  the  core  III  is  given  by 

%  =  ^ro2  "  r2)  (Ro/ro)2  (vR)1/2  • 

The  previously  discussed  requirement  for  compatibility  between  (1) 
the  nonlinear  spin-up  of  the  inviscid  flow  in  I  under  conservation  of 
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angular  momentum  and  (2)  the  pressure  deficit  generated  between  a  column  in 
the  core  and  a  column  in  the  ambient  owing  to  partial  flushing  of  the  core 
with  boundary  layer  air  driven  radially  inward  by  the  outer  inviscid  flow, 
is  now  made  quantitative.  It  is  necessary  to  model  the  radial  and  azimuthal 
velocity  components  in  the  core;  kinematic  and  continuity  considerations 
suggest  that,  at  least  in  the  lower  regions  of  III,  for  this  purpose  it  is 
adequate  to  take 


r  u  =  <t>  =  r  R/R  , 


r  v  »  4,  =  n(RQ2  -  R2)  r2/R2 


(63) 

(64) 


Accordingly  (44)  is  integrated  from  r  =  0  to  r  =  rQ  at  some  fixed  z  above 
the  boundary  layer;  the  profiles  for  u  and  v  are  given  by  (63)  and  (64)  for 
0  <  r  <  R  and  by  (50)  and  (51)  for  R  <  r  <  rQ.  The  term  u  ur  in  (44)  is  a 
perfect  differential  and  yields  precisely  zero;  w  uz  yields  zero  because  u 
is  independent  of  z  and  w  is  small  anyway.  It  follows  that  the  time- 
dependent,  Coriolis,  centrifugal,  and  pressure  gradient  terms  yield,  for 
p  =  const,  and  neglecting  terms  of  0(R2/rQ2)  relative  to  those  of  order 
unity. 


p  a2  R02  <R02 

p(r  ,  z=G,t)  -  p(r=0,  z=0,t)  - 2 — 

U  D1- 


+ 


p  ^  [r  R  In 


(65) 


The  last  term  on  the  right-hand  side  is  characterized  by  the  inverse  square 
of  the  time  for  intensification;  this  time  probably  is  much  larger  than  the 
time  for  equilibration  of  the  boundary  layer.  Thus  normally  the  time- 


124 


derivative  term  can  be  neglected. 

However,  when  inertial  oscillations  associated  with  the  continued 
radial  influx  and  rebound  of  air  in  I  occur,  a  shorter  time-scale  enters 
and  the  time-derivative  term  must  be  retained.  Carrier  (1971b,  p.  157) 
suggests  that  the  time-scale  of  the  oscillation  is  0(4R^  /ftRQ  )  where  R^ 
is  the  final  equilibrium  position  of  the  eyewall  R(t).  For  typical  values, 
with  RQ  =  250  mi,  this  time  scale  is  20  min.  In  fact,  holding  the  convection- 
core  pressure  deficit  from  ambient  fixed  at  about  35  mb,  holding  p  at  its 
ambient  sea-level  value,  letting  RQ  =  250  mi  and  rQ  =  1000  mi,  one  finds 
from  (70)  --  by  dropping  the  temporal  derivative  --  that  R  -  Rf  =  0(35  mi). 
Restoring  the  temporal  term,  adopting  as  initial  conditions  R  =  0(35  mi) 
and  R  =  -n(RQ2  -  R2)/2R  (Carrier,  1971b,  pp.  152  and  155),  and  holding  the 
core  pressure  deficit  relative  to  ambient  fixed  at  35  mb  (near-quasi steady 
conditions),  one  finds  from  (70)  that  R  decreases  to  about  16  mi  in  about 
a  half-hour  before  R  =  0.  Thus,  numerical  values  confirm  the  inertial 
oscillation  period  suggested  by  Carrier.  Such  oscillations  probably  reinforce 
the  eye,  rather  than  create  it,  as  noted  earlier  (Section  4.2).  The  oscil¬ 
lation  of  the  eyewall  is  dwelled  upon  because  it  seems  to  be  the  first 
proposal  of  a  mechanism  i nternal  to  the  hurricane  that  might  explain  a 
crudely  periodic  variability  in  intensity  of  10  to  20%  in  peak  winds,  on 
the  time-scale  of  a  few  hours;  such  variability  appears  to  be  observed 
sometimes,  and  attribution  to  changes  in  the  ambient  atmospheric  and  under¬ 
lying  oceanographic  states  may  not  be  always  satisfactory. 

The  increased  low-level  radial  inflow  toward  the  center  of  the 
developing  cyclone  seems  consistent  with  the  previously  discussed  reorgani¬ 
zation  of  the  statistically  steady  cumulonimbus  activity  which  helps 
maintain  the  ambient  enthalpy  transport  process.  Such  radial  inflow  implies 
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a  stronger  circulation  (or  vorticity),  and  this  implies  a  decreased  mixing 
of  updraft  air  with  surrounding  air;  the  result  would  be  a  decreased 
density  in  the  updraft  region  and  hence  a  decreased  pressure.  As  cumulo¬ 
nimbus  activity  increases  in  the  central  region  of  the  incipient  cyclone, 
there  is  probably  diminution  of  such  activity  iri  the  surrounding  region. 
However,  the  downdraft  in  the  surrounding  region  should  differ  little  from 
ambient  since  the  total  low-level  inflow  into  the  whole  region  will  not 
have  changed  appreciably.  Thus,  the  anticipation  is  that  the  flowfield 
presenttxi  in  this  section  probably  remains  substantially  correct,  within  the 
axisymmetric  approximation. 
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5.  CONCLUDING  REMARKS 


Several  viewpoints  not  universally  accented  have  been  adopted  in  this 
review.  The  existence  of,  and  necessity  for,  a  greatly  augmented  oceanic  heat 
source  in  the  inner  regions  of  a  hurricane  as  the  basic  sustaining  mechanism 
has  been  questioned.  The  validity  of  the  hot-tower  thesis  of  heating  from 
the  tropical  cumulonimbus  has  been  scrutinized.  The  incomplete  under¬ 
standing  of  the  processes  and  time  scales  of  tropical  cyclogenesis  and 
intensification  has  been  emphasized.  The  conviction  that  improved  quanti¬ 
tative  understanding  of  the  tropical  ambient  is  prerequisite  to  further 
progress  on  the  deepening  of  tropical  disturbances  has  been  stated.  Reasons 
for  reservations  concerning  the  probable  success  of  current  hurricane 
seeding  techniques  have  been  off erred. 

This  unconventionality  hopefully  conveys  two  major  points  to  the 
reader.  First,  any  serious  worker  in  hurricanes  would  be  very  remiss  not 
to  consult  the  referenced  works  directly  and  decide  for  himself  whether 
or  not  the  viewpoints  adopted  here  are  in  fact  valid.  Second,  this  review 
on  tropical  cyclones  describes  a  vital  subject  in  rapid  flux  to  which  many 
basic  physical  and  mathematical  contributions  are  still  to  be  made.  This 
is  not  a  field  in  which  only  exquisite  refinements  to  elaborate  computer 
programs  remain.  In  fact,  many  valuable  theoretical  contributions  in  recent 
years  have  been  achieved  by  Carrier  and  Gray  via  approximate  solution  of 
simplified  models  with  only  modest  computing,  as  opposed  to  attempting 
direct  numerical  simulation  of  the  full  conservation  laws.  The  meteorolo¬ 
gist  with  insight  will  hopefully  agree  that  the  tropical  cyclone  is  neither 
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a  closed  subject,  nor  a  subject  to  which  only  those  with  access  to  very 
sophisticated  computing  facilities  can  aspire  to  contribute. 
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APPENDIX  A 


ESTIMATING  THE  KINETIC  ENERGY  AND  WATER  CONTENT  OF  HURRICANES 


To  estimate  very  quickly  and  roughly  the  total  kinetic  energy  con¬ 
tent  of  a  mature  hurricane  of  moderate  intensity,  one  must  evaluate  I  where 

*  2=2.j/2  r=rQ 

I  =  ^  /  p  v2  dV  =  71  f  dz  f  p[v(r)]2  r  dr  .  (A.l) 

V  *z=0  r=r 

e 

The  speed  is  virtually  all  swirl,  and  about  half  the  nine-mile-high, 

axisymmetric  storm  is  rotating  like  (say)  a  potential  vortex.  The  swirl 

is  (suppose)  100  mph  at  the  edge  of  the  (nonrotating)  eye  rg  =  30  mi,  and 

the  storm  extends  out  to  rQ  =  500  mi.  The  density  may  be  set  to  a  con- 

-4  -3 

stant  average  value  of  5  x  10  g  cm  .  Then 


I  = 


TTPZ1 

— 


(v  )‘ 
'max' 


(A. 2) 


Battan  (1961,  p.  21)  characterizes  the  kinetic  energy  of  a  hurricane  as 
in 

10  kilowatt-hours,  or  3.6  x  10  ergs. 

An  upper  bound  on  the  rainfall  per  day  for  a  hurricane  hovering  over 
a  spot  is  now  calculated  according  to  the  Carrier  model.  The  vertical 
velocity  down  into  the  boundary  layer  in  the  mature  stage  is 

R  2 

w(r ,  z  -v  <»,  t  -*.»)  =  -  —^2  (vn)1^2  .  (A. 3) 
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The  volumetric  flux  down  into  the  boundary  layer  is,  for  rQ  >>  R  , 

Qq  =  7i  rQ2  w(r,  z->”,  t  -*■“)=-  n  Rq2  (vfi)1^2  ;  (A. 4) 

this  is  also  the  flux  up  into  the  eyewall.  The  mass  of  air  per  unit  time 
is  given  by  p  Qo  where  p  is  fairly  constant  at  a  fixed  height.  The  water 
content  of  this  mass  of  air,  which  falls  out  as  precipitation,  is  given 
in  mass  per  time  by 

J  =  *p  R02  Y(z  =  zt)  ,  (A. 5) 

where  z is  a  typical  ambient-profile  height  for  characterizing  the  mass 

•  ”3 

fraction  of  water  vapor  for  the  flux  into  the  eyewall.  For  p  =  1.1  x  10 

g  cm"3,  v  =  2  x  105  cm2  sec"1,  a  =  6.25  x  10'2  hr'1,  and  [from  the  data  of 

^  Jordan  (1957)]  Y(zt  =  1000  ft)  -  1.65  x  10"2,  J  =  1.8  x  1011  g  sec'1.  If 

-3 

this  falls  entirely  in  an  annulus  from  20  to  30  miles,  one  gets  4.3  x  10 
g  cm" L  sec"  ,  or  since  1  g  of  liquid  water  occupies  about  1  cm  ,  about 
370  cm  day"1.  The  record  rainfall  over  a  spot  in  one  day  is  associated 
with  a  hurricane:  117  cm  (Baguio,  Phillippines  in  July  1911).  Of  course, 
for  a  hurricane  translating  at  20  mph,  no  point  is  likely  to  lie  under  the 
eyewall  for  more  than  an  hour,  and  (370/24)  =  15  cm  hr  . 

Because  of  the  flood-spawning  torrential  rains  over  Virginia  associ¬ 
ated  with  Hurricane  Camille,  it  is  interesting  tc  consider  whether  the  water- 
vapor  content  of  this  hurricane  at  landfall  is  consistent  with  the  large 
total  rainfall  over  the  Southeast.  From  isohyetal  maps  published  by 
Schwarz  (1970)  it  may  be  crudely  estimated  that  in  eastern  Kentucky, 

Virginia,  and  West  Virginia  an  area  of  150  mi  x  180  mi  received  an  average 
of  five  inches  of  rain  and  an  area  of  100  mi  x  360  mi  received  three  inches. 
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4  2 

In  addition  5  x  10  mi  of  Mississippi,  Alabama,  Louisiana  and  Tennessee, 

1 5 

received  an  average  of  three  inches.  This  totals  2.58  x  10  g,  which  is 
probably  a  conservative  estimate.  From  Jordan's  data  one  can  estimate 
that  the  Caribbean  hurricane  contained 


Q  = 


p  Y  dz  =  tt 


(A. 6) 


where  p  Y  =  py  =  density  of  water  vapor.  For  an  August-September  ambient, 

-  ? 

I,  =  4.71  g  cm  .  Equating  Q  to  the  total  rainfall,  one  finds  rQ  =  260  mi. 
This  is  an  interesting  value  in  that  Camille  was  a  small  intense  hurricane 
with  hurricane  winds  extending  outward  from  the  center  about  60  mi  and 
gale  winds,  110  mi  (Meyer,  1970,  p.  5).  However,  the  results  here  indicate 
merely  that  rQ  =  0(260  mi)  satisfies  a  conservative  estimate  of  the  total 
precipitation  requirement;  it  is  not.  intended  to  imply  rQ  could  necessarily 
be  taken  so  small,  though  if  one  fits  Carrier's  proposed  swirl  profile 
[given  below  (30)]  to  the  data,  rQ  =  0(225  mi). 
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APPENDIX  B 


THE  MOIST  ADIABAT 

Because  of  the  significant  role  it  plays  in  understanding  tropical 
cyclones,  a  derivation  of  the  relation  describing  the  moist  adiabatic  nrocess 
seems  worth  including,  The  following  is  based  on  unpublished  notes  of 
G.  F.  Carrier  of  Harvard  University. 

When  moist  air  expands  adiabatically,  the  temperature  eventually 
decreases  to  a  value  at  which  the  water  vapor  density  reaches  its  saturation 
level.  For  all  further  decreases  in  temperature,  condensation  occurs. 

Here  it  is  taken  that  the  condensed  water  has  zero  volume  (a  good  approxi¬ 
mation),  and  that  the  liquid  precipitates  out  as  it  is  formed. 

An  air-vaper  mixture  is  at  temperature  T,  the  mass  of  air  in  volume 
V  being  ma,  and  the  saturated  vapor  mass  in  V  being  mv-  The  corresponding 
densities  are  =  (ma/V)  and  py  =  (mv/V);  the  pressure  P  =  ba  +  Pv* 

A  relationship  between  T  and  p  (say)  as  T  changes  is  obtained  by  studying 
the  moist  adiabatic  process  in  several  steps  (Fig,  22).  The  mixture  is 
separated  into  two  fluids  (dry  air  and  vapor),  to  be  later  recombined; 
this  step  is  taken  to  require  no  expenditure  of  work. 

Initially,  the  mass  my  of  vapor  occupies  the  volume  V  as  in  the  top 
box  of  Fig.  22.  That  volume  is  decreased  by  a  process  in  which  dmy  of 
the  vapor  is  condensed  while  it  and  the  remaining  vapor  are  at  temperature 
T  and  pressure  py.  Thus,  the  new  volume  (that  of  the  second  box)  is 
V [1  -  (dmv/mv)].  The  work  done  by  external  forces  during  this  process  is 
(Pv  v  dmv/mv),  and  is  indicated  by  the  arrow  to  the  left  of  the  boxes. 
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The  heat  gained  L  dmv  is  shown  at  the  right.  The  vapor  now  expands  adia- 
batically  and  isentropical ly  to  volume  (V  +  dV)  as  shown  in  box  3;  the 

I 

new  temperature  T  is  given  by 


V[1  -  (dmv/mv)])v 

V  +  dV  j  =  1  "  'Yv 


(B.l) 


In  the  third  step,  the  gas  in  box  3  is  heated  to  temperature  (T  +  dT),  as 
indicated  in  box  4.  The  heat  that  must  be  added  during  this  process  is 


(B.2) 


where  (c  )  is  the  heat  capacity  at  constant  volume  of  the  vapor,  (c  )  the 
v  v  p  v 

heat  capacity  at  constant  pressure  of  the  vapor,  R  =  (c  )  -  (c  )  ,  and 

p  v  v  v 

\  “  V  /(cv>  • 
y  v  v 

The  dry  air,  on  the  other  hand,  starts  off  in  the  bottom  box  of 
volume  V  with  mass  m,  at  temperature  T.  It  expands  isentropically  to 

a 

volume  (V  +  dV)  so  that  its  temperature  becomes  T [1  -  (y=  -  1 ) ( dV/V )] . 

a 

Henceforth  y„  =  y4  =  y.  The  dry  air  is  then  heated  at  constant  volume  to 

V  a 

temperature  (T  +  dT)  by  the  addition  of  heat  in  the  amount  m  [(c  )  dT  + 

a  v  a 

RaT(dV/V)]. 

Thus  at  the  end  of  the  process,  a  vapor  mass  (my  -  dmy)  and  a  dry 
air  mass  ma  occupies  a  volume  (V  +  dV)  at  temperature  (T  +  dT).  By 
definition,  the  process  undergone  will  have  been  the  moist  adiabatic  one 
if  the  net  heat  added  is  zero,  and  if  the  final  (pv  +  dpy)  is  the  saturation 
value  for  the  final  temperature  (T  +  dT): 
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dV 


<c.)  <T  *  V  v 

3 


+  m. 


']  *  V 


(cv)  dl  +  R  T  ^|  +  R(T  dm,,  =  -L  dm.,  ; 
v 


(B.3) 


m 


»,  ■  f  =  »,(T) 


(B.4) 


where  p v ( T )  is  the  function  that  describes  saturation  density. 
By  the  equation  of  state  for  dry  air. 


dV  _  \  dT  dPa 

V  '  "  pa  '  T  ‘  Pa  • 


(B.5) 


But  (dpv/pv)  is  comprised  of  two  contributions:  that  due  to  change  of 
volume,  and  that  due  to  the  loss  of  mass  via  condensation.  Thus 


pv  +  dpv 


m.,  -  dm 


m. 


V  +  dV 


(B.6) 


Hence,  after  division  by  V,  upon  use  of  the  equations  of  state, 

M(cv>a  +  V  dT  -  dPa}  +  ♦  «v]  dT  -dP}=  -  7  d(pv  V)  ,  (R.7) 


(B.6) 


where  P(T)  is  the  saturation  vapor  pressure. 

Alternatively,  the  moist  adiabat  may  be  derived  as  a  spec’ll  case 
of  the  equation  for  conservation  of  energy,  written  as  the  sum  of  internal 
plus  kinetic  energy  for  a  multicomponent  mixture.  Here  Cartesian  tensor 
notation  will  be  used  instead  of  vector  notation,  and  a  comna  denotes 
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partial  differentiation.  Furthermore,  the  symbolism  previously  adopted 

in  the  main  text  for  describing  the  dynamics  of  an  effectively  one-component 

2 

gas  will  not  be  used  here.  The  air  velocity  is  denoted  u^ ,  with  u^u^  ~  q 
and  u •  •  o  =  w,  where  the  gravitational  acceleration  g.  =  -g  6..  r  The 

110  I  I  o 

water  vapor  velocity  is  (u.  +  v.).  The  vapor  quantities  p  ,  (c  )  ,  (c  )  , 

11  v  v  p  v 

R  and  the  air  quantities  p  .  (c)  ,  (c)  ,  Ra  have  the  same  designations 

V  a  V  p  a 

a  a 

as  above.  The  specific  internal  energy  of  the  air  is  ea ,  and  of  the  vapor 
e  .  The  specific  heat  of  the  liquid  water  substance  is  c„.  The  heat  lost 

3 

by  radiation  and  all  other  processes  is  0,  with  units  (mass/length-time  ). 

The  conservation  of  energy  for  current  purposes  may  be  written  (upon  neglect 
of  the  kinetic  energy  contribution  relative  to  internal  energy  contributions) 


[oa  ea  +  pv  ev]  .  +  [oa  ea  ui  +  °v  ev(ui  +  vi)]  . 

}t  >1 


=  -(pa  +  pv)  gw  (oij  Ui)  +  (kT,^  .  "  C!J0J  -  Q 

»j 


(B.9) 


where  a.,  is  the  total  stress  tensor,  k  is  the  thermal  conductivity,  and 

1  J 


pv  +  pv(ui  +  vi)  .  = 
,t  ,i 


(B. 10) 


In  general,  w  (the  mass  of  water  vapor  condensed  per  unit  time)  must  be 

specified  explicitly  or  through  a  prescribed  mechanism;  when  the  fluid  is 

saturated  and  (dT/dt)  <  0,  u  is  known. 

For  steady  flow  with  densities,  pressures,  and  temperature  dependent 

on  the  vertical  coordinate  z  only,  with  negligible  diffusion  so  v.  =0  and 

o..  =  -p  in  the  hydrostatic  approximation  p  _  +  pg  =  0,  it  follows 
ij  ij 

from  continuity  that 
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(B.n) 


(p  w)  =  0  ,  {p  w)  =  ; 

,z  ,z 

hence, 

Pa  w  ea  +  Py  w  ey  -  u  ey  +  p  WjZ  +  c;T^  =  0  .  (B.12) 

a  >  z 


With  the  aid  of  the  continuity  equations  and  the  eouations  of  state, 


w  Kvl  +  w  Kv 

-  n  J  _  L  »  «, 


■wp(Z=  „(ev  +  RyT  -  c.T)  .  (B. 13) 


Si  nee 


and 


(B.14) 


(B. 1 5) 


upon  substitution  and  cancellation  of  the  common  factor  w,  one  recovers 
(B.7) . 
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PARTIAL  LIST  OF  SYMBOLS 


English  Symbols 

heat  capacity  of  liquid  water 
c  heat  capacity  at  constant  pressure 

r 

cv  heat  capacity  at  constant  volume 
e  internal  energy 

E  Ekman  number,  v/y 

g  magnitude  of  gravitational  acceleration,  | g ^  | 
acceleration  of  gravity 
H  total  static  enthalpy,  CpT  +  gz  +  LY 

Ht  total  stagnation  enthalpy,  H  +  q  /2 

L  specific  latent  heat  of  phase  transition 
m..  molecular  weight  of  species  i 
n  power  of  algebraic  decay  of  swirl  with  radial  distance 

p  pressure 

P  saturation  vapor  pressure 

q  wind  speed,  | qi  | 

q.j  velocity  vector  for  a  fluid  in  noninertial  frame  rotating  with 
earth 

Q  net  volumetric  flux  per  time;  heat  loss  Der  volume-time 
Q0  steady-state  net  volumetric  flux  per  time 
r  radial  coordinate  in  cylindrical  polar  coordinates 
R  radial  position  of  maximum  azimuthal  velocity  v 
Rj  gas  constant  for  species  i,  W/m. 
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I 

* 

l 

R  universal  gas  constant 

RH  relative  humidity 

t  t  time 

\  T  temperature 

u  radial  velocity  component  of 

velocity  vector  for  dry  air 

U  radial  velocity  component  of  in  the  inviscid  flow  exterior 

to  the  frictional  boundary  layer 

1  v  azimuthal  velocity  component  of  q^ 

v..  velocity  vector  for  water  vapor  relative  to  dry  air;  velocity 
vector  for  a  fluid  in  a  noninertial  system 

V  azimuthal  velocity  component  of  q^  in  the  inviscid  flow  exterior 
to  the  frictional  boundary  layer;  volume 

fi 

9  w  vertical  velocity  component  of 

W  vertical  velocity  component  of  o^  in  the  inviscid  flow  exterior 

to  the  frictional  boundary  layer 
? 

x  fir  /i|<  ;  a  factor  in  the  moist  adiabat  based  on  dH  =  0 

Y  water  vapor  mass  fraction,  p v/p 
z  altitude  above  sea  level 


Greek  Symbols 


I 


*  cp/cv 

p  density 

4)  ru  [sometimes,  r(u  -  U);  sometimes,  nondimensional] 

♦  rU 

i|>  rv  [sometimes,  r(v  -  V);  sometimes,  nondimensional] 
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v  rV 


o 


V 

u 


a 


ratio  of  molecular  weight  of  water  vapor  to  dry  air 
total  stress  tensor 

kinematic  viscosity  (usually  given  turbulent  transfer  values) 
mass  of  water  vapor  condensed  per  volume-time 
Coriolis  parameter  (component  of  the  angular  velocity  of  the 
earth  perpendicular  to  a  plane  locally  tangential  to  the  sea 
surface  in  the  tropics) 

azimuthal  coordinate  in  cylindrical  polar  coordinates;  potential¬ 
like  enthalpy,  c  l  +  gz 

2 

stagnation  potential-like  enthalpy,  0  +  a  /2 


Superscript 

•  total  derivative  with  respect  to  time 

-*  vector  (Gibbs  notation) 


Subscripts 

a  dry  air;  ambient 

c  at  the  axis  of  summetry  at  sea  level 

e  eyewall  at  sea  level 

i  vector  (Cartesian  tensor  notation) 

o  outer  edge  of  hurricane;  initial  value;  typical  value 

s  ambient  at  sea  level 

t  total 

v  vapor 
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FOOTNOTES 


1  For  the  East  Coast  of  the  United  States,  the  hurricane-generated  storm 
surge  depends  on  storm  direction,  but  can  be  roughly  estimated  from 

s  =  2.90-9.70  x  10'2  (v)max  +  1 .33  x  10"3  (v)2ax,  74  <  (v)max  <  149,  where 

(v)  is  the  maximum  wind  speed  in  mph  and  s  is  the  surge  above  normal 

in  ft  [based  on  Savfir,  H.  S .  ( 1 97 3 ) .  Military  Engineer  65,  4].  In 

addition  to  onshore  wind-driven  waves,  with  the  arrival  of  a  hurricane 

and  persisting  often  for  about  a  day  after  its  passing  along  a  coast  with 

a  sloping  bottom  (say,  of  inclination  angle  3),  arise  so-called  edge 

waves.  These  waves  travel  parallel  to  the  coast;  their  crests  are  normal 

to  the  coast.  If  the  component  of  the  hurricane  translational  velocity 

parallel  to  the  coast  is  denoted  U,  the  period  of  the  waves  is  (2iTU/g  sin  B) 

o 

and  the  wavelength  is  ( 2ttU  /g  sin  s);  typically  for  the  East  Coast  of  the 
United  States  the  period  is  fi ve-to-seven  hours,  and  the  wavelength  is 
400  km.  The  amplitude  vanishes  rapidly  from  the  shore  seaward,  and  is 
negligible  at  the  distance  of  one  wavelength;  the  amplitude  is  roughly 
given  by  the  inverse  barometer  rule  (one  centimeter  water  rise  per  millibar 
of  atmospheric  pressure  drop)  so  heights  of  two-to-three  feet  are  typical 
[Munk,  W.,  Snodgrass,  F.,  and  Carrier,  G.  (1956).  Science  123,  127]. 

2 

According  to  data  collected  on  30  typhoons  for  1960-1970,  tropical  storms 
approaching  the  Philippines  with  maximum  winds  in  excess  of  90  kts  leave 
with  winds  reduced  by  40-50%,  while  a  storm  with  peak  winds  less  than  90  kts 
experience  only  a  10-15%  reduction  recovered  within  a  day  after  leaving 
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(and  storms  with  peak  winds  under  60  kts  typically  undergo  an  intensifica¬ 
tion  in  crossing  the  islands).  The  average  crossinq  time  is  14.6  hr,  with 
the  weaker  storms  crossing  more  quickly  [Brand,  S.,  and  Blellock,  J.  VI.  (1972). 
Changes  in  the  characteristics  of  typhoons  crossing  the  Philippines. 
Environmental  Prediction  Research  Facility  Technical  Paper  6-72,  Naval 
Postgraduate  School,  Monterey,  California]. 

3The  convective  instability  of  the  tropical  ambient  explains  how  cumulo- 
nimbi  can  occur,  but  cumulonimbi  explain  how  high-static-enthalpy  levels 
of  the  upper  tropical  troposphere  can  arise.  This  interdependence 
exemplifies  the  earlier  remark  on  the  link  between  the  presence  of  clouds 

and  the  circumstances  regarding  stability. 

^These  authors  emphasize  that  the  ITCZ  consists  of  westward  propagating 
cloud  clusters  that  pass  at  four-to-fiv e-day  intervals. 

^Further  discussion  of  hybrid  (or  semi-tropical)  storms,  which  rely  for 
energy  source  not  only  on  latent  heat  release  but  also  on  baroclinic it.y 
associated  with  positioning  of  warm  and  cold  air  masses,  is  given  by 
Spiegler  (1972). 

6Data  on  this  matter  appears  contradictory .  Using  all  recorded  data  for 
Atlantic  hurricanes  for  1961-1968,  Sheets  (1972)  found  no  correlation 
between  the  size  of  the  radar  eye  and  storm  intensity  indicators  (such  as 
maximum  wind  speed  or  minimum  sea  level  pressure);  there  was  reportedly 
a  slight  tendency  toward  smaller  eye  diameters  at  lower  latitudes. 

Questions  concerning  the  use  of  radar  by  itself  to  determine  the  eve  radius 
will  be  discussed  in  Section  3.4.  In  contrast,  the  eye  diameter  increased 
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typically  from  20  n  mi  to  33  n  ini,  while  the  intensity  typically  decreased 
by  33%,  according  to  data  for  typhoons  crossing  the  Philippines  during 
1960-1970;  furthermore,  for  typhoons  near  the  Philippines  the  eye  diameters 
for  intense  typhoons  (maximum  wind  speed  over  90  kts)  are  typically  13 
smaller  than  for  weaker  typhoons  (maximum  wind  speeds  under  90  kts). 
Incidentally,  if  the  size  is  determined  by  the  average  diameter  of  the 
closed  surface  isobar,  the  size  of  a  sto^m  decreases  by  17%  in  area  in 
crossing  the  Philippines,  and  near  the  Philippines  the  intense  typhoon  has 
a  mean  outer  circulation  diameter  60  n  mi  -  150  n  mi  greater  than  that 
of  a  weaker  typhoon  [Brand,  S.,  and  Blellock,  J.  VI.  (1972).  Changes  in  the 
characteristics  of  typhoons  crossing  the  Philippines.  Environmental 
Prediction  Research  Facility  Technical  Paper  6-72,  Naval  Postgraduate 
School,  Monterey,  California]. 

^ From  "Tropical  Meteorology"  by  Herbert  Riehl.  Copyright  1954  by  McGraw- 
Hill  Book  Company.  Used  with  permission. 

®Most  discussions  emphasize  the  cold  wake  left  in  the  unoer  oceanic  layers 
by  intense,  slowly  moving  hurricanes;  however,  the  reduction  of  surface- 
layer  temperature  by  upwelling  and  mixing  may  be  occurring  directly  under 
the  hurricane.  Black  and  Mallinger  (1972,  p.  74  )  find  evidence  of  radial 
outflow  and  downwelling  of  warm  upper-layer  water  for  Hurricane  Ginger 
(1971);  the  upwelling  of  colder  water  was  found  to  reduce  the  sea-surface 
temperature  by  4°C  on  one  day,  and  by  2.5°C  on  another  day  when  the 
hurricane  was  translating  faster.  Reduced  latent  and  sensible  enthalpy 
transfer  from  sea  to  atmosphere  ensued,  with  a  reduction  in  intensity. 
Thus,  the  interaction  of  a  hurricane  with  the  ocean  may  be  such  that 


144 


reduced  transfer  from  the  sea,  not  augmentation  of  transfer  over  ambient 
level,  may  occur  (cf.  Section  3.9.1). 

^ Kraus  emphasizes  that  +he  Riehl-Malkus  theory  depends  on  an  isothermal 
inflow  layer  and  discusses  whether  frictional  effects  are  fully  accounted 
for  in  that  theory:  "The  importance  of  [an]  additional  energy  supply  to 
air  in  the  hurricane,  before  it  ascends  in  the  warm  core,  was  stressed 
particularly  by  Malkus  and  Riehl  (1960).  The  argument  would  be  weakened 
somewhat  if  a  (verbal)  suggestion  by  Carrier  was  found  to  be  true.  Carrier 
maintains  that  surface  friction  must  cause  a  deceleration  and  compression 
of  a  surface  air  parcel  along  its  trajectory;  this  compression  compensates 
to  some  extent  for  the  expansion  that  would  be  associated  otherwise  with 
the  externally  imposed  pressure  reduction"  [Kraus,  E.  B.  (1972).  "Atmosrheric- 
Ocean  Interaction."  Clarendon  Press,  Oxford,  England,  p.  208]. 

10Black  and  Mallinger  (1972,  p.  64)  characterize  most  models  as  giving  a 
50%  reduction  in  maximum  wind  speed  for  a  2°C  drop  in  sea-surface 
temperature.  Doubt  concerning  such  sensitive  dependence  of  hurricane 
intensity  on  sea-surface  temperature  was  discussed  in  Section  3.8.  Recently 
the  existence  of  a  closed  oceanic  convection  cell  induced  by  hurricane 
winds,  already  discussed,  has  been  further  confirmed  by  additional  measure¬ 
ments  on  Typhoons  Trix  (which  was  situated  over  the  central  South  China 
Sea  on  October  23-24,  1952)  and  Wilma  (which  was  similarly  situated  on 
October  28-29,  1952),  and  on  other  storms  [Ramage,  C.  S.  (1972).  Weather 
27,  484].  Ramage  (1972,  p.  491)  notes:  "When  they  struck  the  Philippines 
Trix  and  Wilma  were  greatly,  equally  intense  typhoons.  Trix  reintensified 
over  the  South  China  Sea  and  Wilma  probably  did.  Paired  aircraft  recon¬ 
naissances  on  the  23rd  and  28th  and  the  24th  and  29th.  .  .  show  little 
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intensity  difference,  in  contrast  to  Ooyama's  numerical  hurricane  model 
•  •  •  in  which  a  smaller  surface  temperature  difference  is  associated 
with  a  significantly  large  intensity  difference."  Furthermore,  the  nath 
for  Wilma  was  much  less  influenced  by  oceanic  coolinq  owinq  to  the  prior 
passage  of  Trix  than  the  work  of  Brand  (1971)  (previously  cited  in  Section 
1.5)  would  suggest,  even  though  Brand  considered  typhoons  generally  more 
separated  in  time.  This  is  one  more  item  of  evidence  suqgestinq  that  ocean 
temperatures  are  not  solely  determining  concerning  where  hurricanes  occur. 
Ramage  (1972,  p.  493)  concludes:  "The  history  of  typhoon  Wilma  sugqests 
that  surface  temperature  may  not  always  control  storm  intensity  nor 
significantly  affect  storm  movement,  for  Hi lma  moved  along  the  cold  water 
wake  of  Trix  and  remained  intense  until  reaching  Indochina." 

]  i 

'Typically  five  runs  of  the  kind  described,  at  two-hour  intervals,  starting 
from  the  radius  of  maximum  winds  and  proceeding  radially  away  from  the 
center,  characterize  the  so-called  eyemod  experiment.  Most  of  the  discus¬ 
sion  is  oriented  toward  this  type  of  experiment.  However,  it  may  be 
noted  that  the  seeding  of  Hurricane  Ginger  (1971)  discussed  below  was  a 
so-called  rainsector  experiment:  all  the  water-containing  clouds  in  a 
45°  sector  from  50  n  mi  to  100  n  mi  from  the  center  are  seeded  at  about 
22,000  ft  to  divert  low-level  inflow  to  premature  eruption  and  high-level 
outflow  outside  the  eyewall,  with  resultant  dispersal  of  latent-heat 

energy  over  a  wider  region  of  the  storm  (Hawkins,  Bergman,  and  Gentry, 

1972).  In  concept,  a  rainsector  experiment  consists  of  four  fifty-minute 
seedings,  at  fifty-minute  intervals.  The  reason  a  rainsector,  rather  than 
eyemod,  experiment  was  attempted  on  Ginger  was  that  it  was  a  large  diffuse 
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storm  (moderate  gale  winds  out  to  250  n  mi,  circulation  out  to  400  n  mi) 
with  only  an  ill -defined  eyewall  (maximum  winds  of  about  60-70  kts  at 
40  n  mi  to  60  n  mi  from  the  center). 

12 

Dr.  R.  Cecil  Gentry,  director  of  the  National  Hurricane  Research  Laboratory 
and  of  Project  Stormfury,  stated  recently  that  budgetary  reductions  had 
grounded  the  four  aircraft  used  in  cloud-seeding  experiments,  and  probably 
no  full-scale  attempts  at  hurricane  modification  would  be  undertaken  in  the 
next  three  years  ("U.S.  Hurricane  Control  Test  Program  Cut,"  Los  Angeles 
Times,  4  February  1973,  part  8,  p.  5). 

13 

Calculation  that  undilute  lifting  of  the  autumnal  Caribbean  ambient 
produces  reduced  temperatures  at  all  tropospheric  levels  was  also  carried 
out  by  Riehl  [Riehl,  H.  (1963) .  Science  141,  1001]. 

14 

Consistent  with  his  concept  that  it  is  the  compressional  heating  of  the 
compensatory  downdraft  that  warms  the  air,  because  heat  directly  released 
by  cumulus  convection  is  almost  entirely  expended  in  raising  the  buoyant 
column.  Gray  (Shea,  1972,  pp.  118-122)  suggests  that  turbulent  diffusion 
and  mixing  conveys  heat  from  the  eye  to  warm  the  eyewall. 

1 5 

From  "Introduction  to  the  Atmosphere"  by  Herbert  Riehl.  Copyright  1972 
by  McGraw-Hill  Book  Company.  Used  with  permission. 
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CAPTIONS  TO  FIGURES 


Fig.  1.  Typical  paths  of  tropical  storms  and  cyclones. 

Fig.  2.  NASA  photograph  of  Hurricane  Gladys  west  of  Naples,  Florida  taken 
from  Apollo  7  on  October  7,  1968.  The  maximum  wind  speed  was 
65  kts  at  this  time,  and  the  tropopause  was  at  54,000  ft. 

Fig.  3.  NOAA  radar  photograph  of  mature  hurricane  revealing  lower-level 
rainbands  outside  the  region  of  overcast  (0/C)  and  upper-level 
cirrus  clouds  of  the  outflow. 

Fig.  4.  Schematic  features  of  the  atmosphere  in  winter.  Latent  hej? 
evaporated  from  tropical  oceans  to  the  lower-leve’  tea  in 
carried  toward  the  ITC7  by  the  sketched  meridional  -*  * •  * 

air  rises  as  the  ITCZ  is  approached,  and  the 
rately  is  released  as  sensible  heat  and  qr»,  rst  . 

energy.  Height  of  the  local  tronopausf  ’s  -  ^fl 
th c  general  circulation  are  less  well  dn,'r‘-d  1 
denotes  tropical  air,  MLA,  rn'dlat i tude  air,  ,  pp’  , 

subtropical  jet  stream,  and  PFJ,  polar  front  jet  *.'»»  , 

permission  from  E.  Palmen  and  C.  W.  Newton,  Atmos  r^rr  1 ecu  la t*on 

Systems ,  Academic  Press,  p.  569.) 

Fig.  5.  Left,  the  equivalent  potential  temperature  p  as  a  function  of 
height  z  and  pressure  p,  from  measurements  near  Barbados  in  the 
Lesser  Antilles  in  July-August,  1968.  Four  characteristic  weather 
types  are  represented:  average  (solid),  suppressed  convection 
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(dotted),  moderately  enhanced  convection  (long  dashes),  and 
strongly  enhanced  convection  (short  dashes).  (From  Harsh, 
Echternacht,  and  Garstang,  1971,  p.  127;  with  permission  of  the 
Americal  Meteorological  Society.)  Right,  soundings  taken  at  Gan 
Island  in  the  Indian  Ocean  (00°41'S;  73°09'E).  The  heavier  curves 
(solid  for  static  temperature,  dashed  for  dew  point)  are  the 
average  for  42  wet  days;  the  lighter  curves  are  the  average  for 
113  dry  days.  Measurements  were  made  in  Jul.ys  1960-1964.  (From 
Johnson,  1969,  p.  122;  with  permission  of  the  Royal  Meteorological 
Society. ) 

Fig.  6.  The  average  total  static  tenperature  (H/cp)  for  the  northern 
Hemisphere  in  winter.  The  curve  marked  A  denotes  the  locus  of 
the  minimum  of  (H/cp);  free  convection  can  most  readily  occur 
below  level  A,  and  only  undilute  ascent  can  continue  much  above 
A;  uridilute  ascent  continues  to  level  B,  where  sea-level  values 
are  recovered.  (By  permission  from  E.  Palmen  and  C.  W.  Newton, 
Atmospheric  Circulation  Systans,  Academic  Press,  p.  574.) 

Fig.  7.  Static  temperature  T,  potential-temperature-like  measure  ( fJ/cD) , 
total  static  temperature  (H/cp),  and  the  total  static  temperature 
for  a  hypothetical  atmosphere  saturated  at  the  actual  local 
temperature  and  pressure  (H/cp)  vs.  height  and  pressure  for  the 
West  Indies  ambient  for  September.  Based  on  data  given  by  Jordan 
(1957).  A  parcel  at  height  z-j  will  become  buoyant  at  zo  where 
(H/Cp)  at  z i  exceeds  (H/cp)  at  z^z-j);  The  ineauality  is  satisfied 
for  air  in  the  planetary  boundary  layer,  and  little  else. 
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The  total  static  temperature  (H/c  )  over  the  ocean  in  the 
equatorial  trough  and  in  the  tropics  near  the  subtropics,  as  a 
function  of  pressure.  Both  typical  cumulus  loci  (cu)  and  also 

cumulonimbus  loci  (cb)  are  noted.  (Based  on,  with  permission, 

* 

E.  Palmen  and  C.  W.  Newton,  Atmospheric  Circulation  Systems, 
Academic  Press,  p.  575.) 

NASA  photograph  taken  in  June,  1965  from  Gemini  IV  of  cumulus, 
cumulonimbus,  and  cirrus  clouds  over  the  Pacific  Ocean  off  the 
western  coast  of  Central  America;  the  view  is  northeast  toward 
Mexico . 

Fig.  10.  This  conjectured  configuration  of  a  mature  hurricane  with  rough 
order-of-magnitude  dimensions  is  not  drawn  to  scale.  The  sub- 
domains  are:  I,  throughput  supply,  a  region  of  rapid  swirl  and 
very  slow  downdraft;  II,  frictional  boundary  layer;  III,  eyewall; 
and  IV,  eye.  Across  the  boundary  layer  there  is  about  a  100  mb 
drop,  and  across  I,  a  further  drop  of  about  200  mb;  the  pressure 
at  the  top  of  the  hurricane  is  about  150  mb,  i.e.,  the  top  is 
near  the  tropopause.  The  eye-eyewall  interface  is  taken  to 
slope  outward,  though  the  effect  may  not  be  so  pronounced  as 
sketched. 

Fig.  11.  Typical  vertical  profiles  of  thermodynamic  state  variables  for 
the  West  Indies  in  September,  based  on  data  by  Jordan  (1957). 

The  relative  humidity  RH,  the  water  vapor  pressure  pv,  the  total 
gas  density  p,  the  water  vapor  density  pv,  and  temperature  T  are 
plotted.  Jordan's  data  for  the  relative  humidity  stop  at  400  mb. 


Fig.  8. 


Fig.  9. 
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Fig.  12.  The  ambient  pressure-temperature  curve  is  based  on  data  for  the 
Caribbean  in  September  given  by  Jordan  (1957).  The  moist  adiabat 
is  based  on  having  sea-level  ambient  air  undergo  dry-adiabatic 
expansion  until  saturation,  then  moist  adiabatic  ascent;  the 
resulting  sea-level  pressure  deficit  from  ambient  qives  a  lower 
bound  on  the  central  pressure  deficit  under  the  adopted  model. 

An  upper  bound  on  the  deficit  is  furnished  by  having  the  air  that 
rose  on  the  so-called  moist  adiabat,  recompressed  drv  adiabatically 
back  down  to  sea  level.  Altitudes  are  associated  with  the  thermo¬ 
dynamic  states  by  use  of  hydrostatics  and  the  equations  of  state 
for  dry  air  and  water  vapor. 

Fig.  13.  The  density  vs.  altitude  relation  for  each  of  the  three  thermo¬ 
dynamic  loci  of  Fig.  12  are  presented.  The  density  variations 
from  ambient  within  the  storm  are  less  than  30%,  and  speeds  are 
far  below  sonic  everywhere. 

Fig.  14.  The  three  curves  are  analogous  to  those  of  Fig.  12,  except  that 
the  moist  adiabat,  while  still  taken  to  be  based  on  ambient  sea- 
level  values  for  relative  humidity  and  temperature,  is  now  com¬ 
puted  From  an  iteratively-determined  sea-level  pressure  consistent 
with  the  definition  of  the  top  of  the  storm  as  an  isobaric, 
isothermal  surface  of  constant  altitude. 

Fig.  15.  In  this  replotting  of  Fig.  14,  explicit  pressure  vs.  altitude 
curves  for  the  ambient,  the  moist-adiabat  taken  to  characterize 
the  eyewall ,  and  the  dry-adiabat  recompression  taken  to  charac¬ 
terize  the  eye  are  presented.  The  moist  adiabat  is  computed  bv 
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the  iterative  procedure  described  in  the  caption  to 
Fig.  14. 

16.  Nondimensional  results  for  the  frictional  boundary  layer  obtained 
by  the  method  of  weighted  residuals,  from  Dergarabedian  and 
Fendell  (1972).  The  divergence  w(x,?  °°)  and  the  volumetric 

r 

flux  6  =  -  I  if)(x,?)  d?  are  presented  for  the  impressed  swirl 
*0 

V  =  1  -  x/x-| ,  x.|  =  20,  believed  pertinent  to  a  hurricane  outside 
the  eyewall.  Except  near  the  axis  where  nonlinear  inertial 
effects  dominate,  the  linear  Ekman  layer  result, 
w (x,?  -*•»)  =  4>  /2^2,  is  an  excellent  approximation  to  the 

A 

numerical  results.  The  volumetric  flux  6 ( x)  is  thus  linearly 
proportional  to  (x-|  -  x)  to  good  approximation.  Normalized 
residuals  indicate  large  errors  for  x  <  3,  and  discount  the  pre¬ 
mature  eruption  as  an  artifact  of  the  method.  The  solution  by 
Carrier  (1971a)  for  small  x  indicates  the  adequacy  of  the  linear 
Ekman  result  for  the  divergence  to  within  a  factor  of  two.  Since 

v  =  (1/75)  mivhr  and  n  =  (1/16)  rad/hr,  dimensionally  the  results 

1/2 

imply  the  boundary  layer  is  of  thickness  0(v/n)  '  =  0(1  mi), 

1  /  n 

the  downflux  into  the  boundary  layer  is  (2vft)  /f'  w(x,^  ->  »)  = 

_3 

0(5  x  10  mi/hr),  and  the  volumetric  flux  erupting  un  the  eye- 
wall  is  (2tt2  yo2  v>/n)^2  s(x)  =  0(7  x  10^  mi3/hr)  where 
characterizes  the  eyewall  relative  angular  momentum  per  unit  mass. 
The  implication  is  that  the  fluid  initially  in  the  boundary  layer 
sustains  the  hurricane  for  about  a  week,  and  the  fluid  lying 
above  the  boundary  layer  (with  supplementary  replenishment  of 
moisture  from  the  ocean)  can  readily  sustain  the  hurricane  for 
more  than  a  week  more. 
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Fig.  17.  Measurements  of  minimum  sea-level  pressure  vs.  local  sea-surface 
temperature  for  several  tropical  cyclones  (from  Sentry,  1969a , 
p.  406;  reproduced  with  permission  of  the  American  Meteorological 
Society.) 

Fig.  18.  Thp  central  sea-level  pressure  is  plotted  as  a  function  of  the 
local  sea-surface  temperature  for  Hurricane  Ginny,  21-28  October 
1963,  which  passed  over  the  Gulf  Stream  on  24  October  (from 
Perlroth,  1967,  p.  266). 

Fig.  19.  In  this  schematic  view  by  Carrier  of  the  flow  configuration 
and  circumferential  velocity  distribution  in  an  intensifying 
tropical  depression  at  some  early  time  t  =  t-j  (say),  the  inter¬ 
face  C-C  between  the  new  and  initial  air  in  the  core  is  idealized 
as  horizontal  for  convenience.  However,  it  is  easy  to  show  that 
for  all  but  very  modest  circulations  (i.e.,  excent  for  angular 
velocities  less  than  three  times  ambient),  a  rigidly  rotating 
core  (which  implies  uniform  updraft  velocity  at  the  top  of  the 
nonlinear  frictional  layer)  cculd  not  accent  all  the  fluid 
pumped  through  the  Ekman  boundary  layer.  Thus  a  uniform  ascent 
in  the  core  is  dynamically  impossible,  and  in  fact  ascent  in  the 
eyewall  is  concentrated  in  cumulonimbi.  For  the  gross  balances 
being  discussed,  many  points  can  be  made  without  accounting  for 
such  refinements  in  eyewall  structure,  though  ultimately  details 
of  the  structure  are  crucial.  (From  Carrier,  1971b,  p.  146.) 

Fig.  20.  Intensification  from  tropical  depression  to  hurricane  has  pro¬ 
gressed  to  a  more  advanced  stage  in  this  schematic  diaqram, 
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holding  at  t  =  tg  >  t-j,  according  to  Carrier's  model.  The  manni- 
tude  of  the  maximum  swirl  is  increased  and  its  position  lies 
closer  to  the  axis  of  symmetry.  (From  Carrier,  1971b,  p.  146.) 

Fig.  21.  Schematic  picture  by  Carrier  of  the  flow  configuration  which 
prevails  when  the  radius  of  maximum  swirl  R  is  increasing  and 
an  incompletely  formed  eye  is  being  filled  with  relatively  dry 
and  motionless  air,  which  sinks  down  from  the  top  of  the  storm 
under  dry-adiabatic  compression.  With  the  formation  of  an  eye 
lighter  in  weight  than  the  eyewall,  the  terminal  stages  of 
intensification  and  the  beginnings  of  auasisteady  mature-stage 
structure  are  realized.  (From  Carrier,  1971b,  o.  150.) 

Fig.  22.  Schematic  diagrams  of  the  seauence  of  thermodynamic  stens  by 
which  the  equation  describing  the  moist-adiabatic,  process  is 
derived.  The  upper  diagram  sketches  operations  on  the  water 
vapor  component,  and  the  lower  diagram  sketches  operations  on  the 
dry  air  component;  the  diagrams  read  from  top  to  bottom.  Arrows 
on  the  left  of  the  diagrams  denote  work  done  on  the  system  (arrow 
in)  and  work  done  by  the  system  (arrow  out).  Arrows  on  the  right 
of  the  diagrams  indicate  heat  removed  from  the  system  (arrow 
out)  and  heat  added  to  the  system  (arrow  in). 
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